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Abstract 

The  Department  of  Energy  is  in  the  process  of  conducting  a  Remedial 
Investigation/Feasibility  Study  for  a  site  contaminated  with  Dense  Non- Aqueous  Phase 
Liquid  (DNAPL)  pollutants  at  their  Paducah  Kentucky  facility.  A  major  part  of  this  effort 
is  the  screening  of  alternative  technologies.  This  thesis  effort  focuses  on  acquiring  insight 
into  a  number  of  remediation  technology  trains  that  are  candidates  for  the  Paducah  site. 
This  insight  is  used  to  recommend  and  justify  the  screening  of  candidate  technology  trains. 
The  research  makes  use  of  two  decision  analysis  models  (one  is  deterministic,  the  other  is 
probabilistic)  built  to  provide  a  quantitative  assessment  of  the  candidate  technology  trains. 
Dominance  considerations  and  multi-attribute  utility  theory  are  utilized  to  make  the 
quantitative  assessments  and  to  gain  insight  into  each  candidate  technology  train. 

The  results  of  the  analysis  provide  the  DOE  with  a  rational  justification  for 
screening  55  of  the  58  candidate  technology  trains  from  further  consideration.  This 
represents  a  95%  reduction  in  the  number  of  technology  trains  that  are  under 
consideration  for  remediating  the  Paducah  site. 
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Analyzing  Remediation  Technologies  for  Department  of 
Energy  Sites  Contaminated  With  DNAPL  Pollutants 


I.  Introduction 


1.1  General  Issue 

One  of  the  most  difficult  environmental  challenges  facing  the  United  States  is  the 
restoration  of  sites  contaminated  with  Dense  Non- Aqueous  Phase  Liquids,  (DNAPLs). 
DNAPLs  are  toxic,  immiscible  fluids  that  have  a  density  greater  than  water.  DNAPLs  are 
typically  segregated  into  four  distinct  categories:  1)  halogenated  hydrocarbon  solvents,  2) 
creosote  based  wood-treating  oils,  3)  Polychlorinated  Biphenyls  (PCBs),  and  4)  pesticides 
[Cohen  and  Mercer,  1993,  2-2], 

The  contaminant  requires  special  attention  because  of  its  physical  properties, 
unique  migration  pathways,  and  toxicity.  Due  to  the  widespread  production,  use  and 
transport  of  these  chemicals  over  the  past  180  years,  it  is  estimated  that  there  are 
thousands  of  sites  contaminated  by  DNAPLs.  These  sites  range  from  coal  tar 
contaminants  that  were  generated  during  the  mid  1800’s  to  degreasing  agents  used  in 
modern  manufacturing  facilities  [Cohen  and  Mercer,  1993,  3-1].  Cohen  and  Mercer  note 
that  in  1988  there  were  55  wood  preserving  sites  listed  on  the  National  Priorities  List 
(NPL).  These  sites  are  contaminated  with  creosote  and/or  coal  tar.  Additionally,  by 
1990,  17%  of  the  NPL  sites  were  classified  as  PCB  sites.  These  sites  (and  many  others 
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yet  to  be  detected )  contain  an  estimated  290  million  pounds  of  PCBs  that  have  not  been 
destroyed  by  incineration  [Cohen  and  Mercer,  1993,  3-19], 

At  some  sites,  DNAPL  presence  has  been  immediately  detected.  At  others, 
because  of  the  limited  volume  of  DNAPL,  the  source  and  amount  of  the  DNAPL 
contaminant  is  difficult  to  determine.  As  a  result,  the  presence  of  DNAPL  is  often 
detected  at  sites  where  DNAPL  contamination  was  not  initially  suspected  [Cohen  and 
Mercer,  1993,  3-1].  Such  a  discovery  can  compromise  a  remediation  process  that  did  not 
anticipate  the  presence  of  DNAPLs. 

At  present,  there  are  fifteen  confirmed  or  suspected  DNAPL  sites  at  various 
Department  of  Energy  (DOE)  complexes  across  the  United  States  [DOE  Technology 
Summary,  August  1996;  55]. 

To  address  the  contamination  at  these  sites,  a  number  of  competing  technologies 
are  being  developed  for  remediation  and/or  containment.  Each  of  these  technologies  is  in 
varying  stages  of  development,  and  must  be  objectively  evaluated  for  application  at  known 
and  projected  DNAPL  sites.  For  this  reason,  the  DOE  Subsurface  Contaminants  Focus 
Area  has  requested  an  analysis  that  will  provide  insight  into  the  selection  of  candidate 
technology  trains  that  may  be  used  to  remediate  DNAPL  sites.  In  order  to  provide  this 
analysis,  the  decision  situation  surrounding  the  selection  of  technologies  must  be 
appropriately  modeled  and  analyzed. 

This  thesis  effort  focuses  on  the  development  of  two  DA  models  that  provide  the 
decision  makers  with  insight  into  the  complex  process  of  selecting  remediation  trains.  The 
two  models  built  were  used  to  1)  screen  candidate  trains,  and  2)  gain  insight  into  the  risk 
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aspects  of  the  reduced  set  of  screened  trains.  The  results  of  these  models  will  provide  the 
DOE  with  a  systematic  analysis  and  insight  into  a  number  of  competing  remediation 
technology  trains.  The  demonstration  analysis  will  focus  on  the  WAG6  site  in  the  C-400 
area  at  the  Paducah,  Kentucky  DOE  complex. 

1.2  Background 

The  Subsurface  Contaminants  Focus  Area  is  tasked  with  developing  technologies 
that  will  characterize,  remediate  and/or  contain,  and  monitor  contaminated  subsurface 
sites.  This  charge  includes  DNAPL  sites.  The  term  “remediate”  is  used  to  indicate  that 
the  contaminated  site  will  be  cleaned  up  or  controlled  to  a  level  that  will  provide  overall 
protection  of  human  health  and  the  environment  as  required  by  the  Comprehensive 
Environmental  Response,  Compensation,  and  Liability  Act  (CERCLA)  [DOE/OR/07- 
1243&D3;  9-3].  The  term  “contain”  is  used  to  describe  any  process  that  is  used  to 
prevent  the  contaminant  from  spreading.  A  typical  containment  technique  is  to  pour  a 
concrete  barrier  around  the  contaminated  site,  creating  a  physical  barrier  that  the 
contaminant  cannot  penetrate.  The  term  “technology”  is  used  to  denote  techniques, 
methods,  or  processes  that  may  be  useful  for  characterizing,  remediating,  containing, 
and/or  monitoring  a  contaminated  site. 

Different  technologies  are  used  during  each  stage  of  the  cleanup  process  of  a 
DNAPL  site.  For  example,  Electro-osmosis  is  a  remediation  technology  sometimes  used 
to  remove  DNAPL  contaminants  from  low  permeability  media.  This  technology  uses 
electrodes  installed  into  the  ground.  An  electrical  charge  flowing  between  electrodes 
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causes  contaminant  migration.  Once  the  contaminant  has  migrated  in  sufficient  quantities, 
it  may  be  either  adsorbed  or  destroyed  in-situ  [DOE  Technology  Summary,  August  1996: 

86], 

A  “technology  train”,  or  simply  a  “train”,  is  used  here  to  identify  a  number  of 
technologies  that  are  grouped  together  to  form  a  complete  technology  package  that  will 
remediate  or  contain  a  site.  The  tactics  that  make  up  a  specific  technology  tram  are  the 

selected  containment,  removal,  or  remediation  technologies. 

At  present,  the  Subsurface  Contaminants  Focus  Area  has  no  standardized  means 
for  selecting  strategies  to  cleanup  DNAPL  contaminants  at  specific  sites.  A  great  deal  of 
uncertainty  surrounds  these  decisions.  It  is  not  certain  which  developing  technologies  will 
succeed;  which  will  fail  (i.e.  a  technology  may  not  work  as  anticipated  and  the  time  and 
money  invested  in  that  particular  effort  will  have  been  wasted);  or  the  true  cost  and/or 
time  each  emerging  technology  will  require  to  become  fully  operational  at  a  specific 
contaminated  site.  Additionally,  there  is  a  substantial  uncertainty  related  to  the  sub¬ 
surface  geology.  Soil  permeability,  geologic  formations,  and  contaminant  quantities  and 
densities  are  a  few  of  the  sources  of  uncertainties  that  must  be  taken  into  account  when 
deciding  upon  a  strategy  for  a  specific  site.  These  uncertainties,  combined  with  the 
combinatorial  nature  of  selecting  sets  of  technologies  to  makeup  a  train  for  a  particular 
site,  make  this  decision  quite  complex. 
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1.3  Problem  Statement 


The  DOE  requires  a  decision  support  system  which  incorporates  operational  risk 
for  assessing  remediation  technology  strategies  at  DNAPL  sites.  The  decision  support 
system  will  be  applied  at  the  WAG6  site  in  Paducah,  Kentucky,  to  assist  the  DOE  in  the 
screening  process  of  remediation  trains.  Other  risk  elements  such  as  health  and  safety,  or 
environmental,  have  not  been  incorporated  in  this  modeling  effort  at  the  request  of  this 
project’s  sponsor  due  to  previous  technology  screens  conducted  for  the  WAG6  site. 

1.4  Research  Objective 

The  primary  objective  of  this  research  is  to  provide  the  DOE  with  insight  into  the 
candidate  technology  trains  that  are  being  considered  for  remediating  sites  contaminated 
with  DNAPL  pollutants.  The  approach  is  to  build  the  decision  support  system  using  a 
screening  model  and  a  probabilistic  model,  with  interaction  from  key  DOE  operational  and 
decision  making  personnel.  The  screening  model  will  be  used  to  analyze  the  decision 
situation  with  a  focus  on  screening  the  candidate  trains  to  a  set  of  top  contenders.  The 
probabilistic  model  will  be  used  to  assess  the  top  contenders  with  respect  to  uncertainty 
and  risks.  The  models  are  demonstrated  by  analyzing  the  WAG6  site  for  a  specified  spill 
volume.  The  main  value  of  the  models  is  to  provide  insight  that  is  useful  in  justifying  the 
screening  of  technology  trains.  This  insight  is  gained  through  the  use  of  utility  theory, 
dominance  considerations,  and  risk  and  sensitivity  analysis.  It  will  assist  the  DOE  decision 
makers  in  their  quest  to  identify  and  select  a  “best”  technology  train. 
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1.5  Scope  of  Thesis 

The  scope  of  this  effort  will  be  restricted  to  the  following  areas: 

.  A  comprehensive  literature  review  to  investigate  the  following  topics: 

.  Decision  analysis  -  with  specific  attention  to  geologic  applications 
.  Multiattribute  Utility  Theory  (MAUT) 

.  The  development  of  DA  tools  that  will  accommodate  life  cycle  cost,  time,  stochastic 
performance  measures,  and  decision  maker  preferences  for  selecting  remediation  or 
containment  technology  trains  for  DNAPL  contaminated  sites 
.  A  thorough  sensitivity  analysis  of  the  results  of  the  WAG6  site  specific  analysis  for  a 
specified  spill  volume 

.  Conclusions  regarding  the  model  and  the  sensitivity  analysis 
.  Recommendations  for  further  research 

1.6  Overview 

Chapter  2  provides  a  discussion  of  the  current  literature  relating  to  this  topic. 
Chapter  3  contains  a  discussion  of  the  methodology  used  to  develop  the  DA  models.  The 
development  of  the  methodology  is  based  on  the  information  contained  in  Chapter  2. 
Chapter  4  discusses  the  output  of  the  models  and  contains  a  thorough  sensitivity  analysis 
of  the  candidate  trains  and  their  key  variables.  Chapter  5  consists  of  the  conclusions  and 
recommendations  that  can  be  drawn  from  the  analysis.  Additionally,  a  detailed  appendix  is 
included  to  aid  in  the  full  understanding  of  the  DA  models  and  the  nuances  of  their 
structure. 
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II.  Literature  Review 


2.1  Introduction 

The  purpose  of  this  literature  review  is  two  fold:  1)  To  provide  the  reader  with  an 
understanding  and  appreciation  of  the  subject  areas  surrounding  the  selection  of 
remediation  trains  for  sites  contaminated  with  DNAPL  pollutants  and  2)  To  investigate 
and  report  on  other  efforts  that  are  similar  to  the  one  undertaken  here.  As  such,  an 
overview  of  DA  is  first  presented.  This  is  followed  with  a  discussion  of  Multiattribute 
Utility  Theory  (MAUT).  Because  of  the  similarities  between  the  gas  and  oil  industry  and 
environmental  restoration/remediation  of  DNAPLs  (both  are  trying  to  remove  specific 
materials  from  the  sub-surface)  a  discussion  of  these  two  areas  and  their  use  of 
DA/MAUT  is  also  presented  . 

2.2  Decision  Analysis 

Decision  analysis,  according  to  Clemen  (page  6),  is  an  iterative  process  used  to 
provide  insight  into  a  decision  situation,  its  pertinent  objectives,  and  its  potential 
alternatives.  The  process  begins  with  a  clear  identification  and  understanding  of  the 
decision  situation.  Once  this  identification  is  made,  the  second  step  is  to  determine  the 
objectives  of  the  decision(s)  and  the  values  of  the  Decision  Maker  (DM).  To  maximize 
the  benefits  of  the  DA  process,  the  objectives  and  values  must  be  stated  as  precisely  as 
possible. 

The  next  step  in  the  process  is  to  identify  the  potential  alternatives  [Clemen,  1996: 
6].  Typically,  the  identification  of  alternatives  makes  use  of  the  previously  stated 
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objectives.  There  are  a  number  of  methods  that  have  been  used  in  the  past  to  create  these 
alternatives;  the  three  most  popular  are  discussed  here. 

Keeney  suggests  an  iterative  process  that  selects  one  objective  at  a  time  while 
ignoring  the  rest,  and  inventing  alternatives  or  scenarios  that  make  that  objective  look  as 
good  as  possible  [Keeney,  1992:  210]. 

Another  possibility,  also  suggested  by  Keeney,  is  to  use  an  approach  that  puts  all 
of  the  stated  objectives  at  their  most  favorable  level.  This  leads  to  the  ideal  alternative. 
This  alternative  can  then  be  analyzed  to  determine  which  constraints  are  holding  one  back 
from  this  ideal  alternative  [Keeney,  1992:  221]. 

Howard  suggests  the  construction  of  a  strategy  generation  table  to  assist  in  the 
identification  of  alternatives.  The  idea  behind  the  strategy  generation  table  is  to  list 
alternatives  with  a  specific  objective  in  a  column  below  that  objective.  Once  each 
objective  has  an  exhaustive  list  of  alternatives,  a  table  is  formed  with  the  columns 
representing  objectives,  and  the  rows  representing  alternatives  for  each  objective.  After 
the  table  is  constructed,  strategies  can  be  selected  and  analyzed  in  an  easily  understood 
and  visual  manner  [Howard,  1988:684]. 

Other  ideas  suggested  by  Clemen  for  generating  alternatives  include  idea 
checklists,  morphological  forced  connections,  means  objective  networks,  brainstorming, 
and  metaphorical  thinking  [Clemen,  1996:  200-207]. 

The  next  step  in  the  DA  process,  as  outlined  by  Clemen,  is  to  decompose  and 
model  the  problem,  decisions,  uncertainties,  and  preferences.  This  step  is  critical  and  must 
be  conducted  with  the  involvement  of  the  DM.  The  primary  decomposition  and  modeling 
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tools  used  are  influence  diagrams  and  decision  trees.  Both  of  these  tools  provide  visual 
insight  into  the  decision(s)  and  problem(s)  at  hand. 

An  Influence  Diagram  provides  a  relatively  simple  way  to  visualize  a  decision 
situation  [Clemen,  1996:  50],  The  Influence  Diagram  makes  use  of  three  simple  block 
shapes  to  represent  specific  pieces  of  the  decision  situation.  Figure  2. 1  summarizes  the 
three  shapes  and  their  significance. 


Decision 

Node 


I  Uncertainty  i 
Node 


Calculation 
or  Value  Node 


The  decision  node  indicates  that  a  decision 
must  be  made:  i.e.  which  alternative,  or 
Yes/No. 

The  uncertainty  node  indicates  that  the 
outcome  of  a  specific  event  is  not  known  with 
certainty. 

The  Calculation  or  Value  node  is  used  to 
perform  calculations  or  to  identify  constants 
within  a  decision. 


Figure  2. 1  Influence  Diagram  nodes 


An  example  of  a  decision  situation  is  shown  in  Figure  2.2.  In  this  figure,  the  decision  is 
which  of  two  remediation  technology  alternatives  to  select.  The  uncertainty  indicates  that 
the  performance  of  the  technology  is  an  uncertain  entity.  The  calculation  or  value  node  is 
used  to  calculate  the  expected  level  of  cleanliness  after  the  site  has  been  remediated. 
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Figure  2.2  Decision  Situation  for  a  Simple  Technology  Selection 


To  expose  the  details  of  the  decision  situation,  a  decision  tree  is  used.  A  decision 
tree  shows  the  mutually  exclusive  and  collectively  exhaustive  outcomes  that  are  possible 
[Holloway,  1979:33].  An  example  of  a  decision  tree  representation  of  the  Influence 
Diagram  shown  in  Figure  2.2  is  shown  in  Figure  2.3. 


Performance 

of 

Technology 

Which 


Figure  2.3  Decision  Tree  Representation 
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In  the  decision  tree  note  that  the  square  represents  the  decision,  and  the  circle  represents 
the  uncertain  node.  The  values  or  outcomes  are  listed  at  the  tail  of  each  branch  of  the 
tree. 

The  final  steps  in  the  DA  process  are  to  select  the  best  alternative,  and  conduct 
sensitivity  analysis  to  assure  that  the  best  alternative  remains  the  preferred  choice  in  the 
presence  of  potential  deviations.  Choosing  the  best  alternative  can  sometimes  be  quite 
difficult,  even  when  an  obvious  best  choice  surfaces  based  on  an  Expected  Value  (EV). 
Risk  profiles,  deterministic  and  probabilistic  dominance,  and  sensitivity  analysis  tend  to 
uncover  subtleties  related  to  the  decision  situation  that  may  not  be  obvious  by  simply 
looking  at  EV  [Clemen,  1996  123]. 

All  of  the  steps  described  in  the  DA  process  are  iterative,  often  requiring  review  at 
each  step.  Decision  Maker  attitudes  and  insights  will  change  as  a  decision  situation  is 
modeled  and  the  fog  of  the  uncertainties  surrounding  a  decision  is  reduced.  It  is, 
therefore,  imperative  that  the  overall  goal  of  the  decision  situation  be  kept  in  view  as  the 
modeling  and  analysis  process  takes  place  [Clemen,  1996:  7]. 


2.3  Multiattribute  Utility  Theory 

Multiattribute  Utility  Theory  (MAUT)  encompasses  a  broad  range  of  topics  that 
are  aimed  at  gaining  insight  into  decisions  that  contain  many,  and  often  times  conflicting, 
objectives  [Clemen,  1996:  532].  To  understand  and  discuss  MAUT,  it  is  important  to 
understand  the  relationship  between  objectives  and  attributes.  Keeney  defines  an  attribute 
of  an  objective  as  “the  degree  to  which  an  objective  is  achieved”  [Keeney,  1992:  100].  He 
goes  on  to  state  that  attributes  “should  be  measurable,  operational,  and  understandable” 
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[Keeney,  1992:  112].  An  example  of  the  relationship  between  attributes  and  objectives  is 
presented  by  Clemen.  If,  for  example,  one  objective  was  to  minimize  costs,  then  the 

natural  attribute  scale  would  be  dollars  [Clemen,  1996.  532]. 

Another  important  concept  that  must  be  understood  when  discussing  MAUT  is 
that  of  additive  utility  functions  versus  utility  functions  with  interactions.  When  an 
additive  utility  function  is  employed,  it  is  assumed  that  each  attribute  has  its  own  utility 
function  and  that  the  attributes  are  additive  and  independent.  The  additive  utility  function 
is  a  weighted  average  utility  function  of  all  of  the  individual  attributes  for  a  given  problem 
[Clemen,  1996:  537].  For  example,  suppose  that  two  objectives  were  stated  as.  1) 
Minimize  the  cost  of  the  project,  and  2)  Minimize  the  time  necessary  to  complete  the 
project.  Clearly  these  are  two  potentially  conflicting  objectives.  MAUT  would  be  a  good 
candidate  for  an  evaluation  technique  for  gaining  insight  into  these  two  objectives.  The 
attributes  of  these  two  objectives  are:  1)  cost  (dollars)  used  to  complete  the  project,  and 
2)  time  (months)  required  to  complete  the  project.  Each  of  these  objective/attribute 

combinations  possess  their  own  utility  fimction. 

The  additive  utility  function  is  simply  a  weighted  average  of  the  combination  of  the 
two  individual  utility  functions.  To  determine  the  weights  assigned  to  each  piece  of  the 
utility  function,  it  is  necessary  to  ascertain  the  relative  importance  of  each  attribute.  For 
example,  is  cost  twice  as  important  as  time  in  completion  of  this  hypothetical  project?  If 
so,  then  the  weight  for  cost  would  be  0.667,  and  that  for  time  would  be  0.333.  By 
convention,  for  the  additive  utility  function,  the  weights  must  sum  to  1.0  [Clemen,  1996: 

537]. 
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According  to  Clemen,  when  interactions  between  the  attributes  exist  (i.e.  attributes 
are  not  independent),  the  use  of  an  additive  utility  function  is  not  justified.  One  must 
consider  a  utility  function  that  captures  the  interactions.  The  utility  function  shown  in 
(2.1)  is  a  general  form  of  a  utility  function  that  contains  interactions. 

U(x,y)  =  Cl  +  C2Ux(x)  +  CaUvCy)  +  C4Ux(x)UY(y)  (2.1) 

To  use  a  utility  function  of  this  form,  it  is  necessary  to  have  seperability  -  i.e.  the  overall 
utility  function  can  be  broken  into  different  pieces  that,  when  put  together,  represent  all  of 
the  different  attributes  and  their  interactions. 

To  assure  that  seperability  exists,  one  must  check  for  two  independence  levels;  1) 
Preferential  Independence,  and  2)  Utility  Independence.  An  attribute  is  preferentially 
independent  of  other  attributes  if  specific  outcomes  of  one  are  not  dependent  on  the 
others.  Likewise,  an  attribute  is  said  to  possess  utility  independence  if  preferences  for 
uncertain  choices  between  different  levels  of  one  attribute  do  not  depend  on  the  level  of 
the  other  attribute(s)  [Clemen,  1996:  581].  Both  conditions  must  be  satisfied  if  one  is  to 
use  a  decomposable  utility  function  like  the  one  shown  in  (2.1). 

However,  because  of  the  difficulty  in  proving  that  these  levels  of  independence 
hold,  most  applications,  especially  first  cut  ones,  assume  complete  additive  independence, 
and  use  an  additive  utility  function  without  any  serious  adverse  results  [Clemen,  1996;  593 
and  Stewart,  1995:  254],  Following  Stewart  and  Clemen’s  suggestion,  the  utility  function 
used  in  this  thesis,  and  described  in  Chapter  3,  assumes  that  the  attributes  are  independent, 
and  makes  use  of  an  additive  utility  function. 
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2.4  Decision  Analysis  and  MAUT  Applied  to  Oil  and  Gas  Recovery 

The  majority  of  the  DA/MAUT  related  research  in  the  oil  and  gas  industry  has 
been  conducted  with  regards  to  oil  exploration.  As  such,  most  of  the  discussion  in  this 
section  is  devoted  to  this  topic.  While  this  is  not  an  exact  1  to  1  mapping  to  conducting 
DA/MAUT  on  the  selection  of  alternative  remediation  technologies,  there  are  similarities 
that  apply  to  this  research. 

2.4.1  Oil  Wildcatting 

An  example  of  the  classic  use  of  DA/MAUT  in  oil  exploration  is  found  in 
Holloway  (pages  195  -  200).  The  decision  situation  presented  is  that  of  choosing  between 
two  sites  for  exploratory  oil  drilling.  One  sight  is  risky  because  of  the  level  of  uncertainty 
surrounding  the  quantity  of  oil  that  may  exist.  The  other  site  is  not  as  risky  because  it  is 
known,  with  a  fairly  high  probability,  that  it  will  produce  at  least  a  small  amount  of  oil. 
Potential  outcomes  (quantity  of  oil  retrieved  and  thus  profits  realized)  are  dependent  on 
the  amount  of  oil  actually  in  the  ground,  the  presence  of  a  dome  structure,  and  the  cost  of 
drilling  each  well.  Only  one  well  (site)  may  be  drilled.  Holloway  demonstrates  the  use  of 
conditional  probabilities  to  gain  insight  into  a  DA  problem,  and  how  to  calculate  the  value 
of  information  and  its  potential  uses  in  gaining  insight  into  a  DA  problem  [1979;  195  - 
200]. 
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2.4.2  Exploration  opportunities  at  Phillips  Petroleum 

Phillips  Petroleum  Company  has  made  use  of  a  DA  based  software  package  called 
DISCOVERY.  The  package  provided  the  Phillips  company  a  means  to  evaluate  risky 
projects  (oil  exploration  opportunities)  in  a  way  that  reflected  the  risk  attitudes  of  their 
company  executives.  Walls,  Morahan  and  Dyer  (pages  39-56)  report  that  Phillips  needed 
to  quantify  the  following  issues: 

.  Evaluate  risk  not  only  as  a  function  of  the  probability  distribution  of  an  oil  well 
outcome,  but  also  as  a  function  of  the  amount  of  capital  exposed  to  the  chance  of  a 
loss. 

.  Trade-off  the  potential  gains  versus  losses  for  individual  projects  (explorations). 

.  Determine  the  appropriate  level  of  involvement  in  other  investments  for  diversification 
purposes. 

Walls,  Morahan  and  Dyer  generate  profit  contribution  based  on  financial  premiums,  cash 
considerations,  risk  preferences,  and  royalties.  A  decision  tree  was  used  to  visualize  the 
decision  situation.  Once  the  profit  contributions  were  established,  an  exponential  utility 
function  was  employed  to  determine  certainty  equivalents,  for  different  levels  of  risk 
aversion.  Certainty  equivalents  are  defined  as  the  certain  value  (amount)  that  a  decision 
maker  is  willing  to  accept  in  lieu  of  an  uncertain  gamble  [Holloway,  1979:101]. 

Walls  et  al.  then  compared  the  impact  of  varying  levels  of  risk  aversion  to  expected 
value  calculations  for  different  levels  of  diversification  involvement.  The  results  indicated 
that  project  rankings  and  their  diversification  share  were  quite  different  when  using  the 
Certainty  Equivalent  basis  versus  the  Expected  Value  basis.  The  researchers  conclude  by 
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stating  that  the  inclusion  of  risk  attitudes  based  on  exposed  capital  and  diversification 
preferences  can  (and  did  in  the  Philips  case)  have  a  dramatic  impact  on  the  preference  for 
project  selections  [Walls  et  al.,  1995:  39  -  56]. 

2.4.3  Case  Study  of  Integrated  Oil  Field  Development  Strategies 

The  development  strategies  for  an  oil  field  are  quite  complex  and  require  an 
analysis  of  uncertain  variables  such  as  seismic  data,  performance  predictions, 
hydrogeology,  water  saturation ,  and  permeability  [Stripe,  Arisaka,  Duandeau,  1993;  155], 
To  better  understand  these  uncertainties  and  their  effect  on  the  selection  of  an  oil  field 
development  strategy,  Stripe  et  al.  conducted  an  analysis  for  the  Hawk  Field  offshore  fi'om 
Abu  Dhabi,  Saudi  Arabia.  In  the  analysis,  optimistic  and  pessimistic  outcomes  were 
evaluated  for  two  different  reservoirs.  The  uncertainties  incorporated  in  the  model  are 
shown  below: 

.  Seismic  data  -  especially  towards  the  boundaries  of  the  reservoir 
.  Water  saturation  levels 
.  Soil  permeability 
.  Aquifer  support  levels 
.  Soil  porosity 

.  Rock  properties  throughout  the  field 

A  total  of  22  different  oil  recovery  methods  were  investigated.  Most  methods 
were  slight  variations  of  a  standard  method.  For  example,  one  recovery  method  studied 
was  to  use  five-spot  water  flooding  with  and  without  peripheral  water  injection.  Another 
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was  to  use  pattern  gas  injection  with  and  without  peripheral  water  injection.  An  initial 
screen  based  on  estimated  field  performance  was  conducted  and  the  list  of  22  removal 
methods  was  reduced  to  12  full  field  production  scenarios  or  strategies. 

The  researchers  then  conducted  an  incremental  economic  analysis  as  well  as  a 
sensitivity  analysis  of  three  of  the  variables.  The  results  were  presented  as  “best 
strategies”  based  on  cost  effectiveness  (from  the  incremental  economic  analysis)  and  level 
of  associated  risk  (from  the  optimistic  and  pessimistic  assessments).  Sensitivities  were 
presented  as  increases  or  decreases  in  incremental  discounted  income  [Stripe  et  al.,  1993: 
155  -  167].  Unfortunately,  the  specific  details  of  the  model  were  not  presented  in  the 
paper,  and  thus  cannot  be  presented  here.  Nonetheless,  the  use  of  DA  to  support 
technologies  for  subsurface  extraction  is  clearly  demonstrated. 

2.5  Decision  Analysis  and  MAUT  Applied  to  Remediation  Projects 

Remediation  projects,  because  of  their  multiple  and  conflicting  objectives,  have 
been  the  source  of  a  number  of  DA  research  opportunities.  This  section  outlines  the 
research  conducted  in  the  past  that  makes  use  of  DA  and  MAUT  principles. 
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2.5.1  Decision  Support  for  CERCLA  Investigations 

Purucker,  Lyon,  Stewart,  and  Nanstad  provides  an  outstanding  overview  and 
introduction  to  DA  and  its  uses  in  support  of  CERCLA  investigations.  The  study,  issued 
in  September  1994,  was  prepared  for  the  DOE’s  Oak  Ridge  National  Laboratory.  In  it, 
the  authors  introduce  the  concepts  and  potential  applications  of  DA  during  a  CERCLA 
process.  The  goal  of  the  introductory  paper  was  to  enlighten  DOE  personnel  about  a 
technique  that  can  be  used  to  supplement  the  normal  Remedial  Investigation/Feasibility 
Study  (RI/FS)  process  and  the  Data  Quality  Objective  (DQO)  process. 

The  paper  outlines  the  DA  process  discussed  in  paragraph  2.2,  but  uses  CERCLA 
language  to  describe  the  steps.  For  example,  the  step  labeled  “identify  the  potential 
alternatives”  is  grouped  into  a  category  known  as  the  feasibility  study.  Purucker  et  al. 
discuss  objective  functions,  deterministic  and  probabilistic  analysis  as  well  as  subjective 
probability  inputs. 

Purucker  et  al.  report  the  specific  application  of  DA  in  the  area  of  DOE  budgetary 
decisions.  The  Environmental  Restoration  Benefit  Assessment  Matrix  (ERBAM)  is  used 
at  the  Oak  Ridge  Operations  -  Environmental  Restoration  Program  to  support  fiscal 
budgetary  decisions.  The  ERBAM  is  a  risk  based  tool  that  uses  prioritization  schemes  and 
multiattribute  utility  analysis  to  rank  order  specific  environmental  restoration  projects.  It 
is  designed  to  measure  the  reduction  of  risk  for  each  project  implemented.  The  output 
from  the  ERBAM  is  a  net  benefit  score  that  is  based  on  the  objective  function  shown  in 
(2.2). 

AProject  Benefits  =  (WiSibLib)  -  ((l-Wi)SiaLia)  (2-2) 
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where; 


Wi  =  the  weight  of  the  impact  -  determined  through  the  use  of  an  objectives 
hierarchy  and  the  establishment  of  performance  scales 
Sib  ~  the  severity  of  the  impact  before  the  project  in  implemented 
Sia  =  the  severity  of  the  impact  after  the  project  is  implemented 
Lib  =  the  likelihood  of  impact  occurring  before  implementation 
Lia  =  the  likelihood  of  impact  occurring  after  implementation 
The  conclusions  drawn  from  the  paper  are  that  DA  techniques  are  being  used  in  the  DOE, 
and  they  should  be  used  more  -  especially  in  the  Remedial  Investigation/Feasibility  Study 
stages  of  a  remediation  project.  The  example  of  ERBAM  was  merely  to  show  that  DA 
techniques  are  employed  by  the  DOE  to  make  budgetary  decisions,  and  that  the  same 
techniques  can  apply  to  more  technical  stages  of  the  remediation  process  [Purucker,  et.  al, 
1994], 

2.5.2  A  Simplistic  Approach  to  Remediation  Alternatives  Selection 

The  November/December  1995  issue  of  The  National  Environmental  Journal 
contained  a  paper  written  by  Christopher  M.  Timm  that  used  a  simplistic  approach  to 
technology  selection  for  a  specific  remediation  project  at  the  DOE  facility  at  Rocky  Flats 
in  Colorado.  According  to  Timm,  there  are  four  steps  necessary  to  making  a  quality 
decision.  These  steps  are:  1)  carefully  define  the  decision  to  be  made,  2)  generate  a  list  of 
achievable  alternatives,  3)  develop  meaningful  and  reliable  information,  and  4)  apply  the 
“logical”  tools  developed. 
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The  site  to  be  remediated  was  a  series  of  five  wastewater  ponds  that  were  used  to 
store  and  evaporate  low  level  radioactive  process  wastes  (sludge).  To  satisfy  step  one,  the 
decision  objective  was  stated  as:  select  “the  remediation  option  that  would  be  the  most 
realistic  in  terms  of  cost  and  schedule.”  In  addition  to  this  primary  decision  objective, 
other  sub-objectives  were  listed  as  well.  These  included  items  that  were  specific  to 
regulatory  requirements  or  future  operational  desires  for  the  site. 

A  list  of  achievable  alternatives  was  generated,  and  reliable  and  meaningful 
information  obtained  for  each  of  them.  This  information  was  probabilistically  based  and 
contained  data  such  as  the  probability  of  a  specific  technology  accomplishing  the 
remediation  in  a  specified  time  period. 

After  the  information  was  gathered,  the  technology  alternatives  were  screened  for 
technical  feasibility.  Two  candidate  technology  alternatives  were  eliminated  by  this 
process.  At  this  point,  nine  remained.  Next,  each  alternative  was  evaluated  to  see  how 
they  met  previously  set  timeline  schedules.  Once  complete,  each  alternative  was  evaluated 
based  on:  1)  critical  elements  that  could  potentially  lead  to  significant  changes  in  either 
cost  or  schedule,  and  2)  life  cycle  costs  and  their  variability. 

Based  on  these  evaluations,  specific  alternatives  (technologies)  surfaced  as  the 
winners  in  a  specific  evaluation  criteria.  For  example,  the  option  labeled  partial  on-site 
disposal  surfaced  as  the  low  cost  technology,  and  the  option  to  rebuild  the  ponds  was  the 
lowest  risk  (least  probable  of  failing)  alternative.  The  end  result  of  the  analysis  was  a  list 
of  technologies  that  could  satisfy  the  remediation  requirements,  and  were  the  best  in  one 
or  more  of  the  evaluation  criteria.  No  attempt  was  made  to  consolidate  or  assign  weights 
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to  each  criteria.  The  information  was  presented  to  the  final  decision  maker  and  a 
remediation  technology  selection  was  made  [Timm,  1995;  46  -  53], 

2.5.3  Decision  Support  System  for  DOE  Test  Pit  9  at  DfEL 

One  example  of  risk  based  DA  research  applied  to  a  specific  site  was  completed  by 
Ralston  in  his  1996  AFIT  thesis  for  the  Department  of  Energy.  In  this  effort,  Ralston  used 
the  DA  process  outlined  in  Section  2.0  to  develop  a  decision  support  system  to  aid  in  the 
selection  of  remediation  strategies  for  a  test  pit  at  Idaho  National  Environmental 
Laboratory  (INEL).  The  DA  model  developed  incorporates  risk  (probability  of  a 
technology  failure),  time,  and  cost  variables  for  each  candidate  remediation  technology. 
Output  from  the  analysis  revealed  that  of  the  approximately  2,100  possible  remediation 
strategies,  less  than  a  dozen  were  viable  alternatives  requiring  further  consideration  given 
the  time  and  cost  constraints  set  by  the  DOE  [Ralston,  1996]. 

2.5.4  Remedial  Alternative  Selection  at  Savannah  River 

A  similar  research  effort  was  conducted  in  a  collaborative  effort  by  Evans, 

Duffield,  Massman,  Freeze,  Stephenson,  and  Buss.  This  DA  effort  focused  on  a  risk 
based  remedial  alternative  selection  at  the  Savannah  River  Site.  In  their  research,  a  DA 
framework  model  was  built  using  three  other  models  as  drivers:  1)  a  risk/cost  based 
economic  model,  2)  a  hydrogeological  simulation  model,  and  3)  a  hydrogeological 
uncertainty  model.  The  purpose  of  the  research  was  to  demonstrate  a  method  to 
determine  the  low  cost  remediation  alternatives  and  to  shed  light  upon  the  cost  drivers  for 
a  waste  site  at  Savannah  River.  The  waste  site  studied  was  an  unlined  earthen  basin 
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disposal  facility  for  low  level  radioactive  effluent  material  and  other  chemicals  used  in 
separation  facilities  on  site. 

The  objective  function  used  for  the  Evans  et  al.  DA  framework  model  is  shown  in 


(2.3). 


(2.3) 


where:  i 

B(t)  =  the  annual  benefits  -  measured  in  dollars 
C(t)  =  the  annual  costs 

R(t)  =  the  annual  risks  -  determined  by  multiplying  a  probability  by  an  expected 
cost  of  a  failure 

This  is  essentially  a  net  present  value  calculation  that  captures  the  three  variables  described 
above.  In  their  research,  the  annual  benefits  were  given  the  value  of  zero,  and  the 
recommendations  were  based  on  the  costs  and  the  risks.  The  costs  were  subdivided  into 
the  following  categories: 

.  Design  capital  cost 

•  Design  Operations  and  Maintenance  (O&M)  Cost 

.  Total  risk  capital  cost  -  an  account  for  corrective  measures  in  the  event  of  failure 
.  Total  risk  O&M  cost  -  an  account  for  corrective  measures  in  the  event  of  failure 
.  Regulatory  penalty  cost 

The  risk  parameter,  R(t),  is  a  function  of  the  probability  of  failure  and  the  cost  of  a  failure 
in  a  given  year.  This  equation  is  shown  in  (2.4). 
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R(t)  =  P.(t)C.(t) 


(2.4) 


where: 

P.(t)  =  the  probability  of  failure  in  year  t 

C.(t)  =  the  cost  of  failure  in  year  t 

Based  on  the  analysis  of  six  different  remedial  alternatives,  the  results  indicated 
that,  in  general,  O&M  costs  for  each  of  the  alternatives  was  the  significant  driver.  These 
costs  ranged  form  54%  to  84%  of  the  total  cost.  Another  important  result  was  that, 
generally,  the  variable  that  exhibited  the  greatest  impact  on  cost  was  the  total  volume  of 
groundwater  extracted  from  the  site.  The  risk  costs  were  typically  not  significant  when 
compared  to  the  O&M  costs  for  this  particular  site.  The  specific  conclusion  drawn  was 
that,  at  the  Savannah  River  site,  the  use  of  extraction  wells  with  surface  water  disposal  of 
treated  effluent  was  the  lowest  cost  alternative  [Evans,  Dufifield,  Massman,  Freeze, 
Stephenson,  and  Buss,  1994], 

2.5.5  Superfund  Decision  Analysis  in  the  Presence  of  Uncertainty 

Jennings,  Mehta,  and  Mohan  describe  a  DA/MAUT  model  developed  for  a 
Superfimd  site  that  pulls  together  the  criteria  necessary  for  assessing  a  site  of  this 
magnitude.  In  Jennings  et  al,  a  modified  version  of  the  National  Contingency  Plan  (NCP) 
Evaluation  Criteria  (see  Appendix  B)  is  used  to  evaluate  alternatives  for  an  old  Wisconsin 
gravel  pit  that  was  turned  into  a  waste  disposal  site.  The  modified  NCP  criteria  were 
assigned  weights  and  each  was  scored  using  a  detailed  worksheet  that  quantifies  each 
criteria  element.  The  results  of  the  composite  score  for  each  technology  alternative  were 
then  used  in  a  normalized  utility  fimction.  The  normalized  utility  function  was  then 
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evaluated  using  Monte  Carlo  techniques  to  simulate  the  uncertainty  surrounding  each  of 
the  variables  and  their  weights. 

Based  on  an  initial  alternatives  screening  using  short  term  and  long  term 
effectiveness,  as  well  as  implementibiltiy  as  screens,  five  remediation  alternatives  were 
subjected  to  the  remaining  selection  process  proposed  by  Jennings  et  al.  The  output  of  the 
analysis  resulted  in  “rating  histograms”  for  each  technology  alternative.  This  provided  a 
visual  means  to  judge  each  alternative  in  a  stochastic  manner.  The  model  also 
accommodates  a  good  degree  of  sensitivity  analysis  that  allows  for  varying  the  weights, 
the  uncertainty  probability  distribution  functions,  and  the  decision  criteria  score  ranges. 

Jennings  et  al’s  results  indicated  that  an  obvious  remediation  alternative  was  not 
available;  however,  it  was  clear  that  some  alternatives  were  better  than  others.  When  one 
of  the  parameters  was  modified  (or  ranged),  the  rating  histograms  shifted.  This  was 
particularly  evident  when  ranging  the  weight  assignments.  For  example,  the  decision 
criteria  “cost”  was  quite  sensitive  to  its  weight  assignment.  By  ranging  this  weight  (the 
other  weights  ranged  accordingly  because  all  weights  summed  to  1.0),  the  “best” 
alternative  selected  would  be  dramatically  different.  In  other  words,  the  weights  of  the 
decision  criteria  were  the  most  sensitive  parameters  in  their  particular  model  [Jennings  et 
al.  1994:  1132-  1150]. 

2.5.6  Modeling  the  Public  as  a  Stakeholder 

To  capture  all  of  the  concerns  related  to  a  remedial  action,  one  must  consider  the 
attitudes  and  objectives  of  three  different  stakeholders:  1)  the  Site  Owner  (SO),  2)  the 
Regulatory  Agency,  and  3)  the  Public  [Apostolakis  and  Bell,  1995:  1021].  To  address 
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these  three  divergent  groups,  Apostolakis  and  Bell  developed  a  methodology  to  include 
risk  attitudes  and  values  into  an  influence  diagram.  The  risk  attitudes  and  values  are 
described  in  terms  of  objectives.  Their  study  indicated  that  the  SO  and  the  Regulatory 
Agency  are  primarily  concerned  with  health  risk  and  cost.  The  Public,  on  the  other  hand, 
is  concerned  with  the  following  issues: 

.  Control;  The  Public  feels  threatened  by  the  contamination  and  would  like  to  somehow 
control  the  outcome  of  any  cleanup  work. 

.  Credibility:  The  Public  feels  more  or  less  at  risk  based  on  the  actions  of  the  SO.  Thus 
the  SO  credibility  is  a  concern  of  the  Public. 

.  Time  to  Completion;  The  Public  is  concerned  about  the  length  of  time  that  a  cleanup 
will  take. 

.  Cost 

To  demonstrate  the  implementation  of  these  concerns  and  objectives  into  a  remedial  action 
decision  situation,  an  Influence  Diagram  was  constructed  based  on  a  theoretical  site  that 
required  a  selection  between  two  different  pump  and  treat  technologies.  The  Apostolakis 
and  Bell  influence  diagram  is  found  in  Figure  2.4. 
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Contaminant 
Level  (Cp 


Figure  2.4  Influence  Diagram  Developed  by  Apostolakis  and  Bell  (page  1024) 
An  explanation  of  the  nodes  shown  in  Figure  2.4  are  discussed  in  Table  2. 1 . 


Table  2. 1  Node  Explanation  for  Figure  2.4 


Node  Identifier 

RA  Contaminant 
Level  (CL) 

Decision 

SO  Action 

Decision 

Cwo 

Chance 

K 

Chance 

PCCw  at  Well 

Chance 

Pcwt 

Chance 

Time  to  Cleanup 

Chance 

Mcwt 

Chance 

Exposure  Factor 

Chance 

Remedial 

Performance 

Chance 

Hazard 

Reduction 

Chance 

Near  Field  Health 

Calculation  or 

risk 

Value 

Credibility  of  SO 

Calculation  or 
Value 

Control  of 

Calculation  or 

Contamination 

Value 

Cost  of  Cleanup 

Calculation  or 
Value 

Value  Public 

Calculation  or 
Value 

Calculation  or 
Value 


Description 

The  decision  here  is  between  MCL  (Maximum 
Contaminant  Level)  and  ACL  (Alternate 

Contaminant  Level).  It  determines  the  amount  of 
contaminant  left  after  cleanup. 

This  decision  node  represents  the  technology 
decision  made  by  the  SO  between  the  two  different 
Dump  and  treat  technologies. 

The  initial  concentration  at  the  site. 

Hvdraulic  Conductivity  of  the  hydrogeology. 

Predicted  Contaminant  Concentration  in  the  water. 

Predicted  Concentration  at  time  t.  Used  to  gauge 
remediation  progress. 

Time  required  to  cleanup  the  site  to  the  level  decided 
upon  by  the  SO. 

Measured  concentration  of  contaminant  in  the  water 
at  a  given  time. 

Represents  the  exposure  assessment  step  in  the 
health  risk  evaluation.  It  is  an  estimate  of  the 
chemical  uptake  from  potential  pathways. 

Used  to  reflect  the  SO’s  ability  to  achieve  the 
established  goals. _ 

Used  to  measure  progress  by  comparing  the  original 
source  concentration  to  the  current  concentration 
after  some  time  t.  (i.e.  (Cwo  -  MCwt)/Cwo  x  100). 

Represents  the  risk  to  human  health  after  cleanup  is 
completed  to  the  SO  specified  level. 

Represents  the  efforts  of  the  SO  and  is  influenced  by 
the  remedial  performance  and  the  communication  of 
risk  to  the  public. 

Represents  the  SO’s  control  over  the  hazardous 
contaminant.  Influenced  by  the  SO  actions  and  the 
hazard  reduction. 

Total  estimated  present  worth  cost  incurred  by  the 

SO  to  cleanup  the  contaminant  to  the  level  agreed 
upon  in  the  decision  node  Level. 

Final  utility  value  of  remediation  technology  as 
perceived  from  the  public. 

Final  utility  value  from  the  Remediation  Authority 
perspective.  Takes  into  account  the  public  utility. 

Final  RA  Value 


Unfortunately,  Apostolakis  and  Bell  do  not  include  a  real  example  application  of  their 
work.  Nonetheless,  it  is  a  useful  representation  of  modeling  stakeholder  preferences  and 
their  attitudes  towards  risk  [Apostolakis  et  al.,  1995;  1021  -  1026], 

A  number  of  attempts  have  been  made  to  implement  DA  and  MAUT  concepts  into 
the  decision  making  required  to  remove  substances  (either  contaminants  or  fuels)  from  the 
subsurface.  However,  there  is  no  known  specific  application  of  remediating  a  site 
contaminated  with  DNAPL  pollutants.  As  stated  in  Chapter  1,  DNAPLs  present  a  unique 
challenge  because  of  their  properties  and  widespread  use.  The  next  chapter  uses  the  ideas 
discussed  in  this  literature  review  to  build  two  models  for  selecting  DNAPL  remediation 
strategies  that  provide  the  greatest  utility  to  the  decision  maker. 
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in.  Methodology 


3.1  Introduction 

Based  on  the  literature  review  and  the  available  documentation,  a  two-phased 
approach  using  multi-attribute  decision  theory  as  a  foundation  was  utilized  in  this  study. 
The  first  phase  constructs  a  MAUT  model  that  was  used  to  screen  the  candidate 
technology  trains.  The  second  phase  uses  the  results  of  the  first  phase  to  build  a  model  to 
conduct  in-depth  risk  analysis  of  the  most  favorable  strategies.  Both  phases  and  their 
results  can  be  used  by  the  DOE  to  complement  other  efforts  that  are  required  during  the 
Feasibility  Study  phase  of  the  C-400  remediation  project. 

This  section  will  cover  the  following  areas; 

.  DOE  requirements  and  model  focus  discussion 
.  WAG6  and  C-400  site  characteristics 

.  The  Screening  (Phase  1)  model  development  and  dominance  discussion 
.  The  Probabilistic  (Phase  2)  model  development 
Due  to  the  iterative  nature  of  modeling  complex  decision  situations,  the  requirements, 
technology  candidates,  and  models  developed  were  continually  modified  and  updated. 

The  major  updates  and  changes  have  been  captured  and  are  discussed  in  Appendix  G. 

3.2  Focusing  the  Model 

At  the  outset  of  this  effort,  the  DOE  expressed  interest  in  the  development  of  a 
generic  decision  model  that  could  be  used  at  any  of  a  number  of  DNAPL  sites.  The 
general  desire  was  to  build  a  modeling  framework  that  uses  hydrogeological 
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characteristics,  cost,  time,  and  performance  data  for  individual  technologies  as  inputs  and 
to  output  technology  train  selections  that  meet  certain  threshold  limits.  At  the  time,  the 
cost  limit  was  stated  as  a  cost  less  than  $50  million  for  the  total  remediation,  with  no  more 
than  $10  million  spent  per  year  during  the  planning  horizon.  The  planning  horizon  time 
constraint  was  6  years.  There  are  no  general  relationships  between  every  possible 
combination  of  cost,  time,  and  performance  for  every  potential  hydrogeologic 
characteristic,  without  prior  knowledge  of  the  site  specific  characteristics.  Because  of 
this,  a  more  site  specific,  tailored  approach  was  used.  The  framework  developed, 
however,  can  be  used  at  other  sites,  provided  it  is  populated  with  appropriate  data  and 
relations. 

Tightening  the  modeling  focus,  two  decision  models  were  developed  for  the  C-400 
site.  The  models  are  explained  in  Section  3.4.  The  requirements  outlined  by  the  DOE 
evolved  as  the  project  progressed.  The  final  set  of  requirements  are  shown  below: 

.  Build  a  decision  model  or  models  that  incorporate  life  cycle  cost,  time,  and 
performance  data 

.  The  models  are  to  use  cost,  time,  and  performance  data  that  has  been  generated  by  a 
life  cycle  cost  model  owned  and  operated  by  Mountain  State  Engineering  (MSE) 

.  The  decision  models  are  to  be  used  to  select  the  best  technology  trains  by  optimizing 
any  one  or  more  of  the  three  attributes  (cost,  time,  performance) 

.  Conduct  sensitivity  analysis  on  the  variables  that  exhibit  the  greatest  degree  of  influence 

•  Include  uncertainty  values  where  appropriate 

It  is  this  final  set  of  requirements  that  the  models  were  built  against. 
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3.3  WAG6  and  C-400  site  characteristics 

This  section  provides  insight  into  the  site  characteristics  for  the  WAG6,  C-400  site 
and  how  they  related  to  the  models. 

The  Paducah  Gaseous  Difiajsion  Plant  (PGDP),  located  in  Paducah,  Kentucky,  is  a 
uranium  enrichment  facility  owned  by  the  DOE  and  operated  under  contract  by  Lockheed- 
Martin  Energy  Systems.  In  May,  1994  the  entire  plant  was  placed  on  the  National 
Priorities  List.  To  address  the  contamination  issues,  a  Remedial  Investigation/Feasibility 
Study  process  was  begun.  During  the  initial  phases  of  this  process,  the  WAG6  site  was 
identified  and  separated  as  an  individual  remediation  project.  The  WAG6  site  is  composed 
of  5  individual  Solid  Waste  Management  Units  (SWMUs).  Four  of  the  five  SWMUs  are 
located  in  an  area  near  the  C-400  building.  The  fifth  SWMU  is  a  transfer  line  located  in 
close  proximity  to  the  others.  These  5  SWMUs  make  up  the  primary  contaminant  release 
points. 

The  two  contaminants  of  concern  at  the  C-400  site  are  Trichloroethylene  (TCE) 
and  Technetium  99  (^^Tc).  Of  these  two,  the  primary  contaminant  of  concern  in  this  study 
is  TCE.  During  the  course  of  this  study,  two  estimates  of  TCE  spill  volume  were 
provided  by  different  DOE  contacts.  One  estimate  ranged  fi-om  2,000  -  50,000  gallons. 
The  other  estimate  ranged  from  100,000  -  500,000  gallons.  Because  of  the  lack  of 
specific  data  and  the  uncertainty  surrounding  the  spill  volume,  the  experts  at  the  DOE 
chose  to  base  the  analysis  on  a  conservative  volume  estimate  of  7,500  gallons.  The  spill 
volume  source  of  7,500  gallons  is  shown  schematically  in  Figure  3.1.  (See  Appendix  J  for 
more  information  regarding  the  spill  volume.) 
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Figure  3.1  Idealized  Schematic  of  Contaminant  Source 
From  Figure  3. 1,  it  can  be  seen  that  there  are  8  different  hydrogeologic  units 
(HUl,  HU2a,  HU2b,  and  so  forth).  Each  hydrogeological  unit  is  delimited  in  Figure  3.1 
by  a  horizontal  dashed  line.  Each  of  these  units  is  composed  of  different  hydrogeological 
characteristics  that  promote  nonhomogeneous  contaminant  migrations.  Based  on 
discussions  with  the  key  decision  makers,  it  was  decided  to  break  the  8  zones  into  3 
operational  zones;  1)  Unsaturated  Zone  (U),  2)  Saturated  Zone  (S),  and  3)  Aquifer  (A). 
This  simplification  was  based  on  the  fact  that  it  is  unlikely  that  more  than  3  different 
contractors,  treating  different  zones,  would  work  on  the  same  site  simultaneously.  The 
logistics  and  coordination  required  to  maintain  more  than  3  contractors  was  considered  to 
be  untenable.  The  use  of  this  simplification  treats  nonhomogeneous  subsurfaces  in  a 
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homogeneous  manner.  Because  of  this,  additional  uncertainty  was  introduced  into  the 
problem.  Nonetheless,  while  this  simplification  aggregates  some  of  the  detailed 
hydrogeological  characteristics,  it  was  felt  by  the  decision  makers  and  the  research  team  to 
be  a  better  reflection  of  the  operational  aspects  of  remediation  projects  that  had  been  done 
in  the  past. 

3.4  Model  Development 

Due  to  the  iterative  nature  of  model  building,  four  different  protot)^e  models  were 
developed  to  varying  degrees  of  maturity.  Each  of  the  models  served  the  purpose  of 
assisting  in  the  screening  process.  Through  the  use  of  these  models,  the  number  of 
potential  technology  trains  was  reduced  from  approximately  16.8  million  to  58.  The 
iterations  required  to  reach  the  requisite  models  are  found  in  Appendix  G.  This  section 
contains  a  complete  description  of  the  requisite  screening  model  (phase  1)  and  the 
probabilistic  model  (phase  2). 

3.4.1  Screening  Model  Development  (Phase  1) 

The  first  model  was  developed  as  a  screening  tool  to  identify  the  most  promising 
trains.  The  decision  situation  is  separated  into  three  distinct  areas;  (1)  Unsaturated, 
Saturated,  and  Aquifer  (USA)  -  a  strategy  that  is  characterized  by  using  trains  that  use  one 
technology  pair  to  treat  all  three  zones,  (2)  Unsaturated/Saturated  and  Aquifer  (U/S  &  A) 

-  a  strategy  that  is  characterized  by  using  trains  that  use  one  technology  pair  to  remediate 
the  unsaturated  zone  and  the  saturated  zone,  and  another  technology  pair  to  remediate  the 
aquifer,  (3)  Unsaturated  and  Saturated/Aquifer  (U  &  S/A)  -  a  strategy  that  is 
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characterized  by  using  trains  that  use  one  technology  pair  to  remediate  the  unsaturated 
zone,  and  another  technology  pair  to  remediate  the  saturated  zone  and  the  aquifer.  Figure 
3.2  shows  the  influence  diagram  of  the  model.  Each  area  represents  a  different  strategy 
that  is  based  on  realistic  trains  developed  through  the  iterative  modeling  process.  The 
trains  are  formed  by  combining  specific  technologies  that  are  used  collectively  to 
remediate  all  3  zones 

The  first  strategy  type,  the  USA  strategies,  are  shown  in  Table  3.2  They  are 
designated  as  USA  trains  because  of  their  ability  to  treat  all  of  the  zones  at  the  C-400  site. 
The  other  two  strategy  types  are:  1)  U/S  &  A  (Unsaturated,  Saturated  &  Aquifer)  and  2) 
U  &  S/A  (Unsaturated  &  Saturated,  Aquifer).  The  U/S  &  A  strategy  is  composed  of  48 
trains  that  are  so  designated  because  of  their  ability  to  treat  the  C-400  site  using  the  same 
technology  pair  for  the  unsaturated  zone  and  the  saturated  zone,  but  require  a  different 
technology  pair  for  the  aquifer.  These  trains  are  shown  in  Table  3.3 .  The  U  &  S/A 
strategy  is  composed  of  8  trains,  and  is  so  designated  because  of  their  use  of  the  same 
technology  pair  for  the  saturated  zone  and  the  aquifer,  but  a  different  technology  pair  for 
the  unsaturated  zone.  These  trains  are  shown  in  Table  3.4.  A  key  for  the  acronyms  used 
in  these  three  tables  is  presented  in  Table  3.5. 
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Figure  3.2  Screening  Model  Influence  Diagram 
An  explanation  of  each  of  the  nodes  shown  in  Figure  3.2  is  found  in  Table  3.1 
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Table  3.1  Node  Explanation  for  Figure  3.2 


Description 

Which  overall 
strategy 

Decision 

The  decision  here  is  between  USA,  U/S  &A,  and  U 
&  S/A  strategy  types 

USA  Trains 

Decision 

The  decision  is  between  2  USA  Trains 

USANPV  Cost 

Value 

Contains  cost  data  for  USA  Trains 

USA  Time 

Value 

Contains  time  data  for  USA  Trains 

USA  Perf 

Value 

Contains  performance  data  for  USA  Trains 

U/S  pairs 

Decision 

The  decision  is  to  select  1  of  4  technology  pairs  to  be 
used  in  the  unsaturated  and  saturated  zones 

Aquifer  pairs 

Decision 

The  decision  is  to  select  1  of  12  technology  pairs  to 
be  used  in  the  aquifer.  It  will  accompany  the  U/S 
pair. 

U/S&  ANPV 
Cost 

Value 

Contains  cost  data  for  U/S  &  A  Trains 

U/S  &  A  Time 

Value 

Contains  time  data  for  U/S  &  A  Trains 

U/S  &  A  Perf  H 

Value 

Contains  performance  data  for  U/S  &  A  Trains 

S/A  pairs 

Decision 

The  decision  is  to  select  1  of  2  technology  pairs  for 

the  saturated  zone  and  the  aquifer 

Unsat  Zone  pairs 

Decision 

The  decision  is  to  select  1  of  4  technology  pairs  to  be 
used  in  the  unsaturated  zone.  It  will  accompany  the 
S/A  pair. 

U&  S/ANPV 
Cost 

Value 

Contains  cost  data  for  U  &  S/A  Trains 

U  &  S/A  Time 

Value 

Contains  time  data  for  U  &  S/A  Trains 

U  &  S/A  Perf 

Value 

Contains  performance  data  for  U  &  S/A  Trains 

Overall  Utility 

Value 

Calculates  the  utility  values  for  each  Trains 

A  technology  pair  is  the  combination  of  a  primary  technology  and  an 
accompaniment  technology.  For  example,  one  technology  pair  shown  in  Table  3.2  is 
DUS,  or  Dynamic  Underground  Stripping  combined  with  Bio  (Biological  Remediation). 
In  this  example,  DUS  is  the  primary  technology,  and  Bio  is  the  technology 
accompaniment.  Together  they  make  up  a  technology  pair  (and  in  this  case  a  tram  too). 
The  other  technology  pair  shown  in  Table  3.2  is  DUS  and  P&T,  or  Pump  and  Treat.  A 
complete  description  of  the  individual  technologies  is  found  in  Appendix  D. 
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Table  3.2  USA  Trains 


for  all  3 
zones 

Number  of 
Combinations 

DUS 

Bio 

P&T 

2 

total  =  2 


Table  3.3  U/S  &  A  (Same  Technology  pair  for  the  Unsaturated 


and  Saturated  Zone,  Different  Technology  for  Aquifer) 


for  unsaturated  & 
saturated  zones 

for  aquifer 

Number  of 
Combinations 

two  phase 

none 

RF 

6  phase 

choose  1  of  12 

3x12=36 

LASAGNA 

none 

choose  1  of  1 2 

1  x12  =  12 

total  =  48 


1  of  12  choice  list 
Oxidation  &  Bio 
Oxidation  &  none 
Cosolvent  &  Bio 
Cosolvent  &  none 
Surfactant  &  Bio 
Surfactant  &  none 
UVB  &  Bio 
UVB  &  None 
P&T  &  Bio 
P&T  &  None 
PTZ  &  Bio 
PTZ  &  None 


Table  3.4  U  &  S/A  (One  Technology  Pair  for  the  Unsaturated  Zone, 
Different  Technology  Pair  for  the  Saturated  Zone  and  Aquifer) 


for  unsaturated  zone 

for  saturated  & 
aquifer  zones 

Number  of 
Combinations 

SVE 

none 

P&T 

Bio 

4 

2  phase 

none 

RF 

6  phase 

none 

4 

total  =  8 
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Table  3.5  Key  for  Tables  3.2,  3.3,  and  3.4 


1  Abbreviation 

Complete  Name 

Heating  technology  that  uses  a  6  point  electrode  scheme 

Bio 

Biological  Treatment 

Cosolvent 

Cosolvent  Flushing 

DUS 

Dvnamic  Underground  Stripping 

LASAGNA™ 

Trade  name  for  specific  heating  technology 

None 

None  -  i.e.  no  accompaniment 

Oxidation 

Makes  use  of  an  oxidant  solution  such  as  H202to  degrade  contaminants 

P&T 

Pump  and  Treat 

PTZ 

Permeable  Treatment  Zone 

RF 

Radio  Frequency  Heating 

Surfactant 

Surface  Active  Agent  Flushing 

SVE 

Soil  Vapor  Extraction 

Two  Phase 

Two  Phase  Heating 

UVB 

UVB 

The  screening  model  shown  in  Figure  3.2  illustrates  the  decision  situation,  but  does 
not  show  the  full  compliment  of  value  nodes  contained  in  the  DPL  model.  Absent  are  the 
value  nodes  that  are  imported  from  an  Excel  spreadsheet.  Because  the  DPL  model  is 
linked  to  an  Excel  spreadsheet,  232  value  nodes  (1  for  each  NPV  cost,  time,  performance, 
and  O&M  cost  value  for  each  train)  are  necessary  to  drive  the  DPL  model.  The  complete 
DPL  model  listing  is  contained  in  Appendix  H. 

The  objective  fimction  contained  in  the  node  “Overall  Utility”  is  composed  of 
weights  and  utility  functions  derived  from,  and  refined  with,  the  Decision  Makers  (DM’s) 
during  a  series  of  meetings.  The  utility  functions  and  weights  were  derived  using  two 
methods.  A  probability  wheel  and  an  Analytical  Hierarchy  Process  technique  were  used 
independently  to  insure  accurate  results.  The  probability  wheel  makes  use  of  an  “eye-ball” 
graphic  technique  that  requires  a  direct  assessment  of  each  attributes  weight.  The 
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Analytical  Hierarchy  technique  involves  the  use  of  specific  trade-offs  between  each 
attribute.  Both  techniques  yielded  similar  results.  The  objective  function  is  shown  in 
(3.1).  The  utility  function  pieces  are  shown  in  Figure  3.3 

Overall  Utility  =  0.44*(cost  utility)  +  0. 12*(time  utility)  +  0.44*(performance  utility)  (3.1) 


1 

S 

0.1 

0 


2000  2006  2012  2018 


Year  Finished  with  Remediation 


Figure  3.3  Utility  Functions 


A  complete  discussion  of  the  results  and  conclusion  derived  from  this  model  are  found  in 
Section  4.2. 


3.4.2  Dominance 

As  part  of  the  screening  process,  an  Excel  spreadsheet  was  built  to  determine  the 
trains  that  were  not  dominated  by  other  trains.  For  this  study,  a  train  is  said  to  be  “strictly 
dominated”  if  each  attribute  value  is  less  than  or  equal  to  that  of  another  train.  For 
example,  if  train  3  had  anNPV  cost  of  $1.4  million,  a  time  value  of  3.5  years,  and  a 
performance  level  of  94%,  and  train  number  144  had  an  NPV  value  of  $2.2  million,  a  time 
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value  of  7.3  years,  and  a  performance  level  of  82%,  then  train  144  is  said  to  be  dominated 
by  train  3.  Since  train  3  is  superior  to  144  in  all  three  of  the  criteria,  there  would  be  no 
need  to  consider  train  144  further  if  the  data  used  for  comparison  were  reliable.  A 
complete  discussion  of  the  results  of  the  dominance  analysis  is  presented  in  Section  4.2. 

3.4.3  Probabilistic  Model  Development  (Phase  2) 

The  probabilistic  model  was  developed  to  gain  in-depth  insight  into  the  most 
promising  candidate  trains  that  were  acceptable  after  the  analysis  based  on  the  screening 
model.  See  Section  4.2  for  an  in-depth  discussion  of  the  results  of  the  screening  model. 
Figure  3.4  shows  the  influence  diagram  of  the  probabilistic  model. 


Figure  3.4  Probabilistic  Model  (Phase  2)  Influence  Diagram 
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An  explanation  of  each  of  the  nodes  shown  in  Figure  3.4  is  found  below  in  Table  3.6. 


Table  3.6  Node  Explanation  for  Figure  3.4 


Node  Identifier 

Description 

Which  Train 

Decision 

Contains  a  list  of  the  seven  trains  chosen  for  further 
analysis.  They  are  no.  3,  4,  94,  119,  135,  148,  and 

156  (156  is  the  baseline). _ _ _ 

Capital  Cost 

Value 

Contains  values  for  the  capital  cost  of  each  train. 

This  data  was  derived  from  the  data  found  in 

Appendix  F. 

O&M  Cost 

Value 

Contains  values  for  the  operating  and  maintenance 
cost  of  each  train.  This  data  was  derived  from  the 
data  found  in  Appendix  F. 

NPV  Cost 

Value 

Combines  the  Capital  Cost  and  the  O&M  Cost.  This 
node  does  not  discount  the  O&M  costs. 

Cost  Utility 

Value 

Converts  NPV  cost  into  a  utility  score  using  the 
functions  found  in  Figure  3.3 

Time  Required 

Uncertainty 

Uses  the  data  found  in  Appendix  I  to  calculate  the 
time  required  for  a  given  performance  level 

Time  Utility 

Value 

Converts  time  used  into  a  utility  score  using  the 
functions  found  in  Figure  3.3 

Interest  Rate 

Value 

Interest  rate  value  used  to  compute  the  present  value 
of  a  stream  of  operating  costs  -  nominal  value;  2.5% 

Performance 

Desired 

Value 

A  single  value  that  can  be  changed  to  analyze  train 
costs  and  time  levels  for  different  performance  levels 

Value 

Converts  the  performance  attained  into  a  utility  score 
using  the  functions  found  in  Figure  3.3 

Overall  Utility 

Value 

Combines  the  individual  utility  scores  into  a 
combined  score  using  (3.1) 

In  addition  to  the  nodes  discussed  above,  there  are  8  other  value  nodes  that  are 
used  to  conduct  sensitivity  analysis.  These  nodes  are  explained  in  Table  3.7. 
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Table  3.7  Node  Explanation  for  Remaining  Nodes 


Node  Identifier 

Description 

Performance 

Utility  Factor 

Value 

Used  to  conduct  sensitivity  analysis  on  the 
performance  utility  function 

Cost  Utility 

Factor 

Value 

Used  to  conduct  sensitivity  analysis  on  the  cost  utility 
function 

import_cost_wt 

Value 

Used  to  conduct  sensitivity  analysis  on  the  cost 
weight  factor  found  in  (3. 1).  This  node  is  imported 
from  an  Excel  spreadsheet 

import_time_ivt 

Value 

Used  to  conduct  sensitivity  analysis  on  the  time 
weight  factor  found  in  (3.1).  This  node  is  imported 
from  an  Excel  spreadsheet 

import_perf_wt 

Value 

Used  to  conduct  sensitivity  analysis  on  the 
performance  weight  factor  found  in  (3.1).  This  node 
is  imported  from  an  Excel  spreadsheet 

export_cost_wt 

Value 

Used  to  conduct  sensitivity  analysis  on  the  cost 
weight  factor  found  in  (3.1).  This  node  is  imported 
from  an  Excel  spreadsheet 

export_time_wt 

Value 

Used  to  conduct  sensitivity  analysis  on  the  time 
weight  factor  found  in  (3.1).  This  node  is  imported 
from  an  Excel  spreadsheet 

export_perf_wt 

Value 

Used  to  conduct  sensitivity  analysis  on  the 
performance  weight  factor  found  in  (3.1).  This  node 
is  imported  from  an  Excel  spreadsheet 

From  Figure  3.4,  it  is  seen  that  one  decision  node  drives  the  model.  This  decision 
node  contains  the  most  promising  trains  that  have  been  selected  by  the  DMs  for  further 
evaluation. 

For  this  in-depth  analysis,  and  to  gain  better  insight  into  the  annual  costs,  the  cost 
figures  need  to  be  separated  into  Capital  and  Operating  and  Maintenance  (O&M).  This  is 
consistent  with  the  work  done  by  Evans  et  al.  In  their  work,  Capital  Costs,  and  Design 
Operations  and  Maintenance  costs  were  broken  out  separately.  Like  the  screening  model, 
these  costs  are  converted  into  utilities  using  the  utility  profiles  presented  in  Figure  3.3. 


3-14 


Note  that  the  Performance  Desired  node  is  not  dependent  on  the  train  decision. 

The  desired  performance  was  modeled  as  an  input  variable  to  determine  the  expected 
value  of  cost  and  time  that  would  be  realized  as  a  result  of  setting  a  specific  performance 
goal.  During  the  analysis,  the  Performance  Desired  variable  was  fixed  for  a  given 
scenario.  It  was  then  manually  changed  to  determine  the  expected  value  of  total  cost  and 
time  required  for  each  train.  For  instance,  by  setting  the  Performance  Desired  value  to  say 
90%,  the  DM  would  be  able  to  determine  the  expected  value  of  total  cost  and  the  time 
required  for  each  train  at  the  90%  performance  level.  This  could  be  compared  to  the 
expected  value  of  cost,  and  time  required  for  95%  and  so  on.  The  performance  value  is 
also  converted  into  a  utility  number  using  the  utility  functions  shown  in  Figure  3.3. 

The  Time  Required  node  is  an  Uncertainty  node  that  represents  the  uncertainty 
surrounding  the  time  that  is  expected  to  be  required  for  a  given  performance  level.  The 
data  was  derived  during  a  technical  discussion  with  DOE  technology  experts.  During  the 
discussion,  three  specific  performance  levels  were  set  to  represent  a  low,  medium,  and 
high  performance  levels  that  might  be  realized  for  each  train.  Then  a  low,  expected,  and 
high  time  value  was  assigned  to  each  performance  level  for  each  train.  The  data  obtained 
during  this  discussion  and  plots  of  this  data  are  presented  in  Appendix  I. 

To  make  the  best  use  of  the  data  that  was  obtained  (see  Appendix  I),  the  Extended 
Pearson  Tukey  Method  was  used.  This  method  was  used  because  of  its  historical  ability 
to  represent  a  multitude  of  distributions  with  only  low,  medium,  and  high  data  points.  The 
method  approximates  continuous  distributions  by  assigning  the  .05  and  the  .95  ffactiles  a 
probability  of  0. 185.  The  median  value  is  assigned  a  probability  of  0.63  [Clemen,  1996: 
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278],  Like  the  other  two  attributes,  the  time  required  is  converted  into  a  utility,  according 
to  the  relationship  given  in  Figure  3.3. 

To  maintain  the  maximum  amount  of  flexibility,  the  model  was  built  so  that  any 
number  of  scenarios  may  be  run  to  maximize  or  minimize  any  of  the  values  listed  below 
for  a  given  Desired  Performance  value: 

•  Capital  Cost  •  O&M  Cost  •  Total  NPV  Cost 

.  Cost  Utility  •  Time  Utility  •  Performance  Utility 

.  Overall  Utility 

By  keeping  track  of  each  of  these  values  during  scenario  runs,  comparisons  and 
specific  insights  were  gained.  These  insights  and  comparisons  are  presented  in  Section  4.3 
The  complete  DPL  model  code  listing  is  contained  in  Appendix  H. 

This  general  methodology,  when  populated  with  site  specific  data  and  relations, 
can  be  used  to  analyze  other  DNAPL  sites. 
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IV.  Results  and  Analysis 


4.1  Introduction 

Analysis  was  conducted  based  on  the  methodology  and  models  described  in 
Chapter  3.  The  results  of  the  analysis  are  presented  here.  Before  presenting  the  results 
though,  it  is  important  to  keep  in  mind  the  purpose  and  objectives  of  this  research.  As 
stated  in  Chapter  1,  a  decision  support  system  which  incorporates  operational  risk  for 
assessing  remediation  technology  strategies  at  DNAPL  sites  was  to  be  built.  This  decision 
support  system,  tailored  for  the  C-400  site,  consists  of  two  models  that  were  developed 
and  used  as  1)  a  screening  model,  and  2)  a  probabilistic  model  which  focuses  on  the  “best” 
technology  trains.  Section  4.2  is  devoted  to  the  first  phase  of  this  analysis  (screening),  and 
Section  4.3  is  devoted  to  the  probabilistic  phase.  Both  models  provide  insight  into  the 
cost,  time,  and  performance  characteristics  of  the  technology  train  candidates.  The 
insights  and  results  of  the  analysis  are  used  to  justify  the  screening  of  technology  trains.  It 
is  this  screening  process  and  the  insight  gained  from  the  results  that  are  the  major 
contributions  of  this  research.  While  the  analysis  is  current  for  the  available  data  and 
assumptions,  care  should  be  exercised  in  drawing  any  final  conclusions  from  the  analysis 
while  questions  of  the  spill  volume  remain  open.  The  analysis  does,  however,  illustrate  the 
type  of  information  that  can  be  made  available  to  the  decision  makers  from  this  project. 

To  fully  understand  the  results  and  analysis,  it  is  necessary  to  review  the 
assumptions  used.  These  assumptions  have  been  discussed  throughout  this  document, 
especially  in  Chapter  3.  The  assumptions  used  are  repeated  below; 
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.  Cost,  Time,  and  Performance  Data  -  The  data  used  as  a  foundation  for  both  models 
(especially  the  screening  model)  was  provided  by  Mountain  State  Engineering 
(MSE).  They  have  built  a  comprehensive  life  cycle  cost  (LCC)  model  that  is 
capable  of  amassing  a  large  volume  of  cost,  time,  and  performance  data,  compiling 
it,  and  providing  output  that  can  be  used  for  analysis  efforts  such  as  this  one.  The 
screening  of  each  train  is  based  on  the  final  data  set  provided  by  MSE  on  15  Jan  97. 
See  Appendix  F  for  a  complete  listing  of  the  data. 

.  Spill  Volume;  -  The  data  set  provided  by  MSE  on  1 5  Jan  97  was  based  on  point 
estimates  for  a  spill  volume  of  7,500  gallons.  See  Appendix  J  for  information  about 
the  uncertainty  of  the  actual  spill  volume. 

•  Probabilistic  Performance  and  Time  Data  -  The  probabilistic  model  (phase  2)  uses 
performance  and  time  relationships  that  were  gathered  on  23  Jan  97  during  a 
meeting  with  DOE  technical  experts.  The  modeling  of  this  relationship  was 
discussed  in  Section  3.4.3.  The  data  is  found  in  Appendix  I. 

.  Cost  Penalty  -  A  cost  penalty  was  assigned  to  trains  that  did  not  meet  a  90% 
performance  level.  The  cost  penalty  was  based  on  an  additional  technology 
application  of  pump  and  treat.  If  a  technology  train  did  not  meet  the  90%  threshold, 
it  was  penalized  by  increasing  its  capital  and  operating  costs  by  assuming  that  a 
pump  and  treat  system  would  be  added  to  the  train. 

.  Cost  Constraint  -  The  Decision  Makers  set  a  cost  constraint  of  $5  million  per  year 
for  O&M  costs. 
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.  Interest  Rate  -  The  Probabilistic  Model  used  a  nominal  interest  rate  of  2.5%  to 
discount  the  O&M  costs.  A  complete  analysis  of  the  sensitivity  of  this  parameter 
indicated  that  the  decision  policies  are  not  sensitive  to  this  choice  of  interest  rate. 

The  remainder  of  this  chapter  is  broken  into  the  following  major  sub-sections: 

.  Section  4.2  -  Results  and  Analysis  of  Screening  Model  (Phase  1) 

.  Section  4.3  -  Results  and  Analysis  of  Probabilistic  Model  (Phase  2) 

4.2  Results  and  Analysis  of  Screening  (Phase  1)  Model 

Because  of  the  deterministic  nature  of  the  data  used  in  the  screening  process,  a 
preliminary  deterministic  based  screening  was  performed.  The  purpose  of  the  modeling 
effort  was  to  screen  technology  trains  from  an  initial  candidate  list  of  58  to  a  smaller  set 
that  would  warrant  more  detailed  analysis.  Table  4. 1  shows  the  ranking  of  the  trains  when 
the  objective  fimction  (3. 1)  is  driven  by  a  single  attribute,  with  the  contributions  of  the 
other  attributes  ignored.  The  reason  for  this  initial  ranking  is  to  determine  which  trains 
exhibited  the  “best”  values  in  each  of  the  attribute  categories.  Column  1  of  Table  4. 1 
gives  the  ranking  with  a  cost  weight  of  1.0;  the  other  weights  (for  time  and  performance) 
are  set  at  0.  Column  2  uses  a  time  weight  of  1 .0;  the  other  weights  (for  cost  and 
performance)  are  set  to  0.  Column  3  of  Table  4. 1  gives  the  ranking  with  a  performance 
weight  of  1.0;  the  other  attribute  weights  (cost  and  time)  are  set  to  0. 
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Table  4.1  Trains  Ranked  Individually  by  Cost,  Time,  and  Performance 


Rank 

Cost 

weight  =  1 

Time 
weight  =  1 

Performance 
weight  =  1 

1 

3 

4 

142 

2 

142 

148 

136 

3 

136 

132 

148 

4 

144 

100 

135 

5 

148 

84 

147 

6 

4 

151 

141 

7 

143 

126 

94 

8 

156 

94 

78 

9 

135 

78 

88 

10 

147 

125 

72 

11 

141 

127 

139 

12 

151 

155 

87 

13 

126 

119 

71 

14 

120 

3 

99 

15 

128 

131 

83 

16 

160 

142 

93 

17 

132 

143 

77 

18 

140 

141 

100 

19 

94 

135 

84 

20 

78 

147 

119 

21 

119 

123 

125 

22 

131 

121 

143 

23 

125 

95 

126 

24 

88 

79 

120 

25 

72 

93 

131 

26 

96 

77 

137 

27 

158 

153 

3 

28 

80 

152 

91 

29 

157 

87 

75 

Rank 

Cost 

weight  =  1 

Time 
weight  =  1 

Performance 
weight  =  1 

30 

139 

71 

132 

31 

100 

156 

140 

32 

84 

128 

4 

33 

127 

160 

95 

34 

155 

96 

79 

35 

95 

80 

89 

36 

79 

124 

73 

37 

87 

122 

123 

38 

71 

92 

138 

39 

99 

76 

92 

40 

83 

90 

76 

41 

93 

74 

127 

42 

77 

99 

121 

43 

137 

83 

144 

44 

138 

120 

90 

45 

153 

139 

74 

46 

152 

137 

124 

47 

124 

144 

153 

48 

123 

140 

152 

49 

122 

138 

96 

50 

91 

136 

80 

51 

75 

91 

122 

52 

92 

75 

128 

53 

76 

89 

155 

54 

121 

73 

151 

55 

89 

158 

158 

56 

73 

157 

157 

57 

90 

88 

160 

58 

74 

72 

156 

From  Table  4. 1,  it  is  seen  that  the  top  5  trains  ranked  by  cost  are  3,  142,  136,  144, 
and  148.  The  top  5  trains  ranked  by  time  are  4,  148,  132,  100,  and  84.  For  performance, 
the  top  5  trains  are:  142,  136,  148,  135,  and  147. 

It  is  interesting  to  note  at  this  point  that,  based  on  the  data  available  and  the  model 
assumptions,  both  of  the  USA  trains  (3  and  4)  are  in  the  top  ranks.  They  both  use 
Dynamic  Underground  Stripping  (DUS)  as  their  primary  technology.  Train  3  couples 
DUS  with  biological  treatment,  and  train  4  couples  DUS  with  pump  and  treat.  The  others 
(84,  100,  132,  135,  136,  142, 144, 147,  and  148)  areU/S  &  A  trains. 

Another  point  of  interest  is  that  of  these  top  U/S  &  A  trains,  3  of  them  (84,  100, 
132)  use  two  phase  vacuum  extraction  as  the  primary  technology  in  the  unsaturated  and 
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saturated  zone.  The  other  6  (135,  136,  142,  144,  147,  148)  use  LASAGNA  as  the 
primary  technology.  In  addition,  note  that  each  of  the  three  two  phase  vacuum  extraction 
trains  use  a  permeable  treatment  zone  (PTZ)  technology  for  the  aquifer.  Conversely,  each 
of  the  6  LASAGNA  trains  uses  a  different  technology  pair  to  treat  the  aquifer.  For 
example,  train  135  uses  LASAGNA  for  the  unsaturated  and  saturated  zone,  and  then  uses 
oxidation  and  biological  for  the  aquifer.  Train  144  uses  LASAGNA  for  the  unsaturated 
and  saturated  zone  as  well,  but  uses  pump  and  treat  for  the  aquifer. 

If  one  were  to  stop  the  analysis  at  this  point,  the  results  thus  far  suggest  that  the 
DUS  (both  USA  trains),  and  LASAGNA  trains  as  well  as  the  trains  that  couple  two  phase 
vacuum  extraction  with  PTZ  should  be  analyzed  further. 

Since  the  decision  makers  are  concerned  with  the  mix  of  benefits  from  these 
attributes,  the  trains  were  ranked  by  overall  utility.  Table  4.2  shows  this  ranking  when  the 
objective  function  is  given  by  (3.1)  -  i.e.  the  weights  for  cost,  time,  and  performance  are 
0.44,  0.12,  and  0.44  respectively.  From  this  ranking,  and  the  use  of  Appendix  F,  it  is  seen 
that  the  top  6  trains  use  LASAGNA  to  remediate  the  unsaturated  and  saturated  zone.  This 
is  another  strong  indication  that,  based  on  the  available  data  and  assumptions,  the 
LASAGNA  trains  are  good  candidates  for  further  evaluation.  The  top  six  trains  were 
chosen  as  a  breaking  point  because  of  the  break  exhibited  in  the  overall  utility  score.  The 
overall  utility  score  drops  from  0.913  on  the  sixth  ranked  train  to  0.899  on  the  seventh. 
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Table  4.2  Trains  Ranked  by  Overall  Utility 


Rank 

Train 

Number 

Overall 

Utility 

Rank 

Train 

Number 

Overall 

Utility 

1 

142 

0.936 

30 

132 

0.597 

2 

148 

0.927 

31 

140 

0.538 

3 

135 

0.917 

32 

4 

0.526 

4 

147 

0.915 

33 

156 

0.522 

5 

141 

0.915 

34 

144 

0.519 

6 

136 

0.913 

35 

151 

0.518 

7 

94 

0.899 

36 

128 

0.510 

8 

78 

0.898 

37 

160 

0.510 

9 

139 

0.888 

38 

96 

0.495 

10 

87 

0.880 

39 

80 

0.495 

11 

71 

0.880 

40 

127 

0.489 

12 

93 

0.878 

41 

155 

0.489 

13 

77 

0.878 

42 

95 

0.487 

14 

99 

0.877 

43 

79 

0,486 

15 

83 

0.876 

44 

153 

0.475 

16 

88 

0.863 

45 

152 

0.474 

17 

72 

0.862 

46 

124 

0.472 

18 

100 

0.802 

47 

158 

0.469 

19 

84 

0.802 

48 

157 

0.468 

20 

119 

0.763 

49 

138 

0.467 

21 

125 

0.761 

50 

123 

0.460 

22 

143 

0.744 

51 

122 

0.458 

23 

126 

0.737 

52 

92 

0.457 

24 

120 

0.728 

53 

76 

0.457 

25 

3 

0.717 

54 

121 

0.446 

26 

131 

0.717 

55 

90 

0.443 

27 

137 

0.697 

56 

74 

0.443 

28 

91 

0.617 

57 

89 

0.427 

29 

75 

0.617 

58 

73 

0.426 

Train  144  (one  of  the  top  ranked  LASAGNA  trains  discussed  previously)  ranked 
no.  34  in  overall  utility  because  it  requires  7.2  years  and  performs  at  a  low  level  of  82%. 
Train  144’s  low  performance  level  is  driven  by  the  use  of  a  pump  and  treat  system  for  the 
aquifer,  whereas  the  other  LASAGNA  trains  use  technologies  for  the  aquifer  that  have 
been  scored  higher  on  performance. 

It  is  interesting  to  note  at  this  point,  that  the  two  USA  trains  (3  and  4)  ranked 
relatively  low  in  overall  utility.  Train  3  ranked  no.  25  and  train  4  ranked  32.  The  reason 
for  this  is  their  relatively  low  performance  level.  Despite  their  low  cost  (train  3  ranked  no. 
1  for  low  cost)  and  their  short  time  horizon  (train  4  ranked  no.  1  for  time),  their  overall 
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utility  score  was  lower  because  of  their  predicted  lackluster  performance.  Train  3  is 
expected  to  attain  a  maximum  performance  level  of  94%,  and  train  4  a  90%  level.  For  a 
complete  listing  and  description  of  these  trains,  see  Appendices  D  and  F. 

To  understand  the  dominance  relationships  between  Net  Present  Value  (NPV) 
cost,  time,  and  performance,  scatter  plots  were  generated  using  the  raw  data  found  in 
Appendix  F.  Scatter  plots  of  cost  versus  time,  cost  versus  performance,  and  time  versus 
performance  were  generated.  They  are  shown  in  the  Figures  4. 1,  4.2,  and  4.3.  The  dotted 
line  represents  the  boundary  of  “best”  values,  the  efficient 
frontier. 


Figure  4.1  NPV  Cost  vs.  Time  Scatter  Plot 
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Performance  I  Performance 


Figure  4.2  NPV  Cost  vs.  Performance  Scatter  Plot 


From  Figures  4.1,  4.2,  and  4.3,  it  can  be  seen  that  trains  3,  4,  142,  and  148  lie  on 
the  boundary  or  frontier  of  at  least  one  of  the  plots.  This  indicates  that  each  of  these 
trains  are  not  dominated  by  the  other  trains.  This  is  a  very  important  finding,  because  if 
we  have  high  confidence  in  the  data  used  to  generate  the  plots,  all  of  the  dominated  trains 
(i.e.  all  trains  except  3,  4,  142,  and  148)  could  be  discarded  from  further  evaluation. 
Another  way  to  visualize  dominance  is  through  the  use  of  a  radar  chart. 


Performance 

(%  contaminant  removed) 


100% 


#3 

#4 

#  142 
#148 

#  121 
dominated 

by  #4  &#  148 


Figure  4.4  Radar  Chart  showing  Dominance  of  Trains 
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Figure  4.4  is  a  graphical  representation  showing  the  nondominated  set  of  trains  and 
one  dominated  train.  The  dominated  train  is  added  to  aid  in  understanding  the  figure. 
When  viewing  a  radar  chart,  one  looks  for  triangles  that  are  completely  encompassed  by 
other  triangles.  These  triangles  represent  trains  that  are  clearly  inferior  to  the 
encompassing  trains;  that  is,  they  are  dominated  by  other  trains.  In  Figure  4.4,  train  121  is 
dominated  by  trains  142  and  148  because  its  triangle  boundaries  are  completely 
encompassed  by  trains  142  and  148.  Asa  result,  this  train  may  be  discarded  from  further 
evaluation  because  its  attributes  of  cost,  time,  and  performance  are  less  desirable  than  the 
same  attributes  for  other  trains  -  in  this  case,  train  142  and  148. 

From  the  results  presented  thus  far,  it  is  clear  that  screening  the  trains  from  58  to  a 
lesser  set  that  exhibit  desired  attributes  is  not  as  straight-forward  as  one  would  hope. 

Based  upon  dominance  considerations  of  the  raw  factors,  trains  3,  4,  142,  and  148  are  the 
clear  choice.  Based  on  utility  considerations  which  recognize  weightings  of  the  factors, 
trains  135,  136,  141,  142,  147,  and  148  are  the  top  contenders.  Table  4.3  shows  the  top 
candidates  based  on  both  dominance  and  overall  utility  scores. 


Table  4.3  Summary  of  Top  Ranked  Trains 


Train 

Notes 

3 

0.717 

Ranked  no.  25  in  utility  and  not  dominated  by  other  trains 

4 

0.526 

Ranked  no.  32  in  utility  and  not  dominated  by  other  trains 

135 

0.917 

Ranked  no.  3  in  utility 

136 

0.913 

Ranked  no.  6  in  utility 

141 

0.915 

Ranked  no.  5  in  utility 

142n 

0.936 

Ranked  no.  1  in  utility  and  not  dominated  by  other  trains 

147 

0.915 

Ranked  no.  4  in  utility 

148 

0.927 

Ranked  no.  2  in  utility  and  not  dominated  by  other  trains 

4-10 


From  Table  4.3,  it  is  seen  that  there  are  2  of  the  4  trains  (#  142  and  #  148)  that  are 
not  dominated  by  other  trains,  and  score  favorably  using  overall  utility  as  a  criteria.  The 
other  two  trains  that  are  not  dominated  (#  3  and  4)  exhibit  relatively  low  utility  scores. 
Their  low  utility  scores  are  due  to  their  relatively  low  performance  value  which  has  a 
weight  of  0.44  in  the  overall  utility  function.  Train  3  performance  level  is  94%.  Train  4 
performance  level  is  90%. 

At  this  point,  it  is  important  to  reflect  on  the  technologies  that  are  represented  by 
the  train  numbers.  Trains  3  and  4  use  Dynamic  Underground  Stripping  as  their  primary 
remediation  technology.  These  trains  are  USA  trains.  The  other  trains  posted  in  Table 
4.3  (135,  136,  141, 142,  147,  and  148)  use  LASAGNA  as  the  primary  technology.  These 
trains  are  U/S  &  A  trains.  From  this  analysis,  it  is  clear  that,  based  on  the  data  provided 
and  the  assumptions  used  to  model  the  decision  situation,  the  trains  shown  in  Table  4.3  are 
the  best  candidates  for  further  evaluation. 

This  initial  set  of  screened  candidates  was  reviewed  by  the  key  decision  makers 
during  a  meeting  at  the  Paducah  DOE  facility.  It  was  decided  to  conduct  further  analysis 
on  trains  3,  4,  94,  1 19,  135,  148,  and  156  (baseline).  Table  4.3a  shows  the  characteristics 
of  each  of  these  trains. 


Table  4.3  a  Summary  of  Trains  Selected  by  Decision  Makers 


Train 

Notes 

3 

0.717 

Ranked  no.  25  in  utility  and  not  dominated  by  other  trains 

4 

0.526 

Ranked  no.  32  in  utility  and  not  dominated  by  other  trains 

94 

0.899 

Ranked  no.  7  in  utility 

119 

0.763 

Ranked  no.  20  in  utility 

135 

0.917 

■gM 

0.927 

mm 

0.522 

Ranked  no.  33  in  utility  1 
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Trains  3  and  4  were  chosen  because  they  use  a  single  technology  pair  to  remediate 
the  entire  site,  reducing  the  complications  involved  in  applying  several  technology  pairs  to 
a  site.  Trains  94  and  1 19  were  chosen  because  they  represent  the  two  phase  vacuum 
extraction  family  of  technologies.  The  DMs  were  interested  in  comparing  these  trains 
with  LASAGNA  trains.  Trains  135  and  148  were  chosen  to  represent  the  LASAGNA 
family  of  technologies.  Train  156  was  chosen  to  represent  a  baseline  train.  It  uses  soil 
vapor  extraction  for  the  unsaturated  zone  and  pump  and  treat  for  the  saturated  zone  and 
the  aquifer. 

4.3  Results  and  Analysis  of  Probabilistic  Model  (Phase  2) 

Based  on  the  screening  model  outputs,  it  was  decided  to  run  the  probabilistic 
model  at  performance  levels  of  80%,  85%,  90%,  95%,  and  99%.  The  data  used  to  drive 
the  calculations  performed  in  the  model  is  found  in  Appendix  I.  The  results  of  these  runs 
are  summarized  numerically  in  Tables  4.4,  4.5,  and  4.6. 


Table  4.4  Time  Required  for  a  Given  Performance  Level 


Performance 

Desired 

Time  Required  (Years) 

Train  Numbers 

3  4  94  119  135  148  156 

80% 

1.1  1.3  13.6  1.2  2.4  42.3  41.2 

85% 

2.3  4.8  17.1  1.9  2.9  51.2  50.4 

90% 

3.4  17.3  21.4  2.8  3.5  60.2  59.9 

95% 

4.5  N/A  N/A  3.6  4.2  69.4  70.0 

99% 

5.4  N/A  N/A  4.3  N/A  77.1  78.5 

\ 
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Table  4.5  Total  Cost  for  a  Given  Performance  Level 


Performance 


Expected  Total  Cost  (Millions  of  Dollars) 
Train  Numbers 


Desired 

3 

4 

94 

119 

135 

148 

156 

80% 

0.6 

2.4 

9.4 

2.4 

2.0 

4.7 

9.1 

85% 

1.0 

3.0 

10.8 

3.2 

2.3 

5.0 

10.0 

90% 

1.3 

4.9 

12.4 

3.9 

2.8 

5.2 

10.7 

95% 

1.6 

N/A 

N/A 

4.8 

3.2 

5.4 

11.4 

99% 

1.9 

N/A 

N/A 

5.5 

N/A 

5.5 

11.8 

Table  4.6  Overall  Utility  for  a  Given  Performance  Level 

Expected  Value  of  Overall  Utility 

Performance 

Train  Numbers 

Desired 

3 

4 

94 

119 

135 

148 

156 

80% 

0.52 

0.36 

0.52 

0.53 

0.37 

0.30 

85% 

la 

0.51 

0.32 

0.51 

0.53 

0.37 

0.29 

90% 

0.54 

0.40 

0.27 

0.50 

0.52 

0.36 

0.28 

95% 

0.75 

N/A 

N/A 

0.71 

0.73 

0.58 

0.49 

99% 

0.92 

N/A 

N/A 

0.87 

N/A 

0.76 

0.66 

Tables  4.4,  4.5,  and  4.6  display  all  of  the  important  values  that  were  extracted 
from  numerous  runs  of  the  probabilistic  model.  Note  that  the  N/ A  entries  denote  that  the 
data  used  to  generate  the  tables  does  not  exist  for  the  specific  conditions.  For  example,  in 
Table  4.4,  in  the  train  4  column  and  95%  performance  row,  N/A  denotes  that  train  4  is  not 
expected  to  reach  a  performance  level  of  95%.  To  gain  better  insight  into  the  data 
presented  in  Tables  4.4,  4.5,  and  4.6  a  number  of  plots  were  developed.  The  plots 
contained  in  Figures  4.5  and  4.6  display  the  data  contained  in  Table  4.4.  The  plots 
presented  in  Figures  4.7  and  4.8  display  the  data  contained  in  Table  4.5.  The  plots  shown 
in  Figures  4.9  and  4. 10  display  the  data  contained  in  Table  4.6. 


443 


Figure  4.5  graphically  shows  the  expected  amount  of  time  that  each  train  will  take 
to  reach  a  given  performance  level. 


Performance  v&  Time 


Desired  Performance  (Percent  Contaminant  Removed) 


- 3 

- 4 

- 94 

- 119 

- 135 

- 148 

. 156 


Figure  4.5  Desired  Performance  vs.  Expected  Time  Required  for  Each  Train 


From  Figure  4.5,  it  is  seen  that  trains  148  and  156  (the  baseline)  are  expected  to  take  a 
much  greater  amount  of  time  to  reach  the  same  performance  levels  as  the  other  trains.  It 
is  also  seen  that  trains  94  and  4  are  expected  to  provide  a  maximum  performance  level  of 
90%.  Train  4  also  has  a  peculiar  knee  in  its  performance/time  relationship.  It  is  expected 
to  reach  80%  performance  in  approximately  1.2  years,  85%  in  5.5  years,  but  is  not 
expected  to  reach  90%  until  approximately  16.3  years.  The  three  trains  bunched  together 
at  the  bottom  of  the  plot  (3,  1 19,  and  135)  are  the  trains  of  most  interest  because  they 
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yield  the  best  performance  in  the  least  amount  of  time.  Figure  4.6  was  developed  to 
further  investigate  these  trains. 


Figure  4.6  Desired  Performance  vs.  Expected 
Time  Required  for  Trains  3,  1 19,  and  135 


From  this  plot  (Figure  4.6),  it  is  seen  that  train  1 19  achieves  the  highest 
performance  level  in  the  least  amount  of  time.  This  is  followed  by  train  3,  and  finally  by 
train  135.  Note,  however,  that  train  135  is  expected  to  reach  a  maximum  of  95% 
performance.  The  others  (3  and  1 19)  reach  their  maximum  at  a  99%  performance  level. 
Of  mild  interest  is  the  point  where  train  3  and  135  cross.  If  one  were  to  choose  between 
these  two  trains  based  solely  on  performance,  train  3  would  be  preferred  unless  it  was 
desired  to  rapidly  reach  performance  in  the  range  of  91%  to  95%.  In  this  range,  train  135 
would  be  the  best  choice  between  the  two.  Above  95%,  train  3  would  again  be  preferred 
over  135  because  135  is  not  expected  to  reach  a  performance  level  beyond  95%. 
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At  this  point  in  the  analysis,  it  is  clear  that  of  the  7  trains  selected  by  the  DMs  for 
fiirther  analysis,  three  have  emerged  as  top  candidates  with  respect  to  the  expected 
amount  of  time  required  for  a  given  performance  level.  The  next  area  investigated  was  the 
cost  of  each  train  at  a  given  performance  level.  Figure  4.7  shows  the  cost  of  each  train  at 
a  given  performance  level.  The  operating  costs  for  each  train  were  discounted  at  2.5%  per 
year. 


Figure  4.7  shows  that  train  94  is  expected  to  be  the  most  costly  due  to  its  operating  costs 
which  are  relatively  high  ($625, 713/year)  and  its  21.4  year  operating  time.  Train  156  (the 
baseline)  also  indicates  very  high  costs  due  to  the  long  operating  time  required  to  reach  a 
desirable  performance  level.  Train  156  takes  approximately  60  years  to  reach  a  level  of 
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90%  performance.  Train  148  also  is  expected  to  take  approximately  60  years  to  reach 
90%  performance,  but  its  operating  costs  are  approximately  1/3  of  train  156.  Train  4 
again  shows  a  distinctive  knee  due  to  its  performance  characteristics  and  the  large  amount 
of  time  required  for  this  train  to  move  from  85%  to  90%  performance.  Trains  3,  1 19,  and 
135  emerge  again  as  the  best  with  respect  to  NPV  costs  and  are  plotted  in  more  detail  in 
Figure  4.8. 


Performance  vs.  Cost  for  Train  3, 119,  and  135 


Figure  4.8  Desired  Performance  vs.  Expected 
Total  Cost  for  Train  3,  1 19,  and  135 

From  Figure  4.8  it  is  clear  that  train  3  is  the  top  performer  with  respect  to  cost, 

135  is  second,  and  1 19  is  third.  Note  that  train  135  is  not  expected  to  perform  better  than 
95%,  whereas  the  other  two  can  potentially  reach  a  99%  performance  level. 

The  final  area  that  was  used  to  investigate  the  seven  trains  was  overall  utility. 
Using  the  utility  functions  and  weights  (cost  weight  =  0.44,  time  weight  =  0. 12, 
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performance  weight  =  0.44)  described  in  Section  3.4.1,  the  overall  utility  was  calculated 
for  each  train.  The  results  are  shown  in  Figure  4.9. 


Performance  v&  Overall  Utility 
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Figure  4.9  Desired  Performance  vs.  Expected  Overall  Utility 


From  this  figure,  it  is  seen  that  the  same  trains  (3,  1 19,  and  135)  are  the  top 
candidates  with  respect  to  overall  utility.  Note  that  train  4  exhibits  its  characteristic  knee 
again.  The  reason  for  this  is  the  same  as  previously  stated.  Note  also  that  the  trains  do 
not  increase  in  overall  utility  until  they  reach  90%  performance.  The  reason  for  this  is  that 
the  contribution  to  overall  utility  from  performance  does  not  “kick  in”  until  a  train  reaches 
the  90%  performance  level. 

To  further  investigate  the  top  candidates,  trains  3,  119,  and  135  are  plotted  in 
more  detail  in  Figure  4. 10. 


85%  90%  95%  99% 


Desired  f^rformance  (Percent  Contaminant  Removed) 
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Performance  v&  Overall  Utility  for  Trains  3, 119,  and  135 


3 

- 119 

-  135 

Figure  4. 10  Desired  Performance  vs.  Expected 
Overall  Utility  for  Train  3,  119,  and  135 

From  this  figure,  it  is  clear  that,  with  respect  to  overall  utility,  the  top  three 
candidates  exhibit  similar  characteristics.  Again,  overall  utility  increases  rapidly  as  the 
performance  increases  past  90%. 

From  the  analysis  thus  far,  it  is  quite  clear  that  there  are  three  trains  that  have 
proven  to  be  the  top  candidates  in  each  of  the  measures.  These  trains  and  the  technologies 
employed  by  each  are  summarized  in  Table  4.7. 


Table  4.7  Top  Three  Trains 


Train 

Unsaturated  Zone 

Aquifer 

3 

USA 

DUS  &  Bio 

DUS  &  Bio 

DUS  &  Bio 

119 

U/S&A 

Oxidation  &  Bio 

135  ^ 

U/S&A 

LASAGNA&None 

LASAGNA  &  None 

Oxidation  &  Bio 
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4.3.1  Sensitivity  Analysis  of  Probabilistic  Model 

A  key  aspect  of  any  analysis  that  can  contribute  to  the  confidence  in  the  overall 
results  is  sensitivity  analysis.  At  this  point,  the  trains  have  been  narrowed  down  to  three 
top  contenders.  The  critical  questions  of  interest  now  are:  Under  what  conditions  will  the 
top  three  trains  remain  the  top  contenders?  If  one  or  more  of  the  parameters  is  varied, 
how  will  the  solution  change?  To  answer  these  questions,  sensitivity  analysis  was 
conducted  on  capital  costs,  O&M  costs,  utility  function  weights,  and  interest  rate  for  the 
top  three  trains  (3,119,  and  135). 

All  of  the  sensitivity  analysis  was  conducted  at  the  93%  performance  level.  This 
level  was  selected  as  a  representative  performance  value  so  that  train  135  would  not  be 
excluded  from  the  analysis.  If  a  value  greater  than  95%  were  chosen,  train  135  would  not 
be  included  in  the  analysis  because  it  can  only  reach  a  maximum  performance  level  of 
95%.  If  a  value  less  than  90%  were  chosen,  the  utility  contribution  from  performance 
would  not  be  included  in  the  analysis  because  of  the  performance  utility  function.  The 
analysis  could  be  conducted  at  any  performance  level  specified  by  the  decision  maker(s). 

Additionally,  all  of  the  sensitivity  analysis  was  conducted  using  overall  utility  as  the 
measure  of  merit.  For  discussion  purposes,  and  comparisons,  this  measure  is  the  most 
straight  forward.  It  is  assumed  that  train  3  is  the  top  candidate,  135  is  the  next  most 
promising  candidate,  and  1 19  is  the  third  ranked  train.  Train  3  has  the  highest  overall 
utility  score,  135  next  highest,  and  119  the  third  highest.  This  is  confirmed  (Figure  4.10) 
by  the  fact  that  the  three  lines  do  not  cross  each  other. 
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The  capital  costs  for  each  of  the  remaining  trains  (3,  1 19,  and  135)  were  separately 
varied  by  ±  50%.  At  the  93%  performance  level,  the  overall  utility  ranking  did  not 
change.  Train  3  held  steady  as  the  best  choice.  This  was  followed  by  train  135  and  then 
train  119.  In  other  words,  the  ranking  of  the  trains  by  overall  utility  at  the  93% 
performance  level  is  not  sensitive  to  changes  in  the  capital  costs  of  the  three  trains.  This  is 
important  because  the  capital  cost  figures  may  have  a  large  margin  of  error.  It  is 
suspected  that  regardless  of  the  performance  level  used  to  conduct  this  sensitivity  analysis, 
the  result  would  be  the  same.  The  main  reason  for  this  is  due  to  the  shape  of  the  cost 
utility  function  which  ranges  from  zero  to  $30  million.  Because  the  costs  used  in  this 
analysis  were  roughly  one  order  of  magnitude  less,  a  change  of  i  50%  in  capital  costs  is 
unlikely  to  aflfect  the  overall  utility  scores  significantly.  A  reduction  of  this  range  may  be 
warranted. 

In  a  similar  manner,  the  O&M  costs  were  varied  by  +  50%  for  the  three  trains  (3, 

1 19,  and  135).  A  policy  change  was  not  realized  as  a  result  of  the  ranging,  however,  there 
were  two  rank  changes  that  resulted  from  ranging  the  O&M  costs  of  trains  1 19  and  135  at 
the  93%  performance  level.  When  the  O&M  costs  for  train  119  was  varied  below 
approximately  $750,000,  train  1 19  ranked  ahead  of  train  135.  Likewise,  when  train  135 
O&M  costs  were  varied  to  greater  than  approximately  $1.1  million,  train  1 19  ranked 
ahead  of  135  in  overall  utility.  Regardless  of  the  ranging,  however,  train  3  always  ranked 
number  1  in  overall  utility.  The  reason  for  the  stability  in  overall  utility  is,  again,  due  to 
the  shape  and  wide  range  of  the  cost  utility  function. 

To  determine  the  sensitivity  of  the  decision  policy  to  the  weights  used  in  the 
overall  utility  function,  each  of  the  weights  were  individually  ranged  from  zero  to  one. 
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Note  that  during  the  ranging  process  the  weights  were  forced  to  always  sum  to  1.0.  The 
results  indicated  that  at  the  93%  performance  level,  the  performance  weight  could  range 
from  zero  to  one  without  instigating  a  policy  change.  The  other  two  weights  (cost,  and 
time)  yielded  a  policy  change  when  they  were  varied.  These  results  are  shown  graphically 
in  Figures  4.11  and  4. 12. 


Sensitivity  to  Cost  Weight 


0  0.1  0.2  0.3  0.4  0.5  0,6  0.7  0,8  0.9  1 

Cost  Weight 


Figure  4.1 1  Sensitivity  to  Cost  Weight 


From  Figure  4. 1 1,  it  is  seen  that,  at  a  93%  performance  level,  the  optimal  policy,  or  best 
train,  changes  when  the  cost  weight  is  set  to  below  approximately  0.05.  At  that  point,  the 
“best”  train,  in  terms  of  overall  utility,  is  1 19.  Otherwise,  train  3  remains  the  best 
selection  in  terms  of  overall  utility. 
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Sensitivity  to  Time  Weight 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0,8  0.9  1 


Time  Weight 

Figure  4.12  Sensitivity  to  Time  Weight 

Figure  4. 12  shows  that  the  best  train  in  terms  of  overall  utility  changes  when  the  time 
weight  is  ranged  above  approximately  0.68.  At  that  point,  train  119  becomes  the  train  of 
choice  in  terms  of  overall  utility  at  the  93%  performance  level. 

The  final  sensitivity  analysis  was  conducted  on  the  interest  rate  used  for  the  NPV 
calculations.  The  interest  rate  was  ranged  from  1%  to  10%  without  causing  a  change  in 
either  the  optimal  decision  policy  (train  3)  or  a  reordering  of  the  top  three  trains.  In  other 
words,  the  overall  utility  score  is  insensitive  to  interest  rate  below  10%.  The  nominal 
interest  rate  used  for  all  of  the  calculations  was  2.5%,  but  any  percentage  between  1  and 
10  could  have  been  used  and  the  relative  results  would  be  the  same. 

In  summary,  at  the  93%  performance  level,  using  overall  utility  as  a  measure  of 
merit,  it  is  evident  that  the  top  three  trains  will  remain  the  top  three  under  a  wide  range  of 
parameter  variations.  It  was  seen  that  the  capital  costs  for  all  three  of  the  trains  could  be 
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varied  by  ±  50%  without  realizing  a  change  in  the  optimal  decision  policy  (tram  3). 
Additionally,  varying  the  O&M  costs  by  +  50%  did  not  yield  a  decision  policy  change, 
however  some  minor  reordering  of  the  top  three  trains  occurred.  It  was  also  shown  that 
by  varying  the  cost  weight  assignment  of  the  utility  function,  the  optimal  decision  policy 
changes  when  it  is  set  below  approximately  0.05.  Likewise,  varying  the  time  weight 
assignment  of  the  utility  function  instigates  a  policy  change  when  it  is  set  above 
approximately  0.68.  Lastly,  it  was  seen  that  the  interest  rate  could  be  ranged  quite 
significantly  without  a  policy  change  or  a  reordering  of  the  top  three  trains.  In  short,  it  is 
quite  evident  that,  at  the  93%  performance  level,  the  top  three  trains  (3,  1 19,  and  135) 
have  been  captured,  and  they  are  quite  robust  to  value  changes.  It  appears  these  three 
trains  remain  quite  robust  over  a  fairly  large  performance  range  as  well,  given  the  initial 

data  and  assumptions. 

4.3.2  Probabilistic  Model  Risk  Discussion 

To  gain  insight  into  the  risks  associated  with  the  selection  of  a  given  train,  it  is 
helpful  to  plot  cumulative  risk  profiles.  Cumulative  risk  profiles  show  the  chances 
associated  with  potential  outcomes.  They  are  used  to  identify  the  relative  riskiness  of 
different  alternatives,  or  in  this  case,  trains  [Clemen,  1996:  119].  One  such  risk  profile  is 
shown  in  Figure  4. 13.  It  is  a  cumulative  risk  profile  based  on  a  93%  performance  level. 
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Figure  4.13  Cumulative  Risk  Profile  for  Overall  Utility  of  Trains 
Number  3,  1 19,  and  135:  Performance  Level  =  93% 


From  Figure  4.13  it  is  seen  that  train  3  stochastically  dominates  the  other  two 
trains  with  respect  to  overall  utility.  In  addition,  train  135  stochastically  dominates  1 19. 
This  indicates  that  at  a  performance  level  of  93%,  if  confidence  was  high  that  the  data  was 
accurate,  a  decision  maker  would  not  select  train  135  or  1 19  over  train  3.  Therefore, 
based  on  the  cumulative  risk  profile,  for  a  93%  performance  level,  and  for  overall  utility  as 
the  measure  of  merit,  it  can  be  concluded  that  trains  119  and  135  can  be  dropped  from 
further  evaluation.  The  remaining  train  (number  3)  stochastically  dominates  the  others. 

Another  insight  can  be  gained  from  the  cumulative  risk  profile.  Focusing  on  train 
3,  in  Figure  4.13,  it  is  seen  that  there  is  approximately  an  18%  chance  that  the  overall 
utility  will  be  0.64  or  less.  In  addition,  there  is  a  81.5%  chance  that  the  overall  utility  will 
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be  approximately  0.67  or  less.  There  is  a  100%  chance  that  the  overall  utility  will  be 
approximately  0.68  or  less.  This  information  provides  the  DM  with  insight  regarding  the 
probabilities  associated  with  a  given  outcome.  Of  course,  these  ranges  are  not  extremely 


wide. 

The  same  analysis  can  be  conducted  for  the  expected  amount  of  time  required  to 
achieve  a  given  performance  level.  One  such  cumulative  risk  profile  is  shown  in  Figure 
4. 14.  This  figure  shows  the  potential  ranges  of  the  amount  of  time  required  to  achieve  a 
performance  level  of  93%  and  the  probabilities  associated  with  each  value. 


Cumulative  Risk  Profile  for  Train  3, 119,  and  135 

Performance  level  =  93% 


Time  Required 

Figure  4. 14  Cumulative  Risk  Profile  for  Time  Required  of  Trains 
Number  3,  1 19,  and  135;  Performance  Level  =  93% 
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From  this  figure,  it  is  seen  that  at  a  performance  level  of  93%,  with  respect  to  the 
time  required  to  achieve  that  performance  level,  train  119  stochastically  dominates  trains  3 
and  135.  These  trains  (3  and  135)  risk  a  higher  number  of  operating  years  to  reach  the 
same  performance  level.  Therefore,  theoretically,  trains  3  and  135  could  be  discarded 
from  further  consideration  if  time  required  were  the  only  criteria  used  as  a  measure  of 
merit  and  one  had  complete  confidence  in  the  data  used  to  generate  Figure  4. 13. 

Figure  4.14  also  shows  the  probabilities  associated  with  specific  time  requirements 
necessary  to  achieve  the  93%  performance  level.  This  information  is  summarized  in  Table 


Table  4.8  Summary  of  Time  Required  and  Associated  Probabilities 


at  93%  Performance  Level 


Probability  Category 

Train  3 

Train  1 19 

Train  135 

18.5%  chance  that  time 
required  will  be  less  than 
or  equal  to: 

1.8 

1.6 

1.8 

81.5%  chance  that  time 
required  will  be  less  than 
or  equal  to: 

3.7 

3.0 

3.5 

100%  chance  that  time 
required  will  be  less  than 
or  equal  to: 

7.5 

5.7 

7.2 

Figures  4.13  and  4.14  are  one  of  many  potential  cumulative  risk  profiles  that  can 
be  generated.  They  were  selected  to  continue  the  discussion  relating  to  93%  performance 
level  that  was  started  at  the  beginning  of  the  Section  4.3.1.  The  expected  time  results  and 
their  associated  probabilities  were  also  determined  for  levels  other  than  93%.  The  results 
of  these  runs  are  shown  in  Table  4.9.  This  table  shows  the  increases  in  times  and  costs 
due  to  desired  performance  level  increases. 
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Table  4.9  Summary  of  Time  Required  and  Associated  Probabilities 
at  Various  Performance  Levels 


Train 

Times  are  in  years 

Costs  are  in  millions  of  dollars 

Performance  Level  I 

80% 

85% 

90% 

95% 

99% 

3 

Exnected  Value  of  time/cos/  required 

1.1/0. 6 

1.311.0 

5.4/7. 9 

18.5%  chance  that  time  required/cosf 
expected  will  be  less  than  or  equal  to: 

0.5/0.4 

1.1/0. 6 

1.5/0. 7 

1.110.9 

2.5/7. 0 

81.5%  chance  that  time  required/cost 
expected  will  be  less  than  or  equal  to; 

0.910.5 

4.111.5 

5. 0/7. 8 

100%  chance  that  time  required/cosf 
expected  viiW  be  less  than  or  equal  to; 

2.5/7. 0 

4.311.6 

6.0/22 

119 

Exnected  Value  of  time/cosf  required 

\. 112.4 

1.913.2 

1.%13.9 

3.614.8 

4.3/5.5 

18.5%  chance  that  time  required/cosf 
expected  vAW  be  less  than  or  equal  to: 

0.5/7.S 

1.0/22 

1.512.6 

1.8/20 

1.013.3 

81.5%  chance  that  time  required/cosf 
cxpecfe</ will  be  less  than  or  equal  to: 

1.0/22 

1.8/20 

1.513.8 

3.314.5 

4.0/5. 7 

100%  chance  that  time  required/cosf 
expected  will  be  less  than  or  equal  to: 

1.512.8 

3.514.7 

5.015.9 

6.517.4 

8.1/29 

135 

Expected  Value  of  time/cosf  required 

1.M2.0 

1.912.3 

3.512.8 

4.112.2 

N/A 

18.5%  chance  that  time  required/cosf 
expected  will  be  less  than  or  equal  to: 

\.lll.2 

1.5/7. 4 

1.111.6 

1.011.7 

N/A 

81.5%  chance  that  time  required/cosf 
expected  wiW  be  less  than  or  equal  to: 

1.111.9 

1.112.2 

3.112.6 

3.8/2.P 

N/A 

100%  chance  that  time  required/cosf 
expected  will  be  less  than  or  equal  to: 

3.913.0 

5.013.8 

6.414.6 

8.0/5.  <5 

N/A 

Note:  Data  convention;  time  value/cosr  value 


Note  that,  in  theory,  a  cumulative  risk  profile  could  be  generated  for  every  conceivable 
performance  level. 
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V.  Summary,  Conclusions,  and  Recommendations 

5.1  Summary  of  Analysis  and  Results 

Based  on  the  results  and  analysis  provided  in  the  previous  sections,  a  clearer 
picture  of  the  decision  situation  has  emerged  for  the  data  and  assumptions  of  the  analysis. 
From  an  initial  deterministic  analysis,  using  utility  and  dominance  theories  as  well  as 
significant  input  from  the  DM’s,  the  number  of  candidate  trains  was  reduced  from 
approximately  13  million  to  7.  Of  these  7  trains,  three  (3,  119,  and  135)  emerged  as  the 
best  in  terms  of  expected  overall  utility,  cost,  and  time  required.  It  was  seen  that  the 
remaining  4  trains  would  take  significantly  more  time  to  reach  desired  performance  levels, 
and  therefore,  because  of  the  resultant  costs,  were  driven  out  of  contention. 

A  sensitivity  analysis  was  conducted  on  capital  and  O&M  costs,  attribute  weights, 
and  the  interest  rate.  This  analysis  indicated  that  the  top  three  trains  are  quite  robust  and 
not  very  sensitive  to  changes  in  any  of  these  parameters  at  the  93%  performance  level. 

By  examining  stochastic  and  deterministic  dominance,  it  was  seen  that  at  the  93% 
performance  level,  train  3  stochastically  dominated  the  others  with  respect  to  overall 
utility.  It  was  also  seen  that  at  the  93%  performance  level,  train  1 19  stochastically 
dominated  the  others  with  respect  to  the  time  required  to  achieve  that  performance  level. 

5.2  Conclusions 

Using  the  data  given,  the  utility  functions  and  weights  derived  from  the  DM’s,  the 
trains  that  emerged  as  the  top  contenders  are  presented  in  rank  order  in  Table  5.1.  The 
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justification  for  this  screening  is  based  on  one  or  more  of  the  following  considerations 
presented  in  Chapter  4; 

.  Dominance  -  Out  of  58  trains,  4  were  not  dominated. 

.  Deterministic  Utility  scores  -  Based  on  the  deterministic  analysis  (phase  1),  6  trains 
were  selected  for  further  evaluation.  The  remainder  exhibited  inferior  utility  scores. 
These  were  combined  with  the  nondominated  set  to  make  up  the  initial  screened 
trains. 

.  Decision  Maker  Selections  -  Based  on  a  presentation  of  the  deterministic  analysis, 
the  key  decision  makers  selected  7  trains  for  further  evaluation.  At  this  point  88% 
of  the  trains  were  screened  from  further  analysis. 

.  Probabilistic  Utility  Scores  -  Based  on  the  probabilistic  modeling  of  time,  given  a 
performance  level,  the  remaining  7  trains  were  further  analyzed.  From  this  analysis, 

3  trains  (3,  119,  135)  emerged  as  top  contenders. 

.  Probabilistic  Performance  Analysis  -  To  further  analyze  the  7  remaining  trains,  an  in- 
depth  analysis  of  the  performance/time  characteristics  was  accomplished. 
Additionally,  performance  ranges  were  investigated.  This  investigation  revealed  the 
potential  time  required  variations  for  the  top  three  trains. 

.  Sensitivity  Analysis  -  To  gain  an  understanding  of  the  sensitivity  of  specific 

parameters  and  their  effects  on  the  outcome  of  train  selection,  a  sensitivity  analysis 
was  conducted  on  the  capital  and  O&M  costs,  attribute  weights,  and  interest  rate.  It 
indicated  that  very  little  sensitivity  existed  in  the  analysis. 
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Table  5.1  Top  Three  Trains 


Train  Number  and  Description 

Discussion 

3  -  This  is  a  USA  Train  which  uses 

Dynamic  Underground  Stripping  and 
Biological  Remediation  Techniques 

Train  with  consistently  high  utility  scores 
across  most  of  the  time  values.  Also  a  train 
that  is  not  dominated  by  any  of  the  other  58 
trains. 

1 19  -  This  is  a  U/S  &  A  Train  which  uses 
Two  Phase  Heating  for  the  Unsaturated 
zone  and  the  Saturated  zone  and  then  uses 
Oxidation  combined  with  Biological 
Remediation  for  the  Aquifer 

Train  with  very  good  utility  scores.  Also, 
this  train  quickly  achieves  a  high  level  of 
performance.  This  train,  however,  was 
dominated  by  trains  no.  78  and  94  during 
the  phase  1  screening  analysis.  It  is  also  the 
most  expensive  of  the  three. 

135  -  This  is  a  U/S  &  A  Train  which  uses 
LAS  AGNA  for  the  Unsaturated  zone  and 
the  Saturated  zone  and  then  uses  Oxidation 
and  Biological  Remediation  for  the  Aquifer 

Train  maintains  good  utility  scores, 
achieves  a  high  level  of  performance  quickly 
and  is  less  expensive  than  train  no.  119. 

This  train,  like  no.  119,  was  shown  to  be 
dominated  by  trains  no.  142  and  148 

It  should  be  noted  that  the  most  salient  quality  of  the  three  trains  shown  in  Table 
5.1  was  performance  potential.  These  trains  were  able  to  perform  at  the  90%  or  better 
level  very  quickly,  and  thus  received  high  utility  scores  and  favorable  performance  ratings. 
The  remainder  of  the  trains  (4,  94,  148,  and  156)  were  found  to  be  much  less  competitive 
than  the  three  shown  in  Table  5.1. 

If  train  3  were  selected,  and  it  was  desired  to  reach  a  95%  cleanup  level,  the  costs 
are  expected  to  range  from  $0.9  million  to  $2.7  million  with  an  expected  value  of  $  1.6 
million.  The  time  required  for  this  remediation  scheme  is  expected  to  be  approximately 
4.5  years,  but  may  be  a  little  as  2. 1  years  or  as  much  as  8.3  years.  (See  Table  4.9  for  more 
in-depth  results  regarding  other  performance  levels). 

If  train  119  were  selected,  and  it  was  desired  to  reach  a  95%  cleanup  level,  the 
costs  are  expected  to  range  from  $3.0  million  to  $7.4  million  with  an  expected  value  of  $ 
4.8  million.  The  time  required  for  this  remediation  scheme  is  expected  to  be 
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approximately  3.6  years,  but  may  be  a  little  as  1.8  years  or  as  much  as  6.5  years.  (See 
Table  4.9  for  more  in-depth  results  regarding  other  performance  levels). 

If  train  135  were  selected,  and  it  was  desired  to  reach  a  95%  cleanup  level,  the 
costs  are  expected  to  range  from  $1.7  million  to  $5.6  million  with  an  expected  value  of  $ 

3.2  million.  The  time  required  for  this  remediation  scheme  is  expected  to  be 
approximately  4.2  years,  but  may  be  a  little  as  2.0  years  or  as  much  as  8.0  years.  (See 
Table  4.9  for  more  in-depth  results  regarding  other  performance  levels). 

Based  on  the  results  summarized  here,  it  can  also  be  concluded  that  the  objectives 
of  this  research  -  specifically  screening  candidate  trains  -  was  accomplished.  The  research 
conducted  provides  the  DOE  with  a  quantitative  justification  and  rationale  for  discarding 
55  out  of  58  (approximately  95%)  candidate  trains. 

5.3  Recommendations  for  Further  Research 

The  recommendations  for  further  research  are  shown  below: 

•  Data  -  More  accurate  data  is  needed.  As  mentioned  in  Chapter  3,  the  data  used  for  this 
analysis  was  derived  from  a  spill  volume  of  7,500  gallons.  At  present,  it  is  believed  that 
the  spill  volume  may  range  from  2,000  to  500,000  gallons  (see  Appendix  J).  The 
results  of  the  analysis  would  be  much  more  meaningfiil  to  DOE  DMs  if  the  data  more 
closely  approximated  reality.  Additionally,  the  cost,  time,  and  performance,  data  used 
in  this  analysis  does  not  reflect  the  difficulties  that  will  likely  be  encountered  due  to  the 
existing  facility  and  underground  wires  and  pipes.  Incorporating  these  factors  into  the 
data  would  make  the  analysis  more  valuable  to  the  DOE. 
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•  Uncertainty  -  This  is  a  side-bar  to  the  bullet  above.  If  data  is  to  be  collected, 

uncertainties  surrounding  the  data  should  also  be  collected  and  reported.  It  will  make 
future  analysis  more  meaningful  to  DOE  DM’s.  For  example,  capital  and  O&M  cost 
ranges  would  be  very  helpful  for  gaining  insight  into  the  annual  costs  expected  for  a 
given  train. 

.  CERCLA  Criteria  -  Rather  than  using  three  attributes  as  decision  criteria,  the  CERCLA 
criteria  that  has  been  established  for  Hazardous  Waste  Sites  (See  Appendix  B).  could 
be  used. 

.  Scaleable  data  -  Because  the  uncertainty  surrounding  the  spill  volume  for  this  site  is  so 
great,  data  showing  performance  levels  and  their  relationship  to  time  for  a  number  of 
potential  spill  volumes  would  be  very  valuable. 

.  Technical  risks  -  An  analysis  of  the  technical  risks  associated  with  each  of  the 

technologies  would  be  quite  helpful.  For  example,  what  are  the  risks  associated  with  a 
given  train  failing  to  meet  its  expected  performance  goal. 

From  this  research  and  analysis,  it  is  quite  clear  that  making  use  of  DA  and  MAUT 
techniques  can  be  a  valuable  approach  for  quantitatively  evaluating  competing  remediation 
trains.  The  DA  approach  assisted  the  DOE  in  their  understanding  of  the  decision  situation 
and  reducing  the  number  of  potential  trains  from  over  13  million  to  the  reasonable  number 
of  58.  The  MAUT  approach  provided  the  DOE  with  a  quantitative  means  to  further 
screen  the  remaining  58  to  a  list  of  three  top  contenders. 

In  addition  to  the  specific  details  provided  for  the  WAG6,  C-400  site,  this  research 
provides  the  DOE  with  a  framework  and  easy  to  follow  approach  to  evaluate  other. 
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similar  DNAPL  sites.  The  analysis  provided  in  this  effort  could  easily  be  adapted  by 
another  analyst  and  applied  to  other  sites  contaminated  with  DNAPLs  or  other  similar 
contaminants. 
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Appendix  A:  Specialized  DA/MAUT  Models 


According  to  Jennings,  Mehta,  and  Mohan,  a  number  of  specialized  MAUT/DA 
models  have  been  built  over  the  past  15  years.  Table  A.  1,  based  on  the  Jennings  et  al. 
work,  provides  a  basic  description  of  these  models  and  identifies  their  specialized  use. 


Table  A.  1  MAUT  Models  Used  for  Specialized  DA  [Jennings  et  al.,  1994:  1133  -  1 135] 


Model 

BS!S9 

Description 

References 

HRS 

(Hazardous 

Ranking 

System) 

EPA 

The  model  is  (was)  used  to  place 
waste  sites  on  the  National  Priorities 
List  (NPL).  It  is  a  simple  model  based 
on  a  composite  weighted  score: 

S  =  Sm  +  Sfe  +  Sdc 
where  Sm  is  the  score  for  potential 
migration,  Sfe  is  the  score  for 
potential  explosion  or  fire  hazards,  and 
Sdc  is  the  score  for  the  potential  to 
encounter  a  direct  contact  with  the 
hazard.  The  model  was  used  to 
generate  the  original  list  of  418  NPL 
sites  in  1982. 

Wu  and  Hilger, 
1984:  797  -  807 

Wells  and 

Caldwell,  1990: 
71-76 

MEPAS 

(Multimedia 

Environmental 

Pollutant 

Assessment 

System) 

DOE 

Used  to  assess  wildlife  endangerment 
problems  at  DOE  sites 

Hartz  and  Whelan, 
1988:  295  -  299 

Droppo  and 

Hopes,  1990:  193 
-205 

RAAS 

(Remedial 

Action 

Assessment 

System) 

DOE 

Used  during  the  Feasibility  Study 
phase  to  provide  alternatives  guidance 

Hartz  and  Whelan, 
1988:  295  -299 

HAZRISK 

DOE 

Used  in  developing  cost  estimates  and 
scheduling  for  projects  involving 
hazardous  waste  cleanup 

Hudson  and 
Shangraw,  1990: 

241  -  244 
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Model 

Description 

References 

POS  (Program 

Optimization 

System) 

DOE 

Used  to  optimize  spending  on 
restoration  projects  and  to  allocate 
remediation  budgetary  resources.  The 
model  is  based  on  the  following  utility 
function: 

U  =  WhsUhsXhs  +  w„u„x„  + 

WpcUpcXpc  -  W„=Crc  -  WfcCfc 

Where  the  X’s  are  performance  scores 
for  health  and  safety  (hs),  regulatory 
responsiveness  (rr),  and  public 
concern  (pc).  The  U’s  represent 
functions  that  define  the  utility  of  each 
X  score.  The  W’s  are  the  importance 
weights,  and  the  C’s  are  the  remaining 
(rc)  and  future  (fc)  costs. 

Merkhofer, 

Cotton,  Longo, 
1988;  39-43 

DPM  (Defense 
Priority  Model) 

DOD 

Used  to  estimate  the  risk  to  human 
health  and  the  environment.  Makes 
use  of  weighted  scores  that  consider 
the  source  materials  and  their 
pathways  to  humans  and  the 
environment. 

Expert,  1990 

Air 

Force 

Predecessor  to  DPM 

Hushon,  1990: 
206-216 

FLEX  (Flexible 
Liner  Expert) 

Used  to  evaluate  chemical 
compatibility  of  liners 

Rossman  and 

Siller,  1987:  113  - 
127 

HERPM 

(Human 

Exposure 

Potential 

Ranking  Model) 

NY  Dept 
of  Health 

Used  as  a  ranking  tool  to  establish 
relative  priorities  for  investigation  and 
remediation  of  sites 

Smith,  Patrick, 
and  Hudson, 

1987:  158-  161 

DRASTIC 

EPA 

Used  to  preliminary  assess  hazardous 
waste  sites.  Acronym  is  from; 

Depth  to  water,  net  Recharge,  Aquifer 
media.  Soil  media,  Topography, 

Impact  of  the  vadose  zone  media,  and 
hydraulic  Conductivity  model 

Allert,  Bennet, 

Lehr,  Petty,  and 
Hackett,  1987 

Model 

Description 

References 

CORA  (Cost  of 

Remedial 

Action) 

Used  as  a  costing  tool  for  remedial 
actions 

Chenu  and 

Crenca,  1990:  162 
-  175 

PAST 
(Potential 
ARAR’s 
Selection  Tool) 

1 

Used  to  assist  in  the  development  of 
ARAR’s  (Applicable  or  Relevant  and 
Appropriate  Requirements  for 
remediation  projects 

Greathouse  and 
Clements,  1991: 
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Appendix  B:  National  Contingency  Plan  Evaluation  Criteria 


There  are  nine  criteria  used  to  evaluate  remediation  projects  [DOE/OR/07-1243&D3]. 

These  nine  are  listed  and  explained  below: 

Threshold  Criteria:  (each  alternative  must  meet  to  be  the  selected  alternative) 

1.  Overall  Protection  of  human  health  and  the  environment 

.  Addresses  if  an  alternative  can  adequately  protect  human  health  and  the  environment  in 
both  the  short  and  long-term,  from  unacceptable  risks  posed  by  hazardous  substances, 
pollutants,  or  contaminants  present  at  the  site  by  eliminating,  reducing,  or  controlling 
exposures  to  contamination. 

2.  Compliance  with  Applicable  or  Relevant  and  Appropriate  Requirements  (ARAR’s) 

.  Assesses  if  the  alternatives  attain  ARARs  under  Federal  environmental  laws  and  State 
environmental  laws. 

Balancing  Criteria:  (used  to  compare  and  weigh  the  major  trade-offs  among  the 

alternatives) 

3.  Long-term  effectiveness  and  permanence 

.  Assesses  the  ability  of  an  alternative  to  provide  long-term  protection  after  remediation 
goals  have  been  met. 

4.  Implementability 

.  Evaluates  the  ease  or  difficulty  of  implementing  the  alternatives. 

5.  Reduction  of  toxicity,  mobility,  or  volume  through  treatment 
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.  Evaluates  the  degree  to  which  alternatives  employ  recycling  or  treatment  that  reduces 
toxicity,  mobility,  or  volume,  including  how  treatment  is  used  to  address  the  principal 
threats  posed  by  the  site. 

6.  Short-term  effectiveness  and  permanence 

.  Addresses  the  time  it  takes  for  an  alternative  to  be  implemented  and  the  potential  effect 
on  human  health  (including  the  community  and  workers)  and  the  environment  during 
implementation. 

7.  Cost 

.  Includes  capital  costs,  including  both  direct  and  indirect  costs,  annual  operation  and 
maintenance  (O&M) 

Modifying  Criteria:  (used  to  modify  an  alternative  or  select  a  different  alternative  based 

on  input  from  the  State  or  the  public) 

8.  State  Acceptance 

.  Considers  the  concerns  of  the  State  on  the  alternatives. 

9.  Community  Acceptance: 

•  Allows  for  a  public  comment  period  when  interested  persons  and  organizations  can 
comment  on  the  proposed  remedy. 


Adapted  from  DOE/OR/07-1243&D3  and  GJPO-TP-11 
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Appendix  C:  Acronym  Listing 


CERCLA 

Comprehensive  Environmental  Response,  Compensation,  and  Liability  Act 

CL 

Contaminant  Level 

DA 

Decision  Analysis 

DM 

Decision  Maker 

DNAPL 

Dense  Non- Aqueous  Phase  Liquid 

DOE 

Department  of  Energy 

DPL 

Decision  Programming  Language 

DQO 

Data  Quality  Objective 

ERBAM 

Environmental  Restoration  Benefit  Assessment  Matrix 

EV 

Expected  Value 

MAUT 

Multiattribute  Utility  Theory 

NCP 

National  Contingency  Plan 

NPL 

National  Priorities  List 

NPV 

Net  Present  Value 

O&M 

Operations  and  Maintenance 

PCB 

Polychlorinated  Biphenyl 

RF 

Radio  Frequency 

RI/FS 

Remedial  Investigation/Feasibility  Study 

SO 

Site  Owner 

SVE 

Soil  Vapor  Extraction 

TCE 

Trichloroethylene 

TIDM 

Technology  Investment  Decision  Model 

TOC 

Total  Organic  Content 

VOC 

Volatile  Organic  Compound 

WAG6 

Waste  Area  Group  6 
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Appendix  D:  Technology  Candidates 


This  appendix  contains  a  description  of  the  candidate  remediation  technologies  that 
have  been  selected  for  consideration  at  the  C-400  site. 

6  Phase  Heating 

This  technology  uses  a  six-pointed  star  configuration  of  electrodes  placed  in  the 
ground  to  enhance  the  mobility  of  contaminants.  When  power  is  applied  to  the  electrodes, 
the  resistive  heat  volatizes  the  contaminants  and  produces  steam  that  in-tum  strips  the 
pollutant  in-situ.  The  volatized  contaminants  are  then  removed  via  SVE. 

The  unique  six  pointed  star  configuration  allows  for  the  uniform  heating  of  the 
soil.  The  patented  process  breaks  the  conventional  3  phase  electricity  into  a  six  phase 
system.  This  provides  for  less  expensive  soil  heating  that  is  estimated  to  be  approximately 
100  kWh  per  cubic  meter  of  soil  [DOE/EM-0248,  1995:  206]. 

Biological  Treatment 

Biological  treatment  or  bioremediation  makes  use  of  existing  microorganisms, 
primarily  bacteria,  to  effectively  destroy  contaminants.  The  key  to  this  technology  is  the 
injection  of  nutrients  and  air  into  the  soil  to  stimulate  the  microorganisms  growth.  When 
the  microorganisms  are  given  the  proper  stimuli  (nutrients),  they  can  transform  the 
contaminants  to  harmless  products  such  as  carbon  dioxide  and  water.  Because  the 
microorganisms  can  be  quite  sensitive  to  changes  in  hydrogeology,  a  good  description  of 
the  subsurface  is  necessary  prior  to  injecting  nutrients.  Additionally,  high  concentrations 
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of  contaminants  may  not  be  affected  by  microbial  growth,  and  may  even  kill  fledgling 
microbes.  It  is  for  this  reason  that  field  tests  have  shown  that  bioremediation  works  best 
around  the  periphery  of  a  highly  concentrated  contaminant  zone  [National  Research 
Council,  1995:  131  -  133], 

Cosolvent  Enhanced  In-situ  Treatment 

Cosolvents,  like  surfactants,  are  used  to  enhance  conventional  pump  and  treat 
systems.  Cosolvents  are  substances  that,  when  mixed  with  water,  can  enhance  the 
solubility  of  the  contaminant.  Typical  cosolvents  are;  methanol,  and  acetone.  The  idea 
behind  this  technology  is  similar  to  that  of  surfactant  flooding.  The  cosolvent  agent  is 
injected  into  the  ground  and  the  area  is  flooded.  Then  the  pump  and  treat  systems  pump 
the  liberated  contaminant  and  the  cosolvent  from  the  ground  for  treatment.  The  cosolvent 
acts  to  increase  the  solubility  of  the  contaminant  as  well  as  decrease  the  retardation  factors 
present  in  the  contaminant,  especially  in  DNAPLs.  An  additional  benefit  of  cosolvents  is 
that  occasionally  (depending  on  the  concentration  levels  of  the  cosolvent  and  the 
contaminant  as  well  as  the  hydrogeology)  bioremediation  is  stimulated  and  microbial 
growth  occurs.  This  Avill  increase  the  in-situ  destruction  of  the  contaminant  [National 
Research  Council,  1995:  148  -  149]. 

Dynamic  Underground  Stripping 

This  technology  is  actually  a  combination  of  three  separate  technologies: 
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1.  Steam  injection  -  Steam  is  injected  around  the  periphery  of  a  contaminated  site  to 
heat  the  subsurface  and  potentially  vaporize  the  contaminants.  This  is  done  in  the  hopes 
that  the  contaminants  will  be  driven  to  a  central  location  and  extracted  via  a  SVE  scheme. 

2.  Electrical  Heating  -  Electrical  probes  are  used  in  the  areas  that  contain  the  highest 
concentration  of  contaminants.  The  heating  will  liberate  the  contaminants  from  the  soil 
and  drive  them  into  the  same  area  that  the  steam  injection  technology  drove  the  vaporized 
contaminants. 

3.  Underground  Imaging  -  Electrical  Resistance  Tomography  is  used  as  an  imaging 
technique  to  provide  assurance  that  the  heated  areas  are  cleaning  up  the  site  as  expected. 
This  technology  can  be  used  above  and  below  the  water  table,  and  is  especially  well  suited 
for  subsurfaces  where  clay  and  sand  are  interlaced  together  [Mather,  1995:  WWWeb]. 

LASAGNA™  or  Electro-Osmosis 

Electro-osmosis  is  a  process  that  uses  electrodes  placed  into  the  soil  to  mobilize 
DNAPL  contaminants.  Once  power  is  applied  to  the  electrodes,  the  DNAPL  contaminant 
will  migrate  in  sufficient  quantities  to  warrant  some  form  of  removal  or  in-situ  destruction. 
Typically,  the  contaminant  is  removed  via  adsorption  or  destroyed  biologically.  The 
technology  appears  to  be  most  beneficial  when  contaminants  are  located  in  the  saturated 
zone  and  where  the  soil  has  a  low  permeability. 

The  term  LAS  AGNA™  was  derived  through  a  consortium  of  private  companies 
that  were  researching  methods  to  speed  VOC  contaminant  cleanup.  The  researchers 
determined  that  a  layered  application  of  Electro-osmotic  probes  would  speed  the 
migration  of  the  contaminants  and  thus  increase  the  speed  at  which  cleanup  could  occur. 
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It  is  because  of  this  layered  application  that  the  term  LASAGNA™  was  bom. 
Subsequently,  Monsanto  (one  of  the  consortium  members)  patented  and  trademarked  the 
process  [Falta  et  al.  1996:  24], 


Oxidation 

Chemical  Oxidation  is  an  in-situ  remediation  technology  that  uses  a  chemical 
oxidant  solution  such  as  hydrogen  Peroxide  (H2O2)  or  Potassium  Permanganate  (KMn04) 
as  an  oxidant  to  degrade  organic  contaminants  (like  DNAPLs).  The  oxidant  solution  (i.e. 
H2O2  or  KMn04)  is  either  injected  into  the  ground  or  is  mixed  with  the  soil  through  a  soil 
mixing  apparatus.  Either  way,  the  oxidant  comes  into  contact  with  the  organic  DNAPL 
contaminant,  and  destroys  it  [West,  1996]. 

Pump  and  Teat 

Conventional  pump  and  treat  systems  operate  by  pumping  ground  water  to  the 
surface,  removing  the  contaminants  (typically  through  some  type  of  air  stripper)  and 
returning  the  water  to  the  ground  or  into  a  municipal  sewage  plant.  This  technology  has 
been  labeled  old  fashioned  and  inefficient  because  of  the  large  volumes  of  water  that  must 
be  pumped  out  of  the  ground.  Also,  due  to  the  chemical  composition  of  DNAPL 
contaminants,  the  fact  that  one  has  pumped  contaminated  water  out  of  the  ground  does 
not  guarantee  the  complete  cleanup  of  the  site.  The  residual  contaminants  that  adhere  to 
subsurface  particles  are  not  as  susceptible  to  pump  and  treat  systems  [National  Research 
Council,  1995:  29]. 
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Permeable  Treatment  Zones 


This  technology  makes  use  of  a  permeable  “wall”  that  is  inserted  into  the  subsurface.  The 
wall  allows  ground  water  to  flow  through  it.  As  the  groundwater  flows  through,  the 
composition  of  the  encourages  a  chemical  treatment  of  the  contaminant.  The  most 
common  type  of  PTZ  wall  is  iron  filings.  The  iron  filings  create  a  chemical  reaction  with 
the  contaminant  (i.e.  TCE)  and  encourage  is  breakdown.  As  the  water  flows  through  the 
permeable  barrier,  the  contaminant  is  trapped  in  the  chemical  reactions  of  the  wall 
[Clayton,  1997]. 

Radio  Frequency  (RF)  Heating 

Radio  Frequency  (RF)  heating  uses  the  heat  energy  induced  by  the  application  of 
RF  energy  into  the  soil  to  enhance  conventional  vapor  extraction  methods.  The  heat 
applied  to  the  soil  through  RF  causes  a  liberation  of  the  contaminants,  especially  Volatile 
Organic  Compounds  (VOCs)  like  TCE.  This  technology  is  most  applicable  to  remediation 
of  the  vadose  zone  [DOE/EM-0248,  1995:  215  -  217]. 

Surfactant  Enhanced  In-situ  Treatment 

Surfactants  (surface  active  agents)  are  used  to  enhance  conventional  pump  and 
treat  systems.  They  are  categorized  as  part  of  the  family  of  soil  flushing.  The  idea  behind 
this  technology  is  to  inject  a  surfactant  into  the  ground  and  flood  the  area  with  a  surfactant 
agent.  Once  the  flood  is  complete,  conventional  pump  and  treat  systems  are  used  to 
recover  both  the  contaminant  and  the  surfactant.  The  surfactant  acts  as  a  loosening  agent 
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to  separate  the  contaminant  from  the  water  saturated  soil  [Falta  et  al.,  1996:26  -  29  and 
National  Research  Council,  1995:  148-  149]. 

Soil  Vapor  Extraction  (SVE) 

This  technology  uses  extraction  wells  to  draw  gases  from  the  ground  in  a  vacuum 
style  configuration.  The  gases  drawn  are  typically  air  bubbles  that  have  picked  up 
contaminants  in  their  journey  to  the  extraction  well  head.  The  air  bubbles,  along  with  their 
molecules  of  contaminant  are  sucked  into  the  vacuum/blower  system,  and  treated  above 
ground  in  an  off-gas  treatment  system.  This  method  is  applicable  above  the  water  table 
(unsaturated  zone)  in  highly  permeable  soils.  It  is  typically  used  in  conjunction  with  an 
enhancing  technology  such  as  air  sparging  or  thermal  heating  technologies.  SVE  has  been 
used  for  over  ten  years,  especially  in  the  remediation  of  petroleum  hydrocarbons  [Falta  et 
al.,  1996:  40-41] 

Two  Phase 

This  technology  makes  use  of  a  powerful  vacuum  system  that  extracts  soil  gases 
and  liquids.  It  is  typically  used  in  low  permeability  subsurfaces.  As  the  vacuum  is  applied 
through  the  screened  well  head,  soil  vapors  are  extracted  and  groundwater  is  entrained  by 
the  extracted  vapors.  Once  the  vapors  and  groundwater  are  above  ground,  they  are 
separated  and  treated  individually.  This  technology  can  be  combined  with  bioremediation, 
air  sparging,  or  bio-venting.  It  is  generally  considered  to  be  more  effective  than  using 
SVE  singularly  [EPA/542/B-94/013,  1994:  4.145  -  4.147]. 
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UVB 


This  technology  is  similar  to  pump  and  treat  except  that  it  does  not  use  a  well  head  like 
pump  and  treat  does.  The  technology  is  used  to  pump  contaminated  groundwater  to  the 
surface,  treat  it  using  some  form  of  air  stripping,  and  then  injecting  it  back  into  the 
subsurface.  Typically  this  technology  is  used  to  affect  a  small  area  no  larger  than  50’  in 
diameter.  This  is  compared  against  pump  and  treat  which  can  affect  an  area  up  to 
approximately  300’  [Clayton,  1997]. 
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Appendix  E:  Phase  1  (Screening  Model)  Decision  Tree 


The  complete  screening  model  decision  tree  is  shown  in  Figure  E.  1.  Figures  E.2 
through  E.6  are  enlarged  views  of  each  subsection  of  the  overall  decision  tree. 


Figure  E.  1  Decision  Tree 


Figure  E.2  Enlarged  Section  #1  of  Decision  Tree 
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Figure  E.4  Enlarged  Section  #3  of  Decision  Tree 
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Figure  E.5  Enlarged  Section  #4  of  Decision  Tree 
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Figure  E.6  Enlarged  Section  #5  of  Decision  Tree 
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Train  descriptions  are  shown  in  Table  F.l 
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Two-Phase  Vacuum  Extraction _ None _ Two-Phase  Vacuum  Extraction  None _ Surfactant  Flushing _ None 

Two-Phase  VacuunrT Extraction  None _ Two-Phase  Vacuum  Extraction  None _ UVB _ Biotogical 

Two-Phase  Vacuum  Extraction  None _ Two-Phase  Vacuum  Extraction  None _ UVB _ None 

Two-Phase  Vacuum  Extraction _ None _ Two-Phase  Vacuum  Extraction  None _ Pump  and  Treat _ Biological 

Two-Phase  Vacuum  Extraction  None _ Two-Phase  Vacuum  Extraction  None  Pump  and  Treat  None 


Table  F.2  Attribute  Values 


Train  NPV  O&M 

No.  Cost_ Cost  Time  Perf 


3  * 

ESI 

mm\ 

4  * 

2.542,477  | 

203,820 

0.50  1 

90%  1 

71 

737,648 

■iiif 

72 

604,910 

lEai 

eh 

73 

gWsjciMi 

74 

M.I.LliliiiliM 

75 

7,771,542 

76 

7,798.676 

700,845 

6.00 

86% 

77 

5,942.740 

99% 

78 

■sESmil 

■cKiIirf 

99% 

79 

5.568,950 

426,315 

80 

609.335 

77%  1 

83 

426.315 

muktM 

eh 

84 

426,315 

wai 

eh 

87 

5,759.190 

737,648 

99%  1 

88 

604,910 

■il-M 

89 

8.840.448 

na 

EH 

90 

8,867,581 

694,365 

KH 

91 

7,749,117 

419,835 

EH 

mm 

7,776,251 

694,365 

EH 

93 

5,920,314 

430,126 

99% 

94  * 

4,347,951 

99% 

95 

5,546,525 

Ba 

96 

602,855 

g*!«M 

Wkkii^M 

99 

■SniEESi 

W^K¥M 

MiWM 

100 

419,835 

119  * 

4,635,890 

1,056,083 

EH 

120 

923,345 

mmM 

EH 

121 

8,647,315 

725,703 

85%  1 

122 

mkLimm 

592,965 

MjMM 

123 

7.555,984 

WESESik 

124 

6,652.951 

125 

748,561 

■Clljll 

126 

3,224,651 

■ci»M 

m:mm 

127 

5.353,392 

725,703 

M»!>1 

128 

■6W;»»klrM 

501,455 

131 

M:WM 

132 

3,892,017 

318,435 

MiWM 

■ESI 

2.949,999 

737,648 

4.56 

99% 

136 

604,910 

99% 

09 

6.058,390 

■eVZI.^ 

DEIll 

4,939,926 

315,758 

lEISi 

183,021 

141 

■tlHKMtM 

MiWM 

142 

1,538,760 

lEE9 

143 

2,737,334 

MEM 

95%  1 

144 

ubesesi 

BEI 

miwm 

147 

I  3,071,842 

■ipciwcr^ 

mem 

Kill 

■ESI 

lEH 

151 

I  3.171,860 

M»!»l 

64% 

152 

IBJUrWJ 

11^ 

80% 

153 

IBB 

IHS 

5,353,392 

lESI 

■^1 

2,888,180 

316,627 

1X3 

IKEH 

157 

4.926,754 

517,825 

11^ 

1  50%  1 

158 

4,904.329 

mm 

160 

3.781,029 

1  501,455 

1  6.00 

1  42%  1 

All  of  the  values  shown  in  table  F.2 
were  acquired  through  Mountain  State 
Engineering  (MSE)  Corporation. 

They  were  generated  using  a  life  cycle 
cost  model  owned  and  operated  by  MSE. 
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Appendix  G  Model  Iterations 


Iteration  Number  1 

At  the  onset  of  the  project,  the  emphasis  was  to  build  a  model  that  contained 
hydrogeological  characteristics  and  their  relationships  to  cost,  time,  and  performance 
parameters  for  each  strategy.  The  strategies  are  divided  into  two  zones  -  the  saturated, 
and  the  unsaturated  zone.  Figure  G.  1  shows  an  influence  diagram  of  this  early  iteration. 

From  Figure  G.  1,  it  is  seen  that  there  are  two  decisions  -  which  technology 
strategy  to  use  for  the  unsaturated  zone,  and  which  one  to  use  for  the  saturated  zone. 
Based  on  these  decisions,  the  cost,  time,  and  technology  effectiveness  (performance)  are 
determined  and  combined  into  a  utility  value.  The  utility  scores  for  the  unsaturated  zone 
and  the  saturated  zone  are  then  combined  into  an  overall  utility  score  and  reported  as  one 
number.  Note  that  the  hydrogeological  considerations  for  both  the  unsaturated  and 
saturated  zone  influence  the  technology  effectiveness.  This  was  an  important  modeling 
consideration  in  this  early  version. 

This  model,  however,  was  discarded  because  of  a  lack  of  relevant  data  relating 
hydrogeologic  characteristics  to  technology  effectiveness  parameters. 
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Figure  G.l  Influence  Diagram,  Iteration  1 


Iteration  Number  2 


The  second  iteration  contained  a  number  of  changes  -  most  notably  the  deletion  of 
influencing  relationships  between  hydrogeological  characteristics  and  effectiveness  levels 
and  the  addition  of  three  more  zones.  Figure  G.2  shows  an  influence  diagram  of  this 
iteration. 


Figure  G.2  Influence  Diagram,  Iteration  2 
From  Figure  G.2  it  is  seen  that  the  decision  situation  is  broken  into  5  separate 
zones.  Each  zone  contains  one  decision  -  which  technology  to  use  in  that  particular  zone. 
The  zones  are  segregated  based  on  hydrogeological  characteristics,  and  cost  (operating 
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and  capital),  time,  and  performance  quantities  are  identified  for  each  technology  for  each 
zone.  These  cost,  time,  and  performance  figures  are  then  combined  to  determine  totals 
(cost,  time,  and  performance)  for  a  complete  strategy. 

The  model  was  discarded  because  of  two  major  issues.  The  first  was  the  number 
of  potential  combinations.  Nine  different  technologies  were  identified  for  the  first  four 
zones,  and  eighteen  were  identified  for  zone  5.  This  meant  that  there  were  1 18,098  (9^*  x 
18)  possible  trains.  This  was  far  too  many  to  evaluate  individually,  so  an  algorithm  was 
developed  to  individually  combine  the  4  parameters  (operating  cost,  capital  cost, 
performance,  and  time).  This  led  to  the  second  issue.  Because  of  the  diversity  of  the 
technologies,  a  reasonable  algorithm  that  combines  individual  technology  parameters 
could  not  be  agreed  upon.  The  decision  makers  at  the  DOE  were  unable  to  accept  the 
simple  algorithms  that  were  developed  and  could  not  provide  alternative  ones. 

Iteration  Number  3 

Iteration  no.  3  was  a  more  accurate  reflection  of  the  decision  situation.  It 
contained  a  more  realistic  set  of  potential  strategies.  The  influence  diagram  of  this  model 
is  shown  in  Figure  G.3 . 
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From  this  Figure,  it  is  seen  that  the  decision  situation  was  been  broken  into  three 
distinct  pieces.  Each  piece  represented  a  strategy  style.  This  was  a  major  change  in  the 
approach  to  this  problem.  The  previous  two  versions  focused  on  remediating  specific 
zones  or  layers  of  subsurface.  This  model  focused  on  realistic  strategies  developed 
through  the  iterative  modeling  process  and  have  been  categorized  into  one  of  three  styles. 
These  styles  were  thought  to  better  depict  the  decision  making  thought  process  used  at  the 
DOE.  The  three  styles  were: 
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.  USA  (Silver  Bullet  -  found  in  the  upper  right  hand  comer  of  Figure  G.3)  -  This 
strategy  style  will  remediate  all  of  the  contaminated  zones.  Note  that  at  the 
time  of  this  iteration,  the  USA  trains  were  designated  as  SB  or  Silver  Bullet 
Trains. 

.  US  &  A  (Unsaturated,  Saturated,  and  Aquifer  -  found  in  the  center  of  Figure 
G.3)  -  This  strategy  style  uses  one  technology  suite  for  the  unsaturated  and 
saturated  zones,  and  another  suite  for  the  aquifer. 

.  S  A  &  U  (Saturated,  Aquifer,  and  Unsaturated  -  found  on  the  right  hand  side  of 
Figure  G.3)  -  This  strategy  style  uses  one  technology  suite  for  the  saturated 
zone  and  the  aquifer,  and  uses  another  suite  for  the  unsaturated  zone. 

An  explanation  of  each  of  the  nodes  shown  in  Figure  G.3  is  found  below  in  Table  G.  1 . 
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Table  G.  1  Node  Explanation  for  Figure  G.3 


1  Node  Identifier 

Description 

Decision 

The  decision  here  is  between  SB,  US  &  A,  and  SA  & 

U 

Which  silver 
bullet 

Decision 

The  decision  is  between  3  SB  primary  technologies 

Which  silver 
bullet 

accompaniment 

Decision 

The  decision  is  between  2  SB  technology 
accompaniments 

SB  cost 

Value 

Contains  cost  data  for  SB  strategies 

SB  time 

Value 

Contains  time  data  for  SB  strategies 

SB  perf 

Value 

Contains  performance  data  for  SB  strategies 

Which  U  and  S 

Decision 

The  decision  is  to  select  one  of  3  primary 
technologies  to  be  used  in  the  unsaturated  and 
saturated  zones 

which  U  and  S 
accompaniment 

Decision 

The  decision  is  to  select  one  of  4  technology 
accompaniments  to  be  used  in  the  unsaturated  and 
saturated  zones 

One  of  sixteen 

aquifer 

technologies 

Decision 

The  decision  is  to  select  one  of  16  different 
technologies  for  the  aquifer 

US  &  A  cost 

Value 

Contains  cost  data  for  US  &  A  strategies 

US  &  A  time 

Value 

Contains  time  data  for  US  &  A  strategies 

Value 

Contains  performance  data  for  US  &  A  strategies 

Which  S  and  A 
tech 

Decision 

The  decision  is  to  select  one  of  two  primary 
technologies  for  the  saturated  zone  and  the  aquifer 

which  S  and  A 
accompaniment 

Decision 

The  decision  is  to  select  one  of  two  technology 
accompaniments  for  the  saturated  zone  and  the 
aquifer 

which  unsat  tech 

Decision 

The  decision  is  to  select  one  of  two  primary 
technologies  for  the  unsaturated  zone 

which  unsat  tech 
accompaniment 

Decision 

The  decision  is  to  select  one  of  four  technology 
accompaniments  for  the  unsaturated  zone 

S  A  &  U  cost 

Value 

Contains  cost  data  for  S  A  &  U  strategies 

SA  &  U  time 

Value 

Contains  time  data  for  SA  &  U  strategies 

Value 

Contains  performance  data  for  SA  &  U  strategies 

total  time 

Value 

Represents  the  total  amount  of  time  required  for  a 
strategy 

Value 

Represents  the  total  performance  provided  by  a 
strategy 

Represents  the  total  cost  required  for  a  strategy 
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Each  technology  suite  (SA,  US,  A,  U,  SB)  is  composed  of  a  primary  technology 
and  a  technology  accompaniment.  For  example,  one  of  the  SA  suites  is  PTZ  and  Bio. 
PTZ  (Permeable  Treatment  Zones)  is  the  primary  technology,  and  bioremediation  is  the 
technology  accompaniment. 

Because  of  the  complete  enumeration  of  the  data  for  each  of  the  variables  (cost, 
time,  and  performance),  the  details  of  the  decision  situation  were  model  in  a  decision  tree. 
The  decision  tree  was  developed  in  conjunction  with  the  influence  diagram.  The  structure 
of  the  decision  tree  is  shown  in  Figure  G.4.  The  purpose  of  Figure  G.4  is  to  show  the 
overall  structure  of  the  decision  tree,  not  the  details  of  the  tree.  A  section  of  the  detailed 
tree  is  found  in  Figure  G.5. 
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Figure  G.4  Decision  Tree  Structure 


From  this  point,  the  requisite  screening  model  was  built.  This  is  completely  described  in 
Chapter  3. 
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Appendix  H  DPL  Listing  Models 


Screening  Model: 


String  Excel_l="c:\\\vinword\\tony\\thesis\\testl07.xls"; 
excel(Excel_l,"Sheetl!costr’)  value  costl=1421 159.08665278; 
excel(ExceI_l /'Sheet  1  Icostl  1”)  value  costl  1=5942739.58435561; 
excel(Excel_l/'Sheetl!costl2")  value  costl 2=4370376.34732296; 
excel(Excel_l/'Sheetl!costl3")  value  costl 3=5568950.23788 13; 
excel(Excel_l/'Sheetl  Icostl 5")  value  costl 5=5903458. 12930399; 
excel(Excel_l/’Sheetl!costl6")  value  costl 6=5037742.04084865; 
excel(Excel_l/'Sheetl  IcostlT')  value  costl7=5759 190. 1278813; 
excel(Excel_l/'Sheetl  IcostlS")  value  costl8=4856 156.89084865; 
excel(Excel_l/'Sheetl  !costl9”)  value  cost 19=8840447. 5 8689769; 
excel(Excel_l/'Sheetl!cost2")  value  cost2=2542476.76712963; 
exceI(Excel_l/'Sheetl  !cost21")  value  cost2 1=77491 17.19963231; 
exceI(Excel_l/'Sheetl!cost23")  value  cost23=59203 14.47435561; 
excel(Excel_l /'Sheet  1  !cost24")  value  cost24=434795 1.23732296; 
excel(Excel_l/'Sheetl!cost25")  value  cost25=5546525. 1278813; 
excel(Excel_l/'Sheetl!cost27")  value  cost27=5881033. 01930399; 
excelCExcefl/'Sheetl  !cost28")  value  cost28=50 153 16.93084865; 
excel(Excel_l /'Sheet  1  !cost29")  value  cost29=4685890. 1278813; 
excel(Excel_l/'Sheetl!cost3")  value  cost3=3 172873.91703265; 
excel(Excel_l/'Sheetl!cost30")  value  cost3 0=3782 8 56. 890 84865; 
excel(Excel_l/'Sheetl  !cost32'')  value  cost32=68641 14.34986504; 
excel(Excel_l/'Sheetl!cost33")  value  cost33=6675817.19963231; 
excel(Excel_l/'Sheetl!cost35")  value  cost35=4847014.47435561; 
excel(Excel_l/'Sheetl!cost36")  value  cost36=327465 1.23732296; 
excel(Excel_l/' Sheet  1  !cost39")  value  cost39=4807733. 01930399; 
excel(Excel_l/'Sheetl!cost4")  value  cost4=2530677.8175; 
excel(Excel_l/'Sheetl!cost40")  value  cost40=39420 16.93084865; 
excel(Excel_l/’Sheetl!cost41")  value  cost41  =4635890. 1278813; 
excel(Excel_  1/' Sheet  l!cost42")  value  cost42=3732856. 89084865; 
excel(Excel_l/' Sheet  1  !cost47")  value  cost47=4797014.47435561; 
excel(Excel_l/'Sheetl!cost48")  value  cost48=322465 1.23732296; 
excel(Excel_l /'Sheet!  !cost5")  value  cost5=5781615.2378813; 
excel(Excel_l /'Sheet  1  !cost5 1")  value  cost5 1=4757733.01930399; 
excel(Excel_l/'Sheetl!cost52")  value  cost52=3892016.93084865; 
excel(Excel_l/'Sheetl!cost53")  value  cost53=2949998.7686409; 
excel(Excel_l/' Sheet  1  !cost54")  value  cost54=2046965.53 160825; 
excel(Excel_l/'Sheetl  !cost55")  value  cost5  5=603 1256.22765729; 
excel(Excel_l/'Sheetl  !cost57")  value  cost57=4939925.84039191; 
excel(Excel_l/'Sheetl  !cost58")  value  cost5  8=4036892.60335926; 
excel(Excel_l/'Sheetl  !cost59")  value  cost59=311 1 123.1 15 1 1522; 
excel(Excel_l/'Sheetl!cost6")  value  cost6=4878582.00084865; 
excel(Excel_l /'Sheet  1  !cost60")  value  cost60=l  538759.87808256; 
excel(Excel_l/’Sheetl!cost6r')  value  cost6 1=2737333. 7686409; 
excel(Excel_l/'Sheetl  !cost63")  value  cost63=307 1841. 660063  59; 
excel(Excel_l/'Sheetl  !cost64’')  value  cost64=2206 125.57160825; 
excel(Excel_l/' Sheet  1  !cost7")  value  cost7=8862872. 69689769; 
excel(Excel_l/’Sheetl!cost9")  value  cost9=777 1542.30963231; 
excel(Excel_l/'Sheetl!export_cost_wt")  value  export_cost_wt=0.44; 
excel(Excel_l/'Sheetl  !export_perf_wt")  value  export_perf_wt=0.44; 
excel(Excel_l /'Sheet  1  lexporttimewt")  value  export_time_wt=0.12; 
excel(Excel_l/'Sheetl!oml")  value  oml=336557.706942313; 
excel(Excel_l,"Sheetl  loml  1")  value  oml  1=430126.057378652; 
excel(Excel_l,"Sheetl!oml2")  value  oml2=632 192.780065968; 
excel(Excel_  1/' Sheet  1  lorn  13")  value  oml 3=4263 14.8; 
excel(Excel_l/'Sheetl!oml5")  value  oml 5=4263 14.8; 
excel(Excel_l/' Sheet  1  !oml6")  value  oml6=4263 14.8; 
excel(Excel_l/'Sheetl!oml7")  value  oml7=737648.0773 12683; 
excel(Excel_l/'Sheetl!oml8")  value  oml8=604910.333333333; 
excel(Excel_l/'Sheetl!oml9")  value  oml9=419834.8; 
excel(Excel_l,"Sheetl  !om2")  value  om2=203819. 962962963; 
excel(Excel_l/'Sheetl  !oni21")  value  om2 1=4 19834.8; 
excel(Excel_l/'Sheetl!om23")  value  om23=430 126.057378652; 
excel(Excel_l,"Sheetl  !om24")  value  om24=6257 12.780065968; 
excel(Excel_l/'Sheetl  !om25")  value  om25=4 19834.8; 
excel(Excel_l/'Sheetl!om27")  value  om27=4 19834.8; 
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excel(Excel_l/'Sheetl!oni28'*)  value  om28'-419834.8; 
excel(ExceI_  1, "Sheet  l!om29")  value  om29=1056082.87731268; 
excel(Excel_l, "Sheet  1  !om3")  value  om3=224248.0773 12683; 
exceI(Excel_l, "Sheet  1  lom30”)  value  om30=923345. 133333333; 
exceI(Excel_l,"Sheetl  !om33")  value  om33=634193.210646017; 
excel(Excel_l,"Sheetl!om35")  value  otn35=748560. 857378652; 
excel(Excel_l,"Sheetl!om36")  value  om36=5243 12.780065968; 
excel(Excel_l,"Sheetl!om39")  value  om39^409945. 133333333; 
excel(Excel_l,"Sheetl!om4")  value  om4=91 5 10.3333333333; 
excel(Excel_l,"Sheetl!om40")  value  oni40=3 18434.8; 
excel(Excel_  1,"  Sheet  l!om41")  value  om41=1056082.87731268; 
excel(ExceI_  1,"  Sheet  l!om42")  value  om42=923345. 133333333; 
exceI(ExceI_  1,"  Sheet  l!om47")  value  om47=748560.857378652; 
exceI(ExceI  l,"Sheetl!oni48")  value  oni48=5243 12,780065968; 
excel(Excel_l,"Sheetl  !om5")  value  om5=737648.0773 12683; 
excel(Excel_l,"Sheetl!om51")  value  om5 1=409945. 133333333; 
excel(Excel_l,"Sheetl  !oni52")  value  om52=3 18434.8; 
excel(Excel_l,"Sheetl!om53")  value  om53=737648.0773 12683; 
excel(Excel_l,"Sheetl!om54")  value  om54=6049 10.333333333; 
excel(Excel_l,"Sheetl!om55")  value  om5 5=3 15 75 8.4 106460 17; 
excel(Excel_E"Sheetl!om57")  value  om57=3 15758.41 06460 17; 
exceI(Excel_l,"Sheetl  !om58")  value  om58=183020.666666667; 
excel(Excel_l,"Sheetl  !om59")  value  oni59=430126.057378652; 
excel(Excel_l,"Sheetl  !oni6")  value  om6=604910.333333333; 
excel(Excel_l,"Sheetl!om60")  value  om60=20 5 877.9 80065968; 
excel(Excel_l,"Sheetl!oni61")  value  oni6 1=3 15758.4 106460 17; 
excel(Excel_l,"Sheetl!om63")  value  om63=103536.823760836; 
excel(Excel_l, "Sheet  1  !om64")  value  om64=  103  536.823760836; 
exceI(Excel_l, "Sheet  1  !om7")  value  om7=426314.8; 
excel(Excel_l,"Sheetl  !om9")  value  om9=4263 14.8; 
excel(Excel_l,"Sheetl  Iperfl")  value  peril  =0.94; 
exceI(Excel_l,"Sheetl  Iperfl 0")  value  peril 0=0.86; 
excel(Excel_l,"Sheetl  !peril  1")  value  peril  1=0.99; 
excel(Excel_l, "Sheet  1  Iperfl  2")  value  perfl2=0.99; 
excel(Excel_l,"  Sheet  1  !perfl3")  value  perfl3=0.9; 
excel(Excel_l,"Sheetl  Iperfl 4")  value  peril 4=0.77; 
excel(Excel_l,"Sheetl  Iperfl 5")  value  peril 5=0.99; 
excel(Excel_l,"SheetlIperfl6")  value  peril 6=0,97; 
excel(Excel_l,"Sheetl  Iperfl 7")  value  peril 7=0.99; 
excel(Excel_l,"Sheetl  Iperfl 8")  value  peril 8=0.99; 
excel(Excel_l,"SheetlIperfl9")  value  perfl9=0.9; 
excel(Excel_l,"Sheetl  Iperf2")  value  perf2=0.9; 
excel(Excel_  1, "Sheet  l!perf20")  value  perf20=0.82; 
excel(Excel_l,"Sheetl  !perf21")  value  perf21=0.94; 
excel(Excel_l,"Sheetl  Iperf22")  value  perf22=0.86; 
excel(Excel_l,"SheetlIperf23")  value  perf23=0.99; 
excel(Excel_l,"Sheetl  Iperf24")  value  perf24=0.99; 
excel(Excel_l,"SheetlIperf25")  value  perf25=0.9; 
excel(Excel_l,"SheetlIperf26")  value  perf26=0,77; 
excel(Excel_l,"  Sheet  1  Iperf27")  value  perf27=0.99; 
excel(Excel_l,"Sheetl  Iperf28")  value  perf28=0.97; 
excel(Excel_l,"Sheetl  !perf29")  value  perf29=0.99; 
excel(Excel_l,"Sheetl  Iperf3")  value  perf3=l; 
excel(Excel_l,"SheetlIperf30")  value  perf30=0.98; 
excel(Excel_l,"Sheetl  Iperf3 1")  value  perf3 1=0.88; 
excel(Excel_l,"  Sheet  1  Iperf32")  value  perf32=0.8; 
excel(Excel  l,"Sheetl  !perf33")  value  perf33=0.92; 
excel(Excel_l,"Sheetl  Iperf34")  value  perf34=0.84; 
exceI(Excel_l,"Sheetl  Iperf35")  value  perf35=0.99; 
excel(Excel_l,"Sheetl  Iperf36")  value  perf36=0.98; 
excel(Excel_l,"Sheetl  Iperf37")  value  perf37=0.88; 
excel(Excel_l,"Sheetl  Iperf38")  value  perf38=0.75; 
excel(Excel_l,"Sheetl!perf39")  value  perf39=0.98; 
excel(Excel_l,"Sheetl  Iperf4")  value  perf4=l; 
excel(Excel_l,"Sheetl  Iperf40")  value  perf40=0.95; 
excel(Excel_l,"Sheetl  Iperf41")  value  perf4 1=0.96; 
excel(Excel_l,"Sheetl  !perf42")  value  perf42=0.95; 
excel(Excel_l,"Sheetl  iperf43")  value  perf43=0,85; 
excel(Excel_l,"Sheetl  Iperf44")  value  perf44=0.77; 
excel(Excel_l,"Sheetl  Iperf45")  value  perf45=0.89; 
excel(Excel_l,"SheetlIperf46")  value  perf46=0.81; 
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exceI(Excel_l, "Sheet  1  !perf47")  value  perf47=0.96; 
exceI(ExceI_l,"Sheetl  !perf48")  value  perf48=0.95; 
excel(Excel_l,"Sheetl  !perf49")  value  perf49=0.85; 
excel(Excel_l,"Sheetl  !perf5")  value  perf5=0.99; 
excel(Excel_l,"Sheetl!perf50")  value  perf50=0.72; 
excel(Excel_l,"Sheetl  !perf5 1 ")  value  perf5 1  =0.95; 
excel(Excel_l,"Sheetl!perf52")  value  perf52=0.92; 
excel(Excel_l,"Sheetl!perf53")  value  perf53=0.99; 
excel(Excel_l,"Sheetl!perf54")  value  perf54=0.99; 
excel(Excel_l,"Sheetl!perf55")  value  perf55=0.95; 
excel(Excel_l,"Sheetl  !perf56")  value  perf56=0.87; 
excel(Excel_l,"Sheetl  !peif57")  value  perf57=0.99; 
excel(ExceI_l,"Sheetl  lperf58")  value  perf58=0.91; 
excel(Excel_l,"Sheetl  !perf59")  value  perf59=0.99; 
excel(ExceI_l, "Sheet  1  !perf6")  value  perf6=0.99; 
excel(Excel_l,"Sheetl  !perf60")  value  perf60=0.99; 
exceI(Excel_l,"Sheetl  !perf6 1 ")  value  perf6 1=0.95; 
excel(Excel_l,"Sheetl!perf62")  value  perf62=0.82; 
excel(Excel_l,"Sheetl  !perf63")  value  perf63=0.99; 
excel(Excel_l,"Sheetl!perf64")  value  perf64=0.99; 
excel(Excel_l,"Sheetl!perf65")  value  perf65=0.64; 
excel(ExceI_l,"Sheetl  !perf66")  value  perf66=0.8; 
exceI(Excel_l,"Sheetl  !perf67")  value  perf67=0.8; 
excel(Excel_l, "Sheet  1  !perf68")  value  perf68=0.75; 
excel(Excel_l,"Sheetl  !peif69")  value  perf69=0.7; 
exceI(Excel_l,"Sheetl  !perf7")  value  perf7=0.9; 
excel(Excel_l,"Sheetl  !perf70")  value  perf70=0.4; 
excel(Excel_l,"Sheetl!perf71")  value  perf71=0.5; 
excel(Excel_l,"Sheetl!perf72")  value  perf72=0.5; 
excel(Excel_l,"Sheetl  !perf73")  value  perf73=0.45; 
exceI(Excel_l,"Sheetl  !perf74")  value  perf74=0.42; 
exceI(Excel_l,"Sheetl  !perf8")  value  perf8=0.82; 
excel(Excel_l,"Sheetl  !perf9")  value  perf9=0.94; 
excel(Excel_l,"Sheetl  Itimel")  value  timel=3.5; 
excel(Excel_l,"Sheetl  ItimelO")  value  tinielO=6; 
excel(Excel_l,"Sheetl  Itimell")  value  timel  1=5.66666666666667; 
excel(Excel_l,"Sheetl  !timel2")  value  timel2=3; 
excel(ExceI_l,"Sheetl!tinie  13")  value  timel 3=5. 66666666666667; 
excel(ExceI_l,"  Sheet  1  !timel4")  value  timel4=6; 
excel(ExceI_l,"Sheetl  !timel5")  value  time  15=6. 16666666666667; 
excel(Excel_l,"Sheetl  !timel6")  value  time  16=2.66666666666667; 
exceI(Excel_l,"Sheetl  !timel7")  value  timel7=5.93060802369355; 
excel(Excel_l,"Sheetl  Itimel  8")  value  timel  8=8.93060802369355; 
excel(Excel_l,"Sheetl  !timel9")  value  timel 9=7.66666666666667; 
excel(Excel_l,"Sheetl  !time2")  value  time2=0.5; 
excel(Excel_l, "Sheet  1  !time20")  value  time20=6; 
excel(Excel_l,"Sheetl!time21")  value  time2 1=7.66666666666667; 
excel(Excel_l,"Sheetl  !time22")  value  time22=6; 
excel(Excel_l,"Sheetl  !time23")  value  time23=5.66666666666667; 
excel(Excel_l,"Sheetl  !time24")  value  time24=3; 
excel(Excel_l,"Sheetl  !time25")  value  time25=5.66666666666667; 
excel(Excel_l, "Sheet  1  !time26")  value  time26=6; 
excel(Excel_l,"Sheetl  !time27")  value  time27=6. 16666666666667; 
excel(Excel_l,"  Sheet  1  !time28")  value  time28=2. 66666666666667; 
excel(Excel_l,"  Sheet  1  !time29")  value  time29=3. 26394 13  5702688; 
excel(ExceI_l, "Sheet  1  !time3")  value  time3=3.5; 
excel(Excel_l,"Sheetl  !time30")  value  time30=6.26394135702688; 
excel(Excel_l,"Sheetl  !time3 1 ")  value  time3 1=5; 
excel(Excel_l,"Sheetl  !time32")  value  time32=6; 
excel(ExceI_l, "Sheet  1  !time33")  value  time33=5; 
excel(Excel_l, "Sheet  1  !time34")  value  time3  4=6; 
excel(Excel_l,"Sheetl  !time35")  value  time35=3; 
excel(Excel_l,"Sheetl  !time36")  value  time36=3; 
excel(Excel_l,"  Sheet  1  !time37")  value  time37=3; 
excel(Excel_l, "Sheet  1  !time38")  value  time38=6; 
excel(Excel_l,"Sheetl  !time39")  value  time39=3.5; 
excel(Excel_l,"  Sheet  1  !time4")  value  time4=6.5; 
excel(Excel_l,"Sheetl  !time40")  value  time40=2.26666666666667; 
excel(Excel_l,"Sheetl!time41")  value  time41=3.26394135702688; 
excel(Excel_l,"Sheetl  !time42")  value  time42=6.26394135702688; 
excel(Excel_l,"Sheetl  !time43")  value  time43=5; 
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excel(Excel_l,"Sheetl  !time44")  value  time44=6; 
excel(Excel_l,"Sheetl  !tune45")  value  time45=5; 
excel(Excel_l,"Sheetl  !tinie46")  value  time46=6; 
excel(Excel_l,"Sheetl  !time47”)  value  time47=3; 
excel(Excel_l,"Sheetl  !time48")  value  tiine48=3; 
excel(Excel_l,"Sheetl  !time49")  value  time49=3; 
excel(Excel_l,"  Sheet  1  !time5")  value  time5=5.93060802369355; 
exceI(Excel_l,"Sheetl  ItimeSO”)  value  time50=6; 
excel(Excel_l,"Sheetl  !tmie5 1 ")  value  time5 1=3.5; 
excel(Excel_l, "Sheet  1  !tinie52")  value  time52=2. 16666666666667; 
excel(Excel_l, "Sheet  1  !time53")  value  time53=4.56438 171941392; 
excel(ExceI_l,"Sheetl  !time54")  value  time54=7.56438171941392; 
excel(Excel_l,"Sheetl  !time55")  value  time55=6.30044036238704; 
excel(Excel_l, "Sheet  1  !time56")  value  time56=7.30044036238704; 
excel(Excel_l, "Sheet  1  !time57")  value  time57=6.30044036238704; 
excel(ExceI_l,"Sheetl  !tmie58")  value  time58=7.30044036238704; 
excel(Excel_l,"Sheetl  !time59")  value  time59=4.30044036238704; 
excel(Excel_l,"Sheetl  !tinie6")  value  time6=8.93060802369355; 
excel(Excel_l,"Sheetl  !time60")  value  time60=4.30044036238704; 
excel(Excel_l,"Sheetl  !time6r')  value  time6 1=4.3004403623  8704; 
excel(Excel_l,"Sheetl  !time62")  value  tinie62=7.30044036238704; 
excel(Excel_l, "Sheet  1  !time63")  value  time63=4.80044036238704; 
exceI(Excel_l,"Sheetl  !tmie64")  value  time64=l. 80044036238704; 
excel(Excel_l, "Sheet  1  !time65")  value  time65=3; 
excel(Excel_l, "Sheet  1  !time66")  value  time66=5. 66666666666667; 
excel(Excel_l,"Sheetl  !time67")  value  time67=5. 66666666666667; 
excel(Excel_l,"Sheetl  !time68")  value  time68=3; 
excel(Excel_l, "Sheet  1  !time69")  value  time69=3; 
exceI(Excel_l,"Sheetl  !time7")  value  time7=7.66666666666667; 
excel(Excel_l,"Sheetl  !time70")  value  time70=6; 
excel(Excel_l,"Sheetl  !time71")  value  time7 1=8.66666666666667; 
excel(Excel_l,"Sheetl  !time72")  value  time72=8.66666666666667; 
excel(Excel_l,"Sheetl  !tinie73")  value  time73=6; 
excel(Excel_l,"Sheetl  !tinie74")  value  tinie74=6; 
excel(Excel_l,"Sheetl  !tiine8")  value  time8=6; 
excel(ExceI_l,"Sheetl  !time9")  value  time9=7.66666666666667; 
excel(Excel_l,"Sheetl  lutill")  value  utill; 
exceI(ExceI_l,"Sheetl  lutillO")  value  utillO; 
excel(Excel_l,"Sheetl  lutill  1")  value  utilll; 
excel(Excel_l,"Sheetl  !utill2")  value  utill2; 
excel(Excel_l,"Sheetl!utill3")  value  utiil3; 
excel(Excel_l,"Sheetl  !utill4")  value  utilM; 
excel(Excel_l,"Sheetl  !utill5")  value  utill5; 
excel(ExceI_l,"Sheetl!utill6")  value  utill6; 
excel(Excel_l,"Sheetl  !utill7")  value  utill7; 
excel(ExceI_l,"Sheetl  !utill8")  value  utill8; 
excel(Excel_l,"Sheetl  !utill9")  value  utill9; 
excel(Excel_l,"Sheetl  !utiI2")  value  util2; 
excel(Excel_l,"Sheetl  !util20")  value  util20; 
excel(Excel_l,"Sheetl!util21")  value  util21; 
excel(Excel_l,"  Sheet  1  !util22")  value  util22; 
excel(Excel_i;'Sheetl  !util23")  value  util23; 
excel(Excel_l,"Sheetl  !util24")  value  util24; 
excel(Excel_l,"Sheetl  !util25")  value  util25; 
excel(Excel_l,"Sheetl  !util26")  value  util26; 
excel(Excel_l, "Sheet  1  !util27")  value  util27; 
excel(Excel_l,"Sheetl  !util28")  value  util28; 
excel(Excel_l,"Sheetl!util29")  value  util29; 
excel(Excel_l,"Sheetl  !util3")  value  util3; 
exceI(Excel_l,"Sheetl  !utiI30")  value  utiBO; 
excel(ExceLl,"Sheetl  !util3 1 ")  value  util3 1 ; 
excel(Excel_l,"Sheetl  !util32")  value  util32; 
excel(Excel_l,"Sheetl  !util33")  value  util33; 
excel(Excel_l,"Sheetl!util34")  value  util34; 
excel(Excel_l,"Sheetl  !util35")  value  util35; 
excel(Excel_l,"  Sheet  1  !util36")  value  util36; 
excel(Excel_l,"Sheetl  !util37")  value  util37; 
excel(Excel_l,"Sheetl  !util38")  value  util38; 
excel(Excel_l,"Sheetl  !util39")  value  util39; 
excel(Excel_l,"Sheetl  !util4")  value  util4; 
excel(Excel_l,"Sheetl  !util40")  value  util40; 
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excel(ExceI_l, "Sheet  1  lutiWl")  value  util41; 
excel(Excel_l,"Sheetl  !util42")  value  util42; 
excel(ExceI_l, "Sheet  1  !util43”)  value  util43; 
excel(Excel_l,"Sheetl  !util44")  value  util44; 
excel(Excel_l, "Sheet  1  !util45")  value  utiI45; 
exceI(Excel_l,"Sheetl  !util46")  value  utiI46; 
excel(Excel_l,"Sheetl  !util47")  value  util47; 
excel(Excel_l,"Sheetl  !util48")  value  util48; 
exceI(Excel_l,"Sheetl!util49")  value  utiI49; 
excel(Excel_l,"Sheetl  !util5")  value  util5; 
exceI(Excel_l,"Sheetl  !util50")  value  util50; 
excel(Excel_l, "Sheet  1  !util5 1 ")  value  util5 1 ; 
excel(Excel_l,"Sheetl!util52”)  value  util52; 
excel(Excel_l,"Sheetl!util53")  value  util53: 
excel(Excel_l,"  Sheet  1  !util54”)  value  util54; 
excel(Excel_i;'Sheetl!util55")  value  util55; 
excel(ExceI_l,"Sheetl!util56")  value  util56; 
excel(ExceI_l,"Sheetl  !util57”)  value  util57; 
excel(Excel_l,"Sheetl  !util58")  value  util58; 
excel(Excel_l,"  Sheet  1  lutil59")  value  util59; 
excel(Excel_l, "Sheet  1  !util6")  value  utiI6; 
exceI(Excel_l,"Sheetl  !util60")  value  util60; 
excel(Excel_l,"Sheetl  !util6 1 ")  value  util6 1 ; 
excel(Excel_l,"Sheetl!util62")  value  util62; 
excel(Excel_l,"Sheetl  !util63")  value  util63; 
excel(Excel_l,"Sheetl  !util64")  value  util64; 
excel(Excel_l,"Sheetl  !util65")  value  util65; 
excel(Excel_l,"Sheetl!util66”)  value  util66; 
excel(Excel_l,"Sheetl  !util67”)  value  util67; 
excel(Excel_l,"Sheetl  !util68")  value  util68; 
excel(Excel_l,"Sheetl!util69")  value  util69; 
excel(Excel_l, "Sheet  1  !util7")  value  util7; 
excel(Excel_l,"Sheetl  lutil70")  value  util70; 
exceI(Excel_l,"Sheetl  !util71")  value  util71; 
excel(Excel_l,"Sheetl  !utiI72")  value  util72; 
excel(Excell,"Sheetl  !util73")  value  util73; 
excel(Excel_l,"Sheetl  !util74")  value  util74; 
excel(Excel_l,"Sheetl  !util8")  value  util8; 
excel(Excel_l,"Sheetl  !util9")  value  util9; 

decision  Which_Strategy.  (USA,US  A,SA _ ^U} ; 

decision  US  A_Priniary_TechnoIogy_2_choices.  {DUS,GTA}; 
decision  US  A_Secondary_Technology_2_choices.  {Bio,Pump_treat} ; 
decision  US_Primary_Technology_2_choices.  {two_phase,Lasagna} ; 
decision  US_Secon^_Technology_4_choices.  {RF, six jphase, hydro, none}; 
decision 

Aquifer_technology_12_choices.{oxid _ bio,oxid _ ^none, cosolvent _ ^bio,cosolvent___none,surf _ bio,surf _ ^none,UVB. 

B _ none.pump  bio, pump _ ^none.PTZ _ bio,PTZ _ none}; 

decision  UnsatPrimaryT  echnology_2_choices.  {  S  VE,Two_phase  } ; 
decision  Unsat_Secondary_Technology_4_choices.{RF,six j)hase, hydro, none}; 

decision  SA_Technology_2_choices.  (Pumpjreat _ ^None,Pump_treat _ Bio} ; 

excel(Excel_l,"Sheetl!costlO")  value  costlO; 
excel(Excel_l,”Sheetl  !costl4")  value  costl4; 
excel(Excel_l,"  Sheet  1  !cost20”)  value  cost20; 
excel(Excel_  1,"  Sheet  l!cost22")  value  cost22; 
excel(ExceI_l,"Sheetl!cost26")  value  cost26; 
excel(ExceI_l, "Sheet  1  !cost3 1")  value  cost31; 
excel(ExceI_l,"Sheetl  !cost34")  value  cost34; 
excel(Excel_l,"Sheetl  !cost37")  value  cost37; 
excel(ExceI_l,"Sheetl  !cost38")  value  cost38; 
exceI(Excel_l, "Sheet  1  !cost43")  value  cost43; 
excel(Excel_l,"Sheetl  !cost44")  value  cost44; 
excel(Excel_l,"Sheetl  !cost45")  value  cost45; 
excel(Excel_l,"Sheetl  !cost46")  value  cost46; 
excel(Excel_l,"Sheetl!cost49")  value  cost49; 
excel(Excel_l,"Sheetl  !cost50")  value  costSO; 
excel(Excel_l,"Sheetl!cost56")  value  cost56; 
excel(Excel_l,"Sheetl  !cost62")  value  cost62; 
excel(Excel_l, "Sheet  1  !cost65")  value  cost65; 
excel(ExceI_l,"Sheetl  !cost66")  value  cost66; 
exceI(Excel_l,"Sheetl  !cost67")  value  cost67; 
excel(Excel_l,"Sheetl!cost68")  value  cost68; 


_bio,UV 
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excel(Excel_l,"Sheetl  !cost69")  value  cost69; 
excel(Excel_l,"Sheetl  !cost70")  value  cost70; 
excel(Excel_l,"Sheetl  !cost71")  value  cost71; 
excel(Excel_l,'’Sheetl  !cost72'‘)  value  cost72; 
excel(Excel_l,"Sheetl!cost73")  value  cost73; 
excel(Excel_l,"Sheetl  !cost74‘')  value  cost74; 
excel(Excel_l,"Sheetl  !cost8”)  value  cost8; 
excel(Excel_l,"Sheetl  limport  cost  vvt")  value  import_cost_wt; 
excel(Excel_l,"Sheetl  !iniport_perf_wt")  value  import_perf_wt; 
excel(Excel_l,"Sheetl!import_time_wt")  value  import  time  wt; 
excel(Excel_l, "Sheet  1  !omlO")  value  omlO; 
excel(ExceI_l,"Sheetl!oml4")  value  oml4; 
excel(Excel_l, "  Sheet  l!om20")  value  om20; 
excel(Excel_l,"Sheetl  !om22”)  value  om22; 
exceI(ExceI_l,"Sheetl!om26")  value  om26; 
excel(ExceI_l,"Sheetl  !om3 1")  value  om3 1; 
excel(ExceI_l,"Sheetl  !om32")  value  om32; 
excel(ExceI_l, "Sheet  1  !om34")  value  om34; 
excel(ExceI_l, "Sheet  1  !om37")  value  om37; 
excel(ExceI_l,"Sheetl  !oni38")  value  om38; 
exceI(Excel_l, "Sheet  1  !om43")  value  om43; 
excel(Excel_l,"Sheetl  !om44")  value  om44; 
exceI(Excel_l,"Sheetl  !om45")  value  om45; 
excel(Excel_l,"Sheetl  !om46")  value  om46; 
excel(Excel_l,"Sheetl!om49")  value  om49; 
excel(Excel_l,"Sheetl!om50")  value  om50; 
excel(Excel_l,"Sheetl  !om56")  value  om56; 
excel(Excel_l,"Sheetl!om62")  value  om62; 
excel(Excel_l,"Sheetl  !om65")  value  om65; 
excel(Excel_l,"Sheetl  !oni66")  value  om66; 
excel(Excel_l,"Sheetl  !om67")  value  om67; 
excel(Excel_  1,"  Sheet  l!om68")  value  om68; 
excel(Excel_l,"Sheetl!om69")  value  om69; 
exceI(Excel_l, "Sheet  1  !om70")  value  om70; 
excel(Excel_l,"Sheetl  !om7r')  value  om71; 
excel(Excell,"Sheetl  !om72")  value  oni72; 
excel(Excel_l, "Sheet  1  !om73")  value  om73; 
exceI(Excel_l,"Sheetl  !om74")  value  om74; 
excel(Excel_l,"Sheetl!om8")  value  om8; 

sequence(  attributes  =  5, 

objective  =  $5  ): 


decide 

to  WhichStrategy.USAthen  decide 

to  USA_Primary_Technology_2_choices.DUS  then  decide 

to  USA_Secondary_Technology_2_choices.Bio  and  get  cost KtimeLperfl,oml, util  1 
to  USA_Secondary_Technology_2_choices.Pump_treat  and  get  cost2,time2,perf2,om2,util2 
to  US A_Primary_Technology_2_choices.GTA  then  decide 

to  US  A_Secondary_Technology_2_choices.  Bio  and  get  cost3,time3  ,perf3,om3,util3 

to  US  A_Secondary_Technology_2_choices.Pump_treat  and  get  cost4,time4,perf4,om4,util4 

to  Which_Strategy.US _ A  then  decide 

to  US_Priniary_Technology_2_choices.two_phase  then  decide 

to  US  Secondary_Technology_4_choices,RF  then  decide 

to  Aquifer_technology_12_choices.oxid _ ^bio  and  get  cost5,time5,perf5,om5,utiI5 

to  Aquifer_technology_12_choices.oxid _ none  and  get  cost6,time6,perf6,om6,util6 

to  Aquifer_technology_12_choices.cosolvent__bio  and  get  cost7,time7,perf7,om7,util7 

to  Aquifer_technology_12_choices. cosolvent _ none  and  get  cost8,tinie8,perf8,om8,util8 

to  Aquifer_technology_12_choices.surf _ ^bio  and  get  cost9,tinie9,perf9,oni9,utiI9 

to  Aquifer_technology_12_choices.surf _ ^none  and  get  costlO,timelO,perflO,omlO,utillO 

to  Aquifer_technology_12_choices.UVB _ ^bio  and  get  costl  Utimel  l,perfl  l,oml  l,utill  1 

to  Aquifer  technologv  12  choices.UVB  none  and  get  costl2,timel2,peifl2,oml2,utill2 
to  Aquifer_technology_12_choices.pump__bio  and  get  cost  1 3, time  1 3, peril 3, om  13, util  13 
to  Aquifer_technology_12  choices.pump _ ^none  and  get 

cost  1 4,time  1 4,peril  4,om  1 4,util  1 4 

to  Aquifer_technology_12_choices.PTZ _ bio  and  get  costl5,timel5,peril5,oml5,utill5 

to  Aquifer Jechnology_12_choices.PTZ _ ^none  and  get  costl6,timel6,peril6,oml6,utill6 

to  US_Secondary_Technology_4_choices.six_phase  then  decide 

to  Aquifer_technology_12_choices.oxid _ bio  and  get  costl7,timel7,perfl7,oml7,utiI17 

to  Aquifer  technology  12  choices.oxid  none  and  get  costl8,timel8,perfl8,oml8,utill8 
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cost  1 9,time  1 9,perfl  9,om  1 9, util  1 9 
cost20,time20,perf20,om20,utiI20 


to  Aquifer_technology_12_choices. cosolvent _ bio  and  get 

to  Aqtiifer  technology  l 2_choices. cosolvent  none  and  get 

to  Aquifer_technology_12__choices.surf _ bio  and  get  cost21,tinie21,perf21,om21,util21 

to  Aquifer  technology  1 2  choices.surf  none  and  get  cost22,tinie22.perf22,oni22,util22 

to  AquiferJeclinology_12_choices.UVB _ ^bio  and  get  cost23,time23,perf23,om23,util23 

to  Aquifer  tcchnology_  1 2  choices.UVB  none  and  get  cost24,tinie24,petf24,oni24,util24 

to  Aquifer_technology_12_choices.pump _ bio  and  get  cost25,time25,peTf25,om25,util25 

to  Aquifer  Jeclinologv  12  choices.pump  none  and  get 


cost26,tittie26,perf26,om26,util26  _ . 

to  Aquifer_technology_12_choices.PTZ _ bio  and  get  cost274inie27,perf27,om27,util27 

to  Aquifer  technology  12  choices. PTZ  none  and  get  cost28,tinie28,perG8,oni28,util28 

to  US_Secondary_Technology_4_choices.hydro  then  decide 

to  Aquifer_technoIogy_12_choices.oxid _ bio  and  get  cost294ime29,perf29,om29,util29 

to  Aquifer_teclinology_12_choices.oxid _ none  and  get  cost30,tinie30,perf30,om30,utiI30 

to  Aquifer  jechnology_12_choices. cosolvent _ bio  and  get 

cost3  l,time3  l,peTf3  l,om3  l,util3 1 

to  Aquifer_technology_  1 2  choices.cosol vent  none  and  get 


cost32,time32,perfi2,om32,util32 

to  Aquifer_technology_12_choices.surf _ bio  and  get  cost33,tinie33,perf33,om33,util33 

to  Aqiiifer  technology  12  choices.surf  none  and  get  cost34,time34,perf34,om34,util34 

to  Aquifer_technology_12_choices.UVB _ bio  and  get  cost35,time35,perf35,om35,util35 

to  Aquifer_technology_12_choices.UVB _ ^none  and  get  cost36,tinie36,perf36,om36,util36 

to  Aquifer_teclmology_12_choices.pump _ ^bio  and  get  cost37,tiiue37,perf37,oni37,util37 

to  Aqiiifer  technology  1 2  choices.pump  none  and  get 


cost38,tinie38.perf38,om38,util38  , 

to  Aquifer_technology_12_choices.PTZ _ bio  and  get  cost39,tmie39,perf39,om39,util39 

to  Aquifer_technology_12_choices.PTZ _ none  and  get  cost40,time40,perf40,om40,util40 

to  US_Secondary_Technology_4_choices.none  then  decide 

to  Aquifer_technology_12_choices.oxid _ ^bio  and  get  cost41,tinie41,perf41,om41,util41 

to  Aquifer_technology_12  choices.oxid__none  and  get  cost42,tinie42.perf42,om42,util42 
to  Aquifer_technology_l 2  choices.cosolvent  bio  and  get 


cost43,tinie43,perf43,om43,util43 

cost44,time44,perf44,om44,util44 


to  Aquifer_technology_12_choices.cosolvent _ ^none  and  get 

to  Aquifer_technology_12_choices.surf _ ^bio  and  get  cost45,tinie45,perf45,om45,util45 

to  Aquifer_technology_12_choices.surf _ ^none  and  get  cost46.tinie46,perf46,om46,util46 

to  AquiferJechnology_12_choices.UVB _ ^bio  and  get  cost47,time4Zperf47,om47,util47 

to  Aquifer  technology  !  2_choices.UVB _ ^none  and  get  cost48,tiiue48jperf48,om48,util48 

to  Aquifer_technology_12_choices.pump _ ^bio  and  get  cost49,tinie49,perf49,om49,util49 

to  Aquifer_technology_12_choices.pump _ ^none  and  get 


cost50,tinie50,perf50,om50,util50 

to  Aquifer_technoIogy_12dioices.PTZ _ ^bio  and  get  cost5  l,time5 l,perD Lom5 l,util5 1 

to  Aquifer  technology  12_choices.PTZ  none  and  get  cost52,tinie52,perf52,oni52,util52 
to  US_Primary_Technology_2_choices.Lasagna  then  decide 

to  US_Secondary_Technology_4_choices.none  then  decide 

to  Aquifer_technology_12_choices.oxid _ bio  and  get  cost53,time53,perf53,om53,util53 

to  Aquifer  Jechnology_12_choices.oxid _ none  and  get  cost544inie54,perf54,om54,util54 

to  Aqiiifer  technology  1 2  choices.cosolvent  bio  and  get 


cost55,time5  5,perf5  5,om55,util55 
cost56,time56,perf56,om56,util56 


to  Aquifer_technology_12_choices,cosolvent _ none  and  get 

to  Aquifer_technology_12_choices.surf _ bio  and  get  cost57,time57,perf57,om57,util57 

to  Aquifer  technology  1 2  choices.surf  none  and  get  cost58,tinie58,perf58,oni58,utiI58 
to  Aquifer JechnologyIl2_choices.UVB_bio  and  get  cost59,time59.per£59,om59,utiI59 
to  Aquifer_technology_12_choices.UVB___none  and  get  cost60,time60,perf60,om60,util60 

to  Aquifer_technology_12_choices.pump _ bio  and  get  cost6l4ime61,perf61,om61,util61 

to  Aquifer  technologv  12  choices.pump  none  and  get 


cost62,time62,perf62,om62,util62 

to  Aquifer  technologv  12  choices.PTZ  bio  and  get  cost63,tinie63,perf63,om63,util63 

to  Aquifer_technology_12_choices.PTZ _ none  and  get  cost64,time64,perf(54,om64,util64 

to  US_Secondary_Technology_4_choices.{RF, six jDhase, hydro}  and  get  99999999,999,0,0,0 

to  Wliich  Strategy.S  A _ then  decide 

to  S  A_Technology_2_choices.Pump_treat _ ^None  then  decide 

to  Unsat  Priniary_T®c^°l®Sy_2_ch^*c®s.SVE  and  get  cost70,tinie70,perf70,om70,util70 
to  Unsat_Primary_Technology_2_choices.Two_phase  then  decide 

to  Unsat  Secondary_Technology_4_choices.RF  and  get  cost71,tinie71,perf71,om71,util71 
to  Unsat_Secondary_TechnoIogy_4_choices.six_phase  and  get 


cost72,tinie72,perf72,om72,util72 
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to  Unsat_Secondary_Technology_4_choices.hydro  and  get 

cost73,tinie73  ,perf73  ,om73  ,util73 

to  Unsat_Secondary_TechnoIogy_4_choices.none  and  get 

cost74,tinie74,perf74,om74,util74 

to  S  A_Technology_2_choices.Pump_treat _ Bio  then  decide 

to  Unsat_Primary_Technology_2_choices.SVE  and  get  cost65,time65,perf65,om65,util65 
to  Unsat_Primary_Technology_2_choices.Two  jjhase  then  decide 

to  Unsat_Secondary_Technology_4_choices.RF  and  get  cost66,time66,perf66,om66,util66 
to  Unsat_Secondary_Technology_4_choices.six_phase  and  get 

cost67,tuiie67,perf67,om67,util  6  7 

to  Unsat_Secondary_Technology_4_choices.hydro  and  get 

cost6  8,tiine6  8,perf6  8,om68,utiI6  8 

to  Unsat_Secondary_Technology  4_choices.none  and  get 

cost69,time69,perf69,om69,utiI69 
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Phase  2  Model: 


string  Excel  l ="c:\\winword\\tony\\thesis\\test202.xls"; 
excel(Excel_l,"Sheetl  !export_cost_wt")  value  export_cost_wt=0.44; 
excel(Excel_l,"Sheetl !export_perf_w1”)  value  export j)erf_wt-0.44; 
excel(Excel_l/'Sheetl!export_time_wt’')  value  exportjime_wt=0,12; 
excel(Excel_i;'Sheetl  !import_cost_\vt")  value  import_cost_wt; 
excel(Excel_l,"Sheetl!import_perf_wt")  value  import_perf_wt; 
exceI(Excel_l, "Sheet  1  limport  jimejwt")  value  import  Jime  wt; 
value  cost_utility_factor=3 0000000; 
value  perfoniiance_utility_factor=.9; 
value  Performance_Desired=.95; 
value  Interest_Rate=.025; 

value  Perfonnance_Utility=@if(Performance_Desired>performance_utility__factor,(Performance_Desired- 


performance_utility_factor)*  1 0,0); 

decision  Wluch_Train.{no_3,no_4,no_94,no_l  19,no_135,no_148,no_156}; 
chance  TimeRequired.  {Low,Med,High}l  Which_Tram= 


{.185,.63,.185}, 
{,185, .63, .185}, 
{.185,.63,185}, 
{.185, .63, .185}, 
{.185,.63,.185}, 
{.185, .63, .185}, 
{.185, .63, .185}, 


//  Wliich_Train.no _ 3 

//  Which_Train.no _ 4 

//  Which_Train.no _ ^94 

//Which_Train.no_119 

//  Which_Tram.no _ 135 

//  Which_Train.no__148 
//Which  Train.no _ 156 


//  Which_Train.no 

@if(Performance_Desired<=.99,-7.93+ 1 0. 52*Performance_Desired,  1000),  // 

@if(Performance_Desired<=.99,- 1 5.86+2 1 .03*Performance_Desired,  1 000),  // 
(^ifTPerformance  Desired<=.99,-29.15+39.39*Performance_Desired,1000),  // 

//  Which_Train.no 

@if(Performance_Desired<=.9,8 1 .22*pow(Performance_Desired, 22.45),  1000),  // 
@if(Performance_Desired<=.9,l  82.77*pow(Performance_Desired,23.03),  1000), 
^ifTPerformance  Desired<=.9,287.93*pow(Performance_Desired,2 1 . 12),  1000), 

//  Which_Train.no_ 

@if(Performance_Desired<=.9,.  16*@exp(Performance_Desired*4.54),  1000),  // 

@if(Perfomiance_Desired<=.9,.33*@exp(Performance_Desired*4.54),1000),  // 
(^ifrPerfonnance_Desired<=.9,.66*@exp(Peiformance_Desired*4.54),1000),  // 

//  Which_Train.no 

@il(Perfbrmance_Desired<=.99,2. 1 3+7. 1  *@ln(Perfonnance_Desired),  1000),  // 

@if(Performance_Desired<=.99,- 1 1 .64+ 1 5.77*Performance_Desired,  1 000),  // 
(SlifTPerformance  Desired<=.99,8.61*pow(Performance_Desired,5.47),1000),  // 

//  Which_Train.no 

@if(Peiformance_Desired<=.95,-2 .75+5 *Performance_Desired,  1 000),  // 

@if(Performance_Desired<=.95,.14*@exp(Performance_Desired*3.47),1000),  // 
(SifTPerformance  Desired<=.95,9.87*pow(Performance_Desired,4.15),1000),  // 

//  Which_Train.no 

@ifi;Performance_Desired<=.99,63.05*pow(Performance_Desired,  1 .99),  1000),  // 
@if(Performance_Desired<=.99,77.16+158.78*@ln(Performance_Desired),1000), 
(oJifTPerformance  Desired<=.99,2.03*@exp(Performance_Desired*3.91),1000), 

//  Which_Train.no 


@if(Performance_Desired<=. 99, 1 .73  *@exp(Performance_Desired*3 . 73),  1 000), 

@if(Performance_Desired<==.99,77.16+158.78*@ln(Performance_Desired),1000), 

@if(Performance_Desired<=.99,2.03*@exp(Perfbrmance_Desired*3.91),1000); 


value  Capital_Cost]Wliich_Train- 
243206, 

2134837, 

2470812, 

1193059, 

323972, 

2019759, 

988418; 

value  0_M_Cost|Which_Train,Time_Required= 


//  Which_Train.no _ ^3 

//  Which_Train.no _ 4 

//  Which_Train.no _ 94 

//Which_Train.no_119 

1 1  Which_Train.no _ 135 

1 1  Which_Train.no _ 148 

//  Which_Train.no _ 156 

//  Which_Train.no 
TimeRequired.  Low 
TimeRequired.Med 
TimeRequired.  High 

//  Which_Train.no 
TimeRequired.  Low 
TimeRequired.  Med 


@pv(336558,Interest_Rate,Time_Required),  // 
@pv(336 5 5  8,Interest_Rate,T ime  Required),  // 
@pv(336558,Interest_Rate,Time_RequiredX  // 

@pv(203 820,Interest_Rate,Time_Required),  // 
@pv(203  820,Interest_Rate,Time_Required),  // 


3 

TimeRequired.  Low 
Time_Required.Med 
Time_Required.  High 
A 

Time_Required.  Low 
//  TimeRequired.Med 

//  Time_Required.High 

94 

T  imeRequired.  Low 
Time_Required.Med 
T  ime_Required.  High 
119 

Time_Required.Low 
T  imeRequired.  Med 
Time_Required.High 
135 

Time_Required.Low 

TimeRequired.Med 

TimeRequired.High 

_148 

T  ime_Required.  Low 
//  T  imeRequired.Med 

//  Time_Required.High 

_156 

//  Time_Required.Low 

//  Time_Required.Med 

//  TimeRequired.High 


H-9 


@pv(203  820,Interest_Rate,Time_Reqiiired),  // 

@pv(62 57 1 3 ,Interest_Rate  Jime_Required),  // 
@pv(6257 1 3 ,Interest_Rate,Time_Required),  // 
@pv(6257 1 3 ,Interest_Rate,Time_Required),  // 

@pv( 1 0 56083,Interest_Rate,Time_Required), 
@pv(  1 0 56083,Interest_Rate,Time_Required), 
@pv(  1 0 56083,Interest_Rate,Time_Required), 

@pv(737648,Interest_RateJime_RequiredX  // 

@pv(737648Jnterest_RateTime_Required),  // 
@pv(737648,Interest_Rate,Time_Required),  // 

@pv(  1 03 537,Interest_Rate,Time_Required),  // 
@pv(  103537,Interest_Rate,Time_Required),  // 
@pv(  1 03  537,Interest_Rate,Time_Required),  // 


TimeRequired.  High 

//Which_Tram.no_94 
Time  Required.Low 
Time_Required.Med 
TimeRequired.HigJi 

//  Which_Train.no _ 1 19 

//  Time_Required.Lo\v 

//  TimeRequired.Med 

//  TimeRequired.High 

//  Which_Train.no _ 135 

TimeRequired.Low 

Time_Required.Med 

TimeRequired.HigJi 

//  Which_Train.no _ 148 

TimeReqiiired.  Low 
TimeRequired.  Med 
TimeReqiiired.  High 

//  Which_Train.no _ 1 56 

TimeRequired.Low 
Time_Required.Med 
TimeRequired.  High 


@pv(3 1 6627,Interest_RateJime_Required),  // 

@pv(3 16627,Interest_RateJime_RequiredX  // 

@pv(316627,Interest_Rate,Time_Required);  // 
value  NPV_Cost=0_M_Cost+Capital_Cost; 

value  Cost_Utility=@if(NPV_Cost<30000000,l-(l/cost_utility_factor)*NPV_Cost,.7-.00000002*NPV_Cost); 
value  Time  Utility|Time_Required= 

@if(Time_Required<6,l,@ifi:Time_Required>18,0,L5-.083333*Time_Required)),  // Time_Required.Low 

@il(Time_Required<6.1,@ifl;Time_Required>18.0,1.5-.083333*Time_Required)),  // Time_Required.Med 

@if!:Time_Required<6.1,@if(Time_Required>18,0,1.5-.083333*Time_Required));  // Time_Required.High 

value  Overall_Utility=import_cost_wt*Cost_Utility+import_time_wt*Time_Utility+import_perf_wt*Perforraance_Utility; 


sequence(  attributes  =  9, 

objective  =  $5  ): 


set  Which_Train.no _ 135  and  get  Capital_Cost,Performance_Desired,Performance_Utility,0, 0,0, 0,0,0  then 

gamble  onTime_Required and  get  0,0,0,O_M_Cost,NPV_Cost,Time_Required,Cost_Utility,Time_Utility,OveralI_Utility 
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Appendix  I  Performance  Data 

Table  1. 1  contains  data  developed  on  23  Jan  97  during  a  trip  to  Paducah 


80% 

90% 

99% 

DUS  for  half  -  Bio  for  the  rest 

Train 

Low 

Med 

High 

Low 

High 

Low 

Med 

High 

3 

0.5 

1.0 

2.5 

1.5 

3.0 

6.0 

2.5 

5.0 

10.0 

series 

80% 

85% 

90% 

DUS  0.5  to  .9,  Pump  and  treat 

Train 

Low 

Med 

High 

Low 

Med 

High 

Low 

Med 

High 

doesn't  help 

4 

0.5 

1.0 

2.5 

2.5 

5.0 

10.0 

7.0 

15.0 

30.0 

series 

75% 

85% 

90% 

UVB  doesn't  go  quick  for 

Train 

Low 

Med 

High 

Low 

Med 

High 

Low 

Med 

High 

3rd  zone  -  good  for  first  two 

94 

5.0 

10.0 

20.0 

7.5 

15.0 

30.0 

10.0 

20.0 

40.0 

parallel 

80% 

90% 

99% 

Oxidation  followed  by  bio  is 

Train 

Low 

Med 

High 

Low 

Med 

High 

Low 

Med 

High 

critical 

119 

0.5 

1.0 

2.5 

1.5 

2.5 

5.0 

2.0 

4.0 

8.0 

parallel 

75% 

85% 

95% 

Bugs  must  wait  for  oxidation 

Train 

Low 

Med 

High 

Low 

Med 

High 

Low 

Med 

High 

to  clear 

135 

1.0 

2.0 

3.0 

1.5 

2.5 

5.0 

2.0 

4.0 

8.0 

Series 

90% 

95% 

99% 

Barrier  must  wait  for  source  to 

Train 

Low 

Med 

High 

Low 

Med 

High 

Low 

Med 

High 

dissolve 

148 

50.0 

60.0 

70.0 

60.0 

70.0 

80.0 

60.0 

75.0 

100.0 

Series 

90% 

95% 

99% 

Train 

Low 

Med 

High 

Low 

Med 

High 

Low 

Med 

High 

156 

50.0 

60.0 

70.0 

60.0 

70.0 

80.0 

75.0 

100.0 

Parallel 

The  data  in  Table  1. 1  was  developed  through  a  session  vrith  technical  experts  at  the  DOE. 
The  high,  low,  and  medium  values  were  determined  through  this  panel  of  experts,  then  the 
data  was  reviewed  independently  by  another  group  of  experts. 
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Plots  of  the  raw  data  follow: 


-  Medium 
-High 


Low 

Medium 


Perf  Perf 


The  data  points  shown  above  were  fit  using  the  following  regression  equations: 


Exponential  equations  are  of  the  form: 
Power  equations  are  of  the  form: 
Linear  equations  are  of  the  form: 
Logarithmic  equations  are  of  the  form: 


y  =  ax’' 

y  =  a  +  bx 
y  =  a+  b(\rcc) 


An  +  bIXi  =  27, 

AIKi  +  =  2(TO 


where: 


Regression 

A 

Xi 

Yi 

Exponential 

ln(a) 

Xi 

Infyi) 

Power 

ln(a) 

ln(yi) 

Linear 

a 

Xi 

Yi 

Logarithmic 

a 

InfXj) 

Yi  _ 

T rain  3 

In  terms  of  time  as  a  function  of  performance: 
time  =  -7.93+10.52*perf  low  f(x) 

time  =  -15.86+21.03*perf  med  g(x) 

tme  = -29.15+39.39*perf  high  h(x) 
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Train  119 


In  terms  of  time  as  a  function  of  performance: 
time  =  2.13+7.10*ln(perf)  low  f(t) 

time  =  -11.64+15.77*perf  med  g(t) 

tme  =  8.61*perf®''’  high  h(t) 


0.8  0.85  0.9  0.95 

performance 


Train  135 

In  terms  of  time  as  a  fiinction  of  performance: 
time  =  -2.75+5*perf  low  f(t) 

time  =  0.14*0^'”"'’''^  med  g(t) 

tme  =  9.87*perf  high  h(t) 
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Train  148 

In  terms  of  time  as  a  function  of  performance: 
time  =  63.05 *perf^'^  low  f(y) 

time  =  77.16+158.78*ln(perf)  med  g(y) 

tme  =  2.03*e^®’’*’®'^  high  h(y) 


performance 


Train  156 

In  terms  of  time  as  a  function  of  performance: 
time  =  I.73*e^''^*P''^  low  f(y) 

time  =  77.16+158.78*ln(perf)  med  g(y) 

tme  =  2.03*e^®'’P''^  high  h(y) 


time 


performance 
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Appendix  J  Spill  Volume  Information 


Because  of  the  uncertain  nature  of  estimating  the  volume  of  a  contaminant  spill, 
many  approximations  have  been  created.  Shown  below  are  two  such  approximations. 
The  first  one,  Figure  J.  1  shows  the  initial  estimate  of  the  spill  volume.  Figure  J.2  shows 
another  estimate.  Because  of  the  uncertainty  surrounding  this  issue,  and  because  of  the 
lack  of  availability  of  scaleable  data,  neither  of  the  two  estimates  were  used.  Instead,  a 
point  estimate  of  7,500  gallons  was  used  to  estimate  the  cost,  time,  and  performance 
values.  It  is  these  values  (the  ones  based  on  7,500  gallons)  that  were  used  as  a  data  basis 
for  the  analysis  contained  in  this  thesis  effort. 


Figure  J.l  Initial  Estimate  of  Spill  Volume 


J-1 


Cumulative  Probability 


Cum.  Distribution  Function  of  Spill  Quantity 


Figure  J.2  Another  Estimate  of  the  Spill  Volume  -  Acquired 
During  23  Jan  97  Trip  to  Paducah 
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